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Executive Summary 

The performance of an axial-flow fan-compressor unit is basically unsteady 
when there is ingestion of water along with the gas phase. The gas phase is a 
mixture of air and water vapor in the case of a bypass fan engine that provides 
thr ust power to an aircraft. The liquid water may be in the form of droplets and 
film at entry to the fan. The unsteadiness is then associated with the relative 
motion between the gas phase and water, at entry and within the machine, while 
the water undergoes impact on material surfaces, centriguging, heat and mass 
t rans fer processes, and reingestion in blade wakes, following peal off from blade 
surfaces. The unsteadiness may be caused by changes in atmospheric conditions, 
and at entry into and exit from rain storms while the aircraft is in flight. In a 
multi-stage machine, with an uneven distribution of blade tip clearance, the 
combined effect of various processes in die presence of steady or time-dependent 
ingestion is such as to make the performance of a fan and a compressor unit time- 
dependent from the start of ingestion up to a short time following termination of 
ingestion. 

The original WINCOF code was developed without accounting for the 
relative motion between gas and liquid phases in the ingested fluid. A 
modification of the WINCOF code has now been developed, named the 
WINCOF-I, which can provide the transient performance of a fan-compressor 
unit under a variety of input conditions. 

A complete description of the modifications introduced in the WINCOF-I 
code is provided. Along with the documentation on the WINCOF code, the 
current description provides the documentation on the WINCOF-I code. 

In order to illustrate the use of the WINCOF-I code for determining the 
performance of a multi-stage compressor, a test case with two stages has been 
calculated, and the details of the input and the output have been presented. 

Finally, a methodology for incorporating the output of the WINCOF-I 
code in a code for the determination of transient performance of a bypass fan 
engine is described. 
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CHAPTER 1 
INTRODUCTION 


Ingestion of water into jet engines during aircraft flight operations in rain 
storms (and also, during take-off over rough runways with puddles of water) has 
been known to cause changes in engine performance, and difficulties in flight 
operations. A knowledge of such performance changes is essential not only to 
improve flight quality and safety but also to enable advanced design, and engine 
testing and certification procedures. 

An approach to developing methods of predicting engine performance is to 
include in an engine transient performance code elements or codes for 
determining the performance of individual components of the engine with water 
ingestion. Over several years, attempts have been made to model axial-flow fan- 
compressor units of fan jet engines operating with air-water mixture flow, and to 
develop a performance code. The output of the code has then been included in an 
engine transient performance code. The fan-compressor performance code with 
water ingestion has been named WINCOF, and References 1-6 provide a 
bibliography of publication that describe the WINCOF code and it application. 

The WINCOF code was developed to provide a tool for determining the 
steady state performance of fan-compressor units with steady ingestion of an air- 
water mixture of fixed quality. Water in the mixture could be in the form of a 
film or in part in film and the rest in droplet form. When determining transient 
engine performance changes, the steady state output of the WINCOF code was 
utilized while the engine performance was time-dependent for other reasons. 

It has since been realized that the performance of an axial-flow fan- 
compressor unit is basically time-dependent whenever there is an ingestion of air- 
water mixture into it. The time-dependence is a consequence of the presence of a 
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relative velocity between air (and water vapor, or the gas phase, in general) and 
water at entry and in the flowpath of the mixture in the machine. It is also a fact 
that die blade geometry, the stage performance, and the blade clearance height do 
not follow any systematic variation in relation to the velocity of flow in the blade 
passages in a multi-stage machine. Therefore, the combination of various 
processes in the fan-compressor unit, such as impact of water on material 
surfaces, centrifugal action, heat and mass transfer, and reingestion of water 
following peal-off from blade surfaces into wakes, lead to a time-dependent 
performance of the unit. Such time-dependent performance extends from the 
instant of start of ingestion up to a short time following termination of ingestion. 

The performance of a fan-compressor unit can also be time-dependent for 
two additional reasons: (i) Due to changes in atmospheric conditions, causing 
changes in the characteristics of the ingested fluid mixture; and, (ii) due to such 
events as entry into and exit from rain storms during flight. 

The WINCOF code has been modified so as to permit taking into account 
the relative velocities between the two phases (water being in film and droplet 
form), and thus obtaining the transient performance of a fan-compressor unit. 
The modified code has been named the WINCOF-I code, and a brief description 
with an illustrative example has been provided for the modified code in Ref. 7. 

The current report provides a more complete documentation of the 
WINCOF-I code with necessary explanations. It also includes the methodology 
for utilizing the results obtained with the WINCOF-I code in determining the 
transient performance of a fan engine. A generic bypass fan engine was chosen 
for illustrative purposes in Ref. 8. 

1.1 Outline of Report 

The methodology for determining the time-dependent performance of a 
fan-compressor unit is described in Chapter 2. The method of incorporating the 
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results generated by utilizing the WINCOF-I code in an engine transient 
performance code is given in Chapter 3. A general discussion on the 
methodology is included in Chapter 4. Various details related to the WINCOF-I 
code are given in a series of five appendices, which are referred to appropriately 
in various chapters. 


CHAPTER 2 

FAN-COMPRESSOR UNIT PERFORMANCE WITH WATER INGESTION: 

METHODOLOGY 

The performance of an axial-flow fan-compressor unit under discussion in 
this chapter is that obtained for the unit when it is operated in isolation by means 
of independent drives as necessary. Most of the discussion pertains to a multi- 
stage compressor unit, while some of it is devoted to the combined fan- 
compressor unit driven by two independent drives, one for the supercharger with 
the fan and the other for the core compressor, at different speeds. 

The operational parameters for a fan or a compressor unit are (a) the 
operating speeds, (b) the total mass flow at entry, (c) the bypass ratio and (d) 
water ingestion parameters, in addition to the ambient conditions of pressure, 
temperature and degree of saturation with respect to humidity. The mass flow at 
a specified value of rotational speed at the radius (in blade span) under 
consideration may be specified in terms of a flow coefficient value. The flow 
coefficient is defined as the ratio of axial flow velocity of air to the rotational 
velocity of the rotor blade row (or stage), both at the same spanwise radius under 
consideration. The range of operating speed of the supercharger, as also that of 
the flow coefficient, differs from that of the core unit by design. The bypass 
ratio is a function of operating speeds of the two units and the power demand. 
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The water ingestion parameters are mass fraction of water taking into account the 
scoop factor, droplet size distribution, difference in velocity between air and 
large droplets and temperature. 

The performance parameters of interest in the case of a fan or a 
compressor unit at specified speed of operation, flow coefficient and water 
ingestion conditions are the following at the outlet of the compressor. 

(i) Pressure ratio, temperature ratio and efficiency (usually adiabatic 
efficiency based on output obtained for given input of work based on 
conservation of energy); 

(ii) mass fraction of water in the air-water mixture; 

(iii) temperature of water in span; and 

(iv) volumetric mean droplet size. 

(v) mass fraction of vapor in the air-water mixture; 

(vi) thickness and mass flux of water in film form in the compressor 
clearance. 

The foregoing may be of interest in any stage of a fan or a compressor or for a 
multistage unit as a whole. 

In the case of operation of a unit with air, it is common practice, in 
simplified approaches, to specify the performances of a compressor in terms of 
one or both of the following: (i) Performance of the unit with reference to a 

specified radius of the unit; and (ii) performance of the unit based on a 
weighted-average of performance data at several chosen radii to encompass the 
span of the unit. In the alternative, one can undertake fully three-dimensional 
predictions for the unit over its entire span. However, for the purposes of 
specifying representative performance parameters, especially for use in engine 
simulation codes, it is again necessary to determine, through some form of 
weighted averaging, the overall performance of the unit. 
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In the current investigation, the performance of an axial-flow fan or 
compressor unit is determined utilizing a modified version of the so-called 
WINCOF code (References 1), developed at Purdue University some years ago. 
A brief description of the WINCOF code follows in Section 2.1. That code has 
been modified in several respects (Reference 7) and the modified code is named 
WINCOF-I code, and is described in Section 2.2. 

In both of these codes, the performance of a fan or a compressor is 
determined with respect to a chosen streamtube. Figure 2.1 shows a set of six 
streamtubes that are illustrative of the type of streamtubes of interest in a fan- 
compressor unit. Among those streamtubes 2 and 5 are considered as mean 
streamtubes and the performance of various units calculated with respect to them 
are utilized as representative for the units. 

The WINCOF-I code has been utilized, for illustrative purposes, to 
determine the performance of various units in the fan-compressor unit of the 
generic engine referred to in Chapter I. 

2.1. Description of WINCOF Code 

The objective in the development of the original WINCOF code has been to 
obtain a means of predicting the steady state performance of an axial-flow 
compression system (single or multi-stage fan, compressor or combination) with 
water ingestion. The WINCOF code was developed to calculate, as stated earlier, 
the performance with respect to a designated streamtube of a single stage of a fan 
or compressor. Thus no account is taken of the radial component of velocity in 
the three-dimensional flowfield that is usually generated in an axial-flow 
turbomachine. The calculation procedure of the WINCOF code is, therefore, 
referred to as "one-dimensional’', implying that the flux of fluid is assumed to be 
in the axial direction along the designated streamtube. However, the input of 
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work absorbed by the fluid is obviously related to the change in moment of 
angular momentum of the fluid; the rotor tangential speed determines the angular 
momentum. The WINCOF code can be utilized for determining the performance 
of a multi-stage machine by the stage-stacking procedure with respect to a 
streamtube passing continuously through the machine. Typical streamtubes for 
such calculation are presented, as referred to earlier, in figure 2.1. A streamtube 
is specified by means of its (a) area of cross-section and (b) radius of location in 
each stage, with necessary attention to continuity and smooth transition from 
stage to stage. 

The WINCOF code can be utilized to determine the performance of a 
single- or multi-staged unit under (i) design conditions for operation with air and 
(ii) off-design conditions with (a) air or (b) air-water mixture. For design point 
calculation of performance during operation with air, the following geometrical 
and aerodynamic details of design are required for the unit. 

i) Hub and casing diameters; 

ii) spanwise distribution of chord and thickness-chord ratio; 

iii) spacing between blade rows; and 

iv) spanwise distribution of metal angles, incidence, deviation, inlet 
Mach number and maximum boundary layer momentum thickness along chord 
for inlet guide vane, rotor and stator, as appropriate, at the design point of the 
unit. 

For the calculation of performance at an off-design point, the following 
procedure is adopted. 

(a) Operation with air : Various rules are incorporated in die code for 
obtaining, based on design point information, the necessary aerodynamic 
performance parameters corresponding, for example, to the specified value of 
rotational speed, ambient conditions and flow coefficient at the designated 
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streamtube location. A typical set of such rules can be found Reference 1, and 
they have been incorporated into the code. Details are repeated for convenience 
in Appendix I. It may be pointed out that another set of rules can be easily 
substituted in place of the set of rules currently adopted. 

(b) Operation with air-water mixture : The air-water mixture modifies 
the aerodynamic parameters as well as introduces several new features that are a 
direct consequence of the two-phase nature of the mixture. They are described 
briefly in Section 2.1.1. along with the procedure for taking into account those 
two-phase fluid-related processes. 

2.1.1. Processes Included for Consideration 

The air-water mixture entering the fan or the compressor unit is 
characterized by (a) composition of air-water mixture and (b) velocities of air 
and water droplets. Regarding (a), further details are of interest concerning (i) 
mass fraction of vapor and water and (ii) droplet size distribution. The latter is 
of critical importance in determining drag, velocity of motion and interphase heat 
and mass transfer processes. However, in the current investigation, some 

simplification is introduced as follows. 

Droplets are considered to be in two categories based on size, namely large 
and small, the latter under 20 microns in mean-volumetric diameter. Small 
droplets are assumed to follow the gas path, while large droplets move 
independently. Also, small droplets absorb work input, while large droplets do 
not. Work is absorbed by water in the course of pealing off from a rotor by 
shearing action. Both sizes of droplets undergo size adjustment based on mass 
transfer and critical Weber number consideration, and both undergo heat and 
mass transfer processes. 
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Both the composition and the velocities may vary radially and 
circumferentially across the plane of ingestion. Radial variation can be accounted 
for by specifying the relevant values in various streamtubes chosen for 
calculation. However, circumferential variation requires consideration of a 
nonaxisymmetric flowfield, and this feature has not been included in die current 
investigation. 

The main two-phase flow processes that are considered to have significant 
effect on the performance of an axial-flow fan or compressor are the following. 

i) Ingestion of water at the machine face, 

ii) droplet impact and rebound from blade surfaces, 

iii) film formation and flow over blade surfaces, 

iv) flowfield and boundary layer modifications based on deviation, 
diffusion factor, and momentum thickness, 

v) centrifugal action on droplets and film over blade surface, and 
corresponding radial movement, 

vi) heat and mass transfer processes between liquid and gas phases, 

vii) reingestion of water into wakes of blades from film flow over 
surfaces, 

viii) droplet size and breakup based on the attainment of a critical value for 
the Weber number, 

ix) the total of work input and its division of work input between the two 
phases, and 

x) film formation and movement at the casing wall. 

In accounting for these processes, several assumptions are introduced as 
follows. 
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(a) Ingestion occurs at a specific plane immediately upstream of the unit 
or the stage under consideration. Initial conditions of air-water mixture, 
including composition and velocities, are specified at that plane. 

(b) Processes (ii), (iii) and (iv) are accounted for in establishing work 
done on the mixture and related losses. 

(c) Processes (v) through (viii) are taken into account at the exit plane of 
a blade row although they occur everywhere along the flow in the blade passage. 

(d) Regarding processes (v) and (vi), which are time rate-governed, it is 
necessary to determine a duration of time for the occurrence of the processes. 
This is done based on characteristic length and velocity scales, namely the chord 
of blades and axial velocity of mean air-water mixture flow through the relevant 
blade passage or stage width.. 

(e) The liquid phase may absorb work input only when the droplets are 
small and follow the gas path or when water is splashed off a blade surface. In 
general, the amount of work done on the liquid phase can be assumed to be small. 

(f) Finally, film formation and movement at the casing wall result from 
centrifugal action on droplets and film over blade surfaces, accumulation of 
water at the casing wall and shearing action due to adjoining air-water mixture 
flow. In the WINCOF code, it has been assumed that the film moves at the same 
speed as the adjoining air-water mixture. 

2.1.2 Single Stage Machine 

To apply the WINCOF code to a stage of a fan or compressor, a specific 
streamtube must be chosen, along which the performance calculations are 
performed. A streamtube is designated by its location along the span, and with a 
specified cross-sectional area. A small value of cross-section area is chosen 
based on local design point mass flow and density and axial velocity data 
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pertaining to the chosen radius of location. But the performance of a blade row 
at the design point is also a function of blade metal angle, incidence angle, and 
deviation angle. Also a stage consists of two or three blade rows. Therefore it 
becomes necessary to use some trial and error in choosing die location and cross- 
section area such that it is compatible with design point data available at discrete 
locations along the stage. 

2.1.3 Multi-Stage Machine 

As mentioned previously, the overall performance of a compressor is 
established by extending the single stage calculation through a "stage-stacking" 
procedure. In this connection, it is worthwhile discussing the manner in which 
centrifugal action on water is determined in various stages. In the WINCOF 
code, in order to calculate centrifugal action on water in a blade row or stage, the 
span of the blade row or stage has been divided into a finite number (ten, for 
example) of streamtubes of equal cross-sectional area. The centrifugal action and 
the resulting displacement of water towards the casing wall are calculated at the 
exit plane of the blade row or stage under consideration over the period of time 
equal to the residence time of the air-water mixture therein. All of the water 
displaced is assumed to accumulate in die casing and not in any of the streamtubes 
themselves. Heat and mass transfer processes as well as droplet size adjustment 
processes are also taken into account at the same exit plane. When the calculation 
proceeds from one row or stage to the next, the water and the water vapor 
contained in any one of these streamtubes is considered to reside in the 
streamtube of the same number in the next stage. The water film at the casing 
wall is also considered to move from one row or stage to the next at the same 
speed as the adjoining air-water mixture. However, in the WINCOF code no 
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account is taken of either (a) the casing clearance space or (b) relative motion 
between film and air-water mixture. 

2.1.4 Weaknesses in the W INCOF Code 

In recent years a detailed examination has been undertaken of the 
assumptions and methodology employed in developing the WINCOF code. A 
number of weaknesses have been recognized in the WINCOF code as follows. 

(i) The WINCOF code was developed and utilized as a code for the 
determination of the steady state performance of both single and multiple stage 
axial-flow machines. In accounting for time-rate-dependent processes such as 
centrifugal action on water droplets and film and inter-phase heat and mass 
transfer processes, a time scale was chosen equal to an estimated value of 
residence time of air-water mixture flow in the blade row under consideration. 
The residence time was defined as the ratio of mean width of a blade row or a 
stage to the mean axial velocity of flow in it. Since such processes were 
accounted for at the exit of a blade row or a stage, the performance so calculated 
was assumed to be the steady state performance of the stage with the given water 

ingestion and other inputs to the stage. 

(ii) Centrifugal action on water droplets and film (over blade surfaces 
or in blade passages) was considered only to establish the amount of ingested 
water that became removed in film form at the casing wall and hence to 
determine the balance of amount of water that remained within the span of the 
blade row. No account was taken of the motion of water film at the casing wall 
due to the shearing action exerted by the adjoining air-water mixture and the 
time-dependent nature of build-up and motion of the film at the casing wall. In 
other words, it was assumed in the WINCOF code, that the blade row or stage 
had attained a state of "equilibrium” or a steady state at the end of the (somewhat 
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arbitrarily selected) step time of calculation. In fact, taking into account the 
nature of centrifugal action over the blade span and the growth and motion of the 
clearance film, the distribution of water in the blade span would become tune- 
dependent. This was left out of account as well as the resulting time-dependent 
changes in the aerodynamic performance of die unit under consideration. 

(iii) Heat and mass transfer processes were accounted for in the WINCOF 
code at the exit of a blade row or a stage utilizing the residence time of air- water 
mixture as the time interval for the processes to occur. The interactive time- 
dependent changes in the composition and state of air-water mixture between 
centrifugal action, film formation and motion in casing wall on the one hand and 
heat and mass transfer processes on the other were again left out of account. 

(iv) No account was taken of the possible difference in velocity between 
water droplets and air at entry to a compressor. While small droplets were 
assumed (correctly) to move with air, large droplets were considered also as 
entering the machine with the same velocity as that of air. Thus the relative 
velocity between air and (all of) the droplets was assumed to be zero at entry to a 
fan or a compressor. 

(v) At entry to a fan or a compressor, two other features of air-water 
mixture that are of interest are (a) nonuniformities due to scoop factor and other 
causes and (b) presence of a film at the casing wall, for example, due to flow 
over the inlet wall. Neither of these was included in the utilization of the 
WINCOF code. 

(vi) Finally, interphase heat and mass transfer processes in a 
conglomeration of droplets are strong functions of the size distribution and 
number density of droplets. There are not enough reliable data for including the 
details of shielding of heat and mass transfer processes in a cloud of large 
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droplets; however, a parametric study could be conducted to determine the effect 
of shielding. This was not included in the WINCOF code. 

There are several implications of (i) to (v) in the use of the WINCOF code 
for determining the performance of a multistage machine. In actual calculation, a 
single sweep was made through a fan or a compressor, over all of its stages. In 
determining centrifugal action and heat and mass transfer, a calculation time was 
introduced for each stage equal, as stated earlier, to the residence time of air- 
water mixture in the blade row or stage under consideration. Then the 
performance of the machine obtained at the end of the single calculation sweep 
of time was assumed to be the steady state performance of the unit. The 
performance determined for any chosen streamtube pertained to the air-water 
mixture in it, taking into account centrifugal action and heat and mass transfer 
processes. Meanwhile, no further account was taken of the water at the casing 
w all or the consequences of its motion relative to the air-water mixture in the 
span. 

2.2 The WTNCOF-I Code 

The WINCOF-I code has been developed with the objective of removing 
some of the aforementioned weaknesses and limitations in the WINCOF code. 
The modifications introduced pertain to the following. 

i) Time-dependent nature of performance of axial-flow fans and 
compressors during water ingestion (Section 2.2.1.); 

ii) modeling centrifugal action and film formation and flow as a function 

of time (Section 2.2.2.); 

iii) modeling heat and mass transfer (Section 2.2.3); 

iv) accounting for droplet-relative velocity with respect to air velocity 

(Section 2.2.4); and 
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v) accounting for desired entry conditions to the machine with respect to 
air temperature and thickness of film at the casing wall. 

Each of these is discussed in some detail in the following. 

2.2.1. Time-Dependence in the Performance of a Fan or a Compressor 

A single or multiple stage fan or compressor cannot be expected to operate 
under steady state conditions with water ingestion may be shown as follows. 

Water entering a stage and moving over the surface of a blade or in the 
passage between blades becomes centrifuged to the casing wall on account of the 
rotational field. At the casing wall, a film is formed. The film moves along the 
casing wall due to the shearing force exerted on it by the adjoining air-water 
mixture in relatively high speed motion. The shearing force is obviously small 
and there arises a large difference between the velocities of the film and the air- 
water mixture. Since water is being centrifuged continuously, one can expect a 
net growth in the thickness of the film at any blade row location, while some of it 
moves along the casing wall to the exit or a following blade row location. This is 
equivalent to an accumulation of water in the clearance at any blade row. The 
accumulation ceases only when (a) the local, instantaneous clearance between the 
blade row and the casing wall becomes filled up with water or (b) the film 
thickness is such that the adjoining air-water mixture flow is no longer able to 
shear it to a recognizable velocity. In either case, any water centrifuged 
subsequently may only be splashed back into the span of the blade row. Such a 
state of operation of any blade row (or a stage, when a complete stage is under 
consideration) is named the ’’equilibrium” state of operation. Beginning with the 
instant of ingestion, the first or any other stage of a machine and also the entire 
machine may reach ’’equilibrium", if only instantaneously. It corresponds to a 
limiting condition of operation beyond which any water centrifuged may only be 


14 


splashed back into the span and cannot accumulate in the local clearance space. 
The performance of the stage or the machine is time-dependent during the 
attainment of equilibrium condition. It turns out that it will continue to be time- 
dependent, perhaps periodic, beyond that instant of time. 

In discussing redistribution of water due to centrifugal action, it is 
necessary to distinguish and to treat individually the following: (i) a stage, 

whether it is the first or a subsequent stage in a multistage machine, and (ii) the 
nature of ingestion, whether (a) steady, uniform, (b) steady, nonuniform, or (c) 
nonsteady uniform or nonuniform. The case of steady, uniform ingestion is 
considered first, and some remarks are added about the other cases at the end. 

The reasons for time -dependent performance during water ingestion may 
be summarized as follows. 

(i) Considering first a single stage machine or the first stage of a multi- 
stage machine and steady ingestion of water, one can then visualize a period of 
time-dependent performance during which the film at the casing wall is growing 
steadily and also moving downstream from the stage. Finally, at the end of a 
certain number of calculation time step intervals, an equilibrium state is reached. 
In the next calculation time step interval, any water centrifuged becomes splashed 
back into the span. If it is assumed that the water splashed back into the span 
becomes uniformly distributed over all of the calculation streamtubes in the span, 
the one or more streamtubes at the hub of the row that had originally been 
depleted of water due to centrifugal action would now again contain water due to 
the splash back just described. Meanwhile the steady, uniform ingestion at the 
front of the stage continues to bring in air-water mixture across the entire span. 
Hence the process in the blade row remains time-dependent so long as there is 
water ingestion. 


15 



(ii) During such a period of time, the effects of heat and mass transfer 
processes - in fact, the processes themselves - become time-dependent. Thus 
there arises a further time-dependent effect. 

(iii) Considering next the second or any subsequent stage of a multi-stage 
fan or compressor, there are two inputs of water at the casing wall, one from the 
previous stage and the other due to centrifugal action in its own span. Both of 
those are time-dependent. Along with the influence of heat and mass transfer, the 
processes in any but the first stage of a multi-stage machine remain time- 
dependent. 

One can next ask, if steady state conditions are ever likely to be reached in 
a multi-stage machine. A multi-stage machine has a reducing area of cross- 
section from inlet to outlet in view of pressure and density increases along the 
compressor. Individual stages, in general, have different values of casing 
clearances, blade-heights, geometrical features and aerodynamic characteristics. 
In any practical machine, in general, those quantities do not vary according to 
any organized, simple, functional relationship from stage to stage. On the other 
hand, it is clear that, in order that a steady state may be attained, at least the 
following quantities must vary linearly with respect to each other along the 
compressor: (a) residence time of air-water mixture, equal to stage width 

divided by mean axial velocity over the stage, (b) span height, hub radius and 
blade clearance, (c) aerodynamic performance, (d) heat and mass transfer and (e) 
droplet break-up and reforming characteristics. Such an occurrence is extremely 
unlikely, if not impossible, in any practical machine. Thus, for a steady, uniform 
ingestion of a given type of water ingestion, one invariably obtains a time- 
dependent performance for a multistage fan or compressor. 

There is one other feature of a multistage fan or compressor that must also 
be recognized in regard to its time-dependent performance with water ingestion. 
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This relates to the fact that in a given multistage machine, succeeding stages may 
not necessarily reach their "limiting" conditions one after the other, that is in the 
same order as they appear, the third after the second, the (n-l)th after the (n-2)th 
stage. By "limiting" conditions, it is implied that a steady set of conditions has 
come into existence with respect to performance of the stage including the growth 
and motion of water film in the local casing clearance. If such a state is not 
reached in succeeding stages in succession, then the conditions at the exit plane of 
the machine will definitely have a time-dependent character. However, at the end 
of a sufficiently long period of time, all of die stages will have attained limiting 
conditions for an instant of time. But at the very next instant of time, further 
ingestion (of the same steady type and rate) will cause a change in some of the 
stages, while not in others, and therefore the time dependence will persist 
although with a long-time, periodic behavior or pattern. 

The foregoing applies to steady, uniform ingestion of water. If the 
ingestion is steady but nonuniform, the extent to which different calculation 
streamtubes of the stage receive water becomes different. This leads to changes, 
compared to the case of uniform ingestion, in (a) the distribution of water across 
the span at the end of each calculation step, (b) the interval of time required for 
the attainment of equilibrium conditions for the first time and (c) the periodic 
nature of nonsteady performance. 

In the case of nonsteady ingestion of water, die processes become further 
complicated depending upon whether ingestion is (a) continuously or discretely 
nonsteady and (b) simultaneously, uniform or nonuniform. No general 
statements can be made in regard to such situations, especially because of the 
complicated influence of stage aerodynamic characteristics and because of the 
presence of various values of clearance in different stages. 
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In summary, the time-dependent nature of performance of the first or any 
other stage of an axial-flow compression machine or of the machine as a whole 
depends upon the following. 

(a) Type of water ingestion; 

(b) geometry of various stages including the casing clearance; 

(c) aerodynamic characteristics; and 

(d) heat and mass transfer processes. 

It can be stated that it is extremely unlikely that steady state performance becomes 
possible in any machine with water ingestion. 

2.2.2 Model for Centrifugal Action and Film Formation and Flow 

The WINCOF code does not take into account the motion of the film due to 
the shearing action of the adjoining air-water mixture, the relative motion 
between the film and the air-water mixture, and the resulting time-dependent 
processes until equilibrium conditions become set up in the stage. In order to 
account for these, the following procedure is adopted in die WINCOF-I code. 

(a) Single stage machine or the first stage of a multi-stage machine. 

(i) Calculation of centrifugal action: The blade span is divided into 
a finite number of parts, for example 10, based on (a) size of droplets present, 
(b) increment in thickness of film at casing wall resulting from radially displaced 
water during a calculation step in time and (c) depletion in film thickness on 
account of motion of water in the film caused by the shearing action of the 
adjoining air-water mixture, again in the same calculation step in time. The 
number of parts into which the span is divided has to be chosen by a trial-and- 
error procedure. Taking the local rotational component of motion into account, 
the amount of water displaced from each calculation streamtube to the next is 
determined as well as the net water content in itself. 
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(ii) Film formation and motion: Water displaced from the 

calculation streamtube in the tip region of the blade is allowed to accumulate in 
the casing clearance. The film so formed becomes subject to the shearing action 
of the adjoining air-water mixture. Based on the momentum transferred to the 
film, its velocity is determined. 

(iii) Net state of film: By knowing the velocity of film motion and 
the resulting depletion of mass of water in the film, the final thickness of the film 
is determined. 

(b) The second or other stage of a multi-stage machine. 

In carrying out a sweep of calculation across the machine, each stage of the 
machine is dealt with assuming the following inputs. 

(i) The state of air-water mixture at the exit of the previous stage 
after all of the processes in that stage have been taken into account; 

(ii) the state of air-water mixture across the span in the various 
calculation streamtubes at the exit of the previous stage taking into account any 
depletion that has occurred in water in any of the streamtubes, including the one 
near the hub; and 

(iii) the state of the film in the clearance that is entering the current 
stage from the exit of the previous stage. 

The rest of the procedure for the calculation of centrifugal action and film 
formation and motion is the same as for the first stage. 

Details of the analysis and procedure are given in Appendix n. There are 
also some implications of the foregoing discussion that are of interest in the 
numerical procedure adopted. These are discussed in Sections 2.2.6 and 2.3. 
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2.2.3 Model for Heat and Mass Transfer 

In the WINCOF code, interphase heat and mass transfer was calculated 
based on (a) fixed coefficients of heat and mass transfer and (b) a time interval 
equal to the width of a blade row or a stage divided by the mean axial velocity 
therein. These calculations were preformed at the exit of a stage as a correction 
to the aerodynamic performance obtained for the stage. 

In the WINCOF-I code two modifications have been introduced as follows. 

(i) Shielding of heat and mass transfer processes due to multiplicity of 
droplets: In a conglomeration of droplets, heat and mass transfer are affected, 
along with other parameters, by (a) the number density, (b) the size distribution, 
and (c) the vapor clouds around the droplets, which grow in size with time. In 
order to take these into account in a parametric form, it is assumed that (a) there 
is a maximum separation of droplet clouds beyond which droplets do not interact 
with each other; for example, 10 volumetric mean diameters, and (b) there is a 
minimum separation of droplets less than which no heat and mass transfer is 
possible; for example, 2 volumetric mean diameters. Figure 2.2 illustrates 
schematically the nature of the assumptions. Based on such conditions, a linear 
curve has been constructed to relate a so-called heat and mass transfer weighting 
parameter, z, to the mean separation between droplets utilizing mean volumetric 
droplet diameter as an additional parameter. The curve is presented in figure 2.3 
It is then possible to utilize this parameter in the calculation of heat and mass 
transfer at any location when droplets of a given mean volumetric diameter are 
present therein at a specified number density. 

(ii) Characteristic length of time for heat and mass transfer processes: 
The residence time of droplets in a blade passage depends upon the velocity of 
droplets therein. In general, the residence time of droplets should be several 
times larger than the residence time of the gas phase. In order to take this into 
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account parametrically, a time factor has been introduced in WINCOF-I, whose 
value may vary over the range 1 - 100. The length of residence time of air- 
water mixture in a blade row or a stage is then multiplied by the time factor and 
the resulting time utilized in the calculation of interphase heat and mass transfer. 

Details concerning the calculation of heat and mass transfer in the 
WINCOF-I code are given in Appendix HI. 

It may be pointed out here that heat and mass transfer processes are 
confined entirely to the span region of a compressor unit. Thus no heat and mass 
transport are included for the film in the casing clearance space; this is based on 
the reasoning that the film, growing continuously in thickness with centrifuged 
water, does not undergo transport processes. 

2.2.4 Relative Velocity of Droplet with Respect to Air 

There are two situations where the relative velocity between the droplets 
and the air become of interest: (i) at entry to a turbomachine and (ii) between 
two blade rows. In case (i), as stated earlier, the velocities of large droplets may 
be quite independent of the velocity of air. Regarding (ii), water droplets are 
present in the blade passage, and are also generated from water pealed off from 
blade surfaces into blade wakes. In either case, the attainment of mechanical 
equilibrium and, thereby, a particular size distribution requires a length of time 
and distance. In other words, there is a lapse time of residence in which water is 
essentially in an indefinite state and, therefore, during this time water can be 
assumed to be moving at a substantially lower velocity compared to the velocity 
anywhere else. 

The difference in velocity between droplets and air was entirely neglected 
in the WINCOF code. In the WINCOF-I code it has been taken into account 
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parametrically by introducing a so-called relative velocity factor, The relative 
velocity factor is defined by writing 

$ - (V a - V D )/V a 

where V a and Vd denote velocities of air and droplet, respectively. Thus £=1 
represents the condition Vd = 0; and ^=0 corresponds to Vd = V a . 

2.2.5 Other Parameters 

Other quantities that have been taken into account as parameters in the 
WINCOF-I code are as follows. 

(i) Ambient water temperature: The effects of large differences in 

temperature between air and water have been explored. 

(ii) Film thickness at entry to the machine: The effects of the presence of 
a film of finite thickness at entry to the machine have been explored in addition to 
nonuniformities in water ingestion parameters. 

(iii) Ingestion as a function of time: The mass fraction of water ingested 
may vary as a function of time in the form of a series of telegraphic signals of 
varying magnitude. It is then of interest to investigate the nature of changes in 
performance as a function of the manner in which the mass fraction of ingested 
water changes with time, both on short time and long time basis. 

2.2.6 Numerical Procedure and Modifications Introduced in WINCOF-I Code 
No modifications of the WINCOF code have been required in regard to the 

parameters mentioned in Section 2.2.5. The initial conditions are specified as 
necessary at entrance to the first stage of the machine under consideration. 
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Modifications and additions made in various subroutines, in the WINCOF-I 
code, relative to those in the WINCOF code, are discussed in Appendix IV . 

Concerning the numerical procedure, the main parameter of interest is the 
set of calculation times in determining the transient performance of the machine. 
The problem is formulated as follows: air-water mixture at given state conditions 
and composition enters the machine at the instant of time t = to • The residence 
time of the air-water mixture in a blade row or a stage is assumed, as stated 
earlier, to be equal to the width of the element of machine divided by the mean 
axial velocity therein, AtRS secs. Centrifugal action and film build-up and 
motion are calculated over the interval of time AtRS secs. At the exit of the 
element of machine under consideration, one obtains (a) the state conditions of 
air-water-vapor mixture in the calculation streamtube under consideration, (b) 
the extent of span from the hub over which water has become depleted due to 
centrifugal action, and, therefore, the extent of span over which water remains up 
to the tip of the blades, (c) the extent of clearance that has become filled up with 
water film in the element of machine under consideration and (d) the film flow 
conditions. Utilizing those as entry conditions to the next element of machine 
(blade row or stage) the calculations are repeated. In this way, the entire 
machine is swept once over an internal of time that is equal to the sum of the 
residence times of air-water mixture in all of the elements of machine, the so- 
called sweep time, AtRM* The calculation must now be repeated utilizing given 
ingestion conditions and starting with the entry section of the machine in order to 
obtain "equilibrium" conditions of film formation and flow in all of the stages. 
This requires, in general, a large number of sweeps across the machine until, at 
the end of an interval of time, AtEM > tdl of the stages indicate operation under 

equilibrium conditions. 
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It may be useful to recall here that, depending upon the geometrical and 
aerodynamic characteristics of the machine under design and off-design 
conditions, various successive stages may not attain equilibrium in succeeding 
order. Then at the end of time t = to + AtEM . there arises equilibrium in all of 
the stages at the end of a calculation sweep time equal to AtsM; however, when 
the calculation is extended further in time, starting with entry to the machine, one 
again observes time-dependence in the exit conditions from the machine. The 
exit conditions of interest are the same for any stage, including the last stage, as 
those given for the first stage. When the calculation is carried out over an 
interval of time equal to several times the AtEM. one obtains a periodic pattern of 
performance with equilibrium conditions arising instantly at the end of various 
successive values of AtEM • It must be pointed out here that the foregoing 
calculation intervals of time and the resulting performance estimates may vary 
substantially in different machines of differing aerodynamic and geometrical 
design, rotational speed, and physical size and also, with different air-water- 
vapor mixture conditions and their distribution with respect to span at entry to 
the machine. 

In summary, the calculation of transient performance of a multistage 
machine involves several characteristic calculation times, namely AtRS , AtRM > 
AtEM • In any machine three types of performance of a stage may be found even 
under equilibrium conditions: (i) clearance partially filled but film unable to be 
sheared by the adjoining air-water mixture flow, (ii) clearance completely filled 
and film itself able to be sheared, and (iii) clearance completely filled and film 
also unable to be sheared. In each case, any water centrifuged beyond setting in 
of that state is splashed back and redistributed uniformly across the span of 
blading, by assumption. The manner in which various stages with different kinds 
of performance as described happen to become stacked in a given machine under 
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a particular set of operating and ingestion conditions determines how AtEM varies 
in the long time as calculations are continued in time. Such variation is affected 
further when there is radial nonuniformity in ingested water and also when 
ingestion itself involves an unsteadiness. 

An important observation here pertains to the calculation procedure based 
on using a finite difference approach in the calculation of transient performance. 
The discrete intervals involved in introducing AtRS and AtSM l ea d to changes in 
performance of various stages and the machine with a periodicity corresponding 
to those intervals. Such periodicity or unsteadiness is entirely due to the finite 
difference-based calculation procedure and has no physical significance. The 
periodicity of interest is only that associated with the instantaneous attainment of 
equilibrium, at the end of AtEM , in all of the stages during a long time operation 
of the machine with water ingestion. The machine operates, (even) with a steady 
set of conditions at entry to it, in a time-dependent maimer with a periodic 
behavior in the long time, the periods being of the order of AtEM • 

Finally, the procedure and the nature of results described in the preceding 
pertain to the case of steady conditions of air-water ingestion into the machine. 
There are various situations in practice in which the conditions may vary as a 
function of time, either continuously or, more often, over discrete intervals of 
time, in telegraphic signal form. It is then necessary to introduce changes in the 
entry conditions at the appropriate time and proceed as before. Such changes will 
invariably give rise to changes in all of the three characteristic calculation time 
intervals and also, in the long time performance of the machine. For example, 
AtRS Is a function of mean axial velocity in the element of machine under 
consideration and therefore, is a function of local aerodynamic performance 
which, in turn, is influenced by the local state of air-water vapor mixture 
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composition. Once there is a change in AtRS in any part of the machine, there 
also arise changes in AtRM and AtEM • 

2.3 Prediction Methodology 

The methodology employed in determining the transient performance of an 
axial-flow compression machine may be summarized as follows for a given set of 
input and operational data. It may be recalled that (a) the WINCOF-I code is 
operated in the transient mode considering a series of successive, discrete 
intervals of time; (b) centrifugal action, heat and mass transfer processes and 
droplet size adjustment are considered at the exit plane of a blade row or a stage 
following the determination of aerodynamic performance through the element of 
machine; and (c) performance of a multi-stage unit is obtained by adopting a 
stage-stacking procedure. 

(i) The streamtube for which the calculation has to be carried out is 
identified and its area of cross-section and radius of location are selected through 
the element of machine under consideration. 

(ii) The WINCOF-I code is utilized to determine the aerodynamic 
performance of the element of machine with respect to the chosen streamtube. 

(iii) In order to determine the effects of centrifugal action, the span of the 
blading is divided into a finite number of parts, each part being sufficiently small 
in order to obtain the desired levels of accuracy and, at the same time, adequately 
large with respect to the largest droplet size and radial displacement of water 
expected. The selection of the number of parts into which the blading is divided 
involves a certain amount of trial-and-error. However, in most of the problems 
dealt with to date in relation to the generic engine unit, division into 10 parts 
have proved adequate. 
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(iv) A time interval is chosen for the calculation of the effect of 
centrifugal action. As stated earlier, this interval of time has been chosen as 
equal to AtRs , the residence time of air- water mixture in the element of 
machine. Centrifugal action is calculated over this interval of time in all of the 
ten parts of the span. As a result, a part of the water in the element of span 
closest to the casing moves into the casing clearance space. At the same tune, in 
each of the elements of span starting from the hub, one obtains values for the next 
amount of water taking into account what is removed from it into the adjoining 
element towards the tip and what is added on to it from the adjoining element on 
the other side towards the hub. Insofar as the streamtube under consideration is 
concerned, one obtains the net amount of water present in it following centrifugal 
action. 

(v) Utilizing the same interval of time and taking into account the net 
amount of water present following centrifugal action, the heat and mass transfer 
processes are calculated. 

(vi) Finally, the size of droplets that can be present in equilibrium with the 
local state of gas phase is determined for the streamtube under consideration at 
the exit plane of the element of machine. At this stage, one has all of the data that 
are needed as input for the determination of the performance of the succeeding 
element of the machine. 

(vii) The one remaining parameter of interest is the net amount of water 
that is present in film form in the casing clearance. One can establish this taking 
into account the amount of water centrifuged into that space from the nearest 
element of span and the amount of that water which is sheared by the adjoining 
air-water mixture over the same period of time, namely AtRs- When the 
clearance space at the start of calculation is already filled up to the extent that (a) 
shearing of the film is no longer possible or (b) the clearance space itself is filled 
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up by the film, then any centrifuged water must become splashed back into the 
span and the thickness of the film remains the same as the input value. At the 
same time the excess amount of water that is splashed back is assumed to be 
redistributed uniformly across the span. This affects the value of (a) the mass of 
water that is calculated as present in the streamtube (b) the interphase heat and 
mass transfer and (c) the droplet size adjustment, calculated according to (iv), (v) 
and (vi) above, respectively. Thus, some iteration is indicated in establishing the 
final state of air-water-vapor mixture in the streamtube at the exit of the element 
of the machine. 

(viii) One can then proceed to calculations for other elements of the 
machine through a simple stage-stacking procedure for the streamtube under 
consideration. The streamtube must be located in each element of the machine in 
terms of its area of cross-section and its radius of location. When, according to 
this procedure, one has reached the exit plane, one has, by definition, completed a 
sweep over a period of time equal to AtRM- 

(ix) The next step is to repeat the sweep starting from entry into the 
machine utilizing the input data concerning ingestion and operating conditions, 
but taking into account the current state of the air-water mixture in the machine 
in various elements of it, as given by the results from the previous calculation 
sweep. The sweeps are continued until all elements of the machine reach 
equilibrium condition. 

(x) A long time calculation is now initiated over a period of time equal 
to several times the period AtRM » in order to determine the periodic nature of 
the time-dependent performance of the machine. It may be pointed out that this 
calculation may be of interest for certain selected stages of the machine or for the 
machine as a whole. In practice, mainly in order to reduce computational costs, 
calculations have been performed in the current calculation ranging over 100 - 
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1,000 AtRM* Generally calculations over 500 - 3,000 sweeps have proved 
adequate to uncover the periodic nature of performance in most of the cases 
investigated. It may be pointed out that in the multistage compressor chosen for 
illustration the period of time for 500 sweeps corresponds to about 3 seconds of 
actual operation time at rotational speed of 100 per cent of design speed, flow 
coefficient of 0.450 and ingestion of 4.0 per cent mass fraction of water in 
saturated air. 

(xi) When it is desired to take into account a non-zero value of velocity 
between water and the gas phase (air-vapor mixture) at entry to the first or any 
other element of the machine, the input as utilized in (i) above needs to be 
modified. 

(xii) When it is desired to take into account modifications in heat and mass 
transfer coefficients based on number density of droplets, the calculation scheme 
as given in (v) above needs to be modified taking into account the heat and mass 
transfer weighting factor. 

(xiii) If the entry input conditions to the machine vary as a function of 
time, such changes are taken into account in the input as utilized in (i) above at 
the instant of time of appropriate step advance in time. 

2.4. Application of WINCOF-I Code 

The application of the WINCOF-I code for the determination of a 
compressor performance with water ingestion is illustrated in Appendix V of the 
current report. 
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Figure 2.2. Schematic of a conglomeration of droplets. 
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Figure 2.3. Relation between mass transfer parameter (Q and mean separation 
between droplets. 



CHAPTER 3 

BYPASS ENGINE PERFORMANCE WITH WATER INGESTION: 

METHODOLOGY 

The transient performance of engines is of great interest in all cases when 
operational conditions are a function of time. It is also of interest when the input 
conditions to the engine are substantially different from those assumed in the 
design of the engine. During water ingestion the input conditions become 
different from those for which most engines are designed, namely operation with 
air. Water ingestion may occur under various steady and changing operating 
conditions, the latter in ambient conditions, and power demand as a function of 
time. Water ingestion parameters themselves may become affected by such 
changes through modification of scoop factor, for example. Finally water 
ingestion may be a function of time in certain cases due to changes in atmospheric 
conditions. In all such cases, it is of importance to determine the extent to which 
the performance, operation and handling characteristics become affected as a 
function of time. The setting-in of critical and unstable conditions needs 
attention. The surge margin and the margin in combustor performance with 
respect to the occurrence of flame-out are examples of critical conditions. 
Instability may lead to such and other undesirable situations. 

Early work on simulation of bypass engine performance under a variety of 
water ingestion conditions has been discussed in References 4 and 5. It was 
undertaken to determine the influence of several parameters related to water 
ingestion on the transient performance of a typical bypass engine, illustrated in 
figure 3.1. The parameters pertained to (a) compressor performance with water 
ingestion, determined utilizing the WINCOF code, (b) combustor performance 
and (c) instrument response during flooding with ingested water. An engine 
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performance simulation code was utilized to determine the transient performance 
of a generic bypass fan engine. 

In the following, the manner in which the output from the WINCOF-I code 
can be incorporated into the engine performance code is discussed briefly. 

3.1. Generation of Parameterized Compressor Maps for Use in an Engine 

Simulation Code 

The engine simulation code has provision to incorporate the performance 
of the fan-compressor unit (fan, lower pressure (supercharger) compressor and 
high pressure (core) compressor) in the form of a set of parameterized maps. 
The maps consist of the following for each component unit. 

(i) Loss as a function of work coefficient; 

(ii) minimum loss and minimum loss work coefficient as functions of 
rotational speed of machine; 

(iii) (loss - minimum loss) as a function of (work coefficient - minimum 
loss work coefficient) squared; 

(iv) minimum loss coefficient as a function of rotational speed of machine; 
and 

(v) pseudo-Mach number of flow as a function of (work coefficient - 
minimum loss work coefficient). 

The first three maps provide a representation of efficiency and the latter 
two, a representative of mass flow. Thus mass flow, work done and efficiency 
are represented as functions of speed in a combined fashion such that (iii) and (v) 
become available in a multi-linear form. This makes it relatively simple to 
introduce and to use the maps in the engine simulation code. Those maps and the 
methodology for determining them have been described in detail in References 4 
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and 5 and in summary form in Appendix VI. It may be pointed out that the 
engine performance characteristics presented in References 6 and 8 are based on 
the use of performance results for the fan-compressor unit that were obtained 
utilizing the WINCOF code. 

3.1.1 Use of Fan-Compressor Unit Performance Maps in the Engine Simulation. 

Code 

In generating performance of a fan-compressor unit for use in engine 
simulation with water ingestion, it is necessary to cover a sufficiently wide range 
of operating speeds and mass flows. In general, it is not feasible to determine a 
priori the ranges required since ingestion may cause the engine to operate in 
unusual ranges of operation parameters. Some guidance may be obtained from 
calculations performed for operation with air, and this is the basis employed in 
generating maps in the current investigation. 

Typical parameterized performance maps for a multi-stage compressor 
unit with water ingestion are illustrated in figures 3.2. It may be observed 
therein that the performance is given for a series of discrete values of mass 
fraction of water entering the unit. It is clear that engine simulation is feasible 
only for steady or discretely changing (from one value to another) amount of 
water ingestion at the front of the fan-compressor unit. In view of the nonlinear 
changes in performance with water ingestion, no simple interpolation is possible 
between performance values calculated for even closely separate values of 
ingested parameters. 

There arises yet another important consideration in the use of performance 
maps for the fan-compressor unit, namely that they refer to quasi-steady state 
operation during water ingestion. As stated earlier, at the end of a certain interval 
of time beginning with the inception of ingestion, a quasi-steady is reached in 
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which the aerodynamic performance of the unit (as a whole) remains nearly 
steady. However, it is necessary to recognize two features in regard to 
performance: (i) the distribution of water across the span and in the casing 
clearance continues to vary with a periodic feature and (ii) during the time 
required to attain quasi-steady state the performance of the fan-compressor unit 
can be appreciably worse than on reaching that state. 

On one hand, the periodic variation of water distribution has to be taken 
into account in the prediffuser, and combustor and also, if necessary, in the 
turbine and the thrust nozzle. 

On the other hand, the variation of performance during the attainment of 
quasi-steady state is not available in the performance maps corresponding to the 
quasi-steady state. In determining the transient performance of an engine, the use 
of "constant" (or non-time variant) quasi-steady maps, therefore, does not 
account for the variation of performance during the setting-in of the quasi-steady 
state. This method of determining performance becomes questionable in, at least, 
two cases considered in the current investigation as follows: (i) prediction of 
engine performance desired from the instant of commencement of ingestion; and 
(ii) prediction of engine performance from the instant of a change in the amount 
of water being ingested. In both of those cases, it is clear that the fan-compressor 
unit is in a non-steady mode of operation over an initial period of time during 
which conditions are tending towards the setting-in of a quasi-steady state. 

The performance of a fan-compressor unit with water ingestion is, 
unfortunately, highly nonlinear. It is, therefore, not possible to rationalize 
changes in performance during the period required for attaining the quasi-steady 
state in the form of functional relationship governing performance with respect 
to time even for a single unit. Even for a single unit, it is necessary to take into 
account rotational speed of operation, mass flow of air and mass fraction of water 
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being ingested as three (primary) variables governing performance; each of them 
has effects intimately intertwined with the other two. 

Two ways of overcoming this difficulty may be considered; (i) to 
introduce transient performance of the fan-compressor unit directly into the 
calculation of transient performance of the engine; and (ii) introducing some 
form of approximation for the performance over the initial period of time 
beginning with ingestion or change in amount of ingestion. Concerning (i), it has 
already been stated that current computer capability does not make it practicable. 
In fact, the setting up of quasi-steady performance maps is in itself a consequence 
of such limitations. Regarding (ii), the following method has been adopted in the 
current investigation. 

The change in performance with water ingestion, relative to the base 
performance for operation with air, is expressed in terms of adders in the engine 
simulation program. They are called in as necessary with reference to the mass 
fraction of water entering the compressor. The adders are obtained from 
performance maps of the type illustrated in figures 3.2 - 3.4. In using the adders 
for calculations immediately following the start of ingestion or a change in it, the 
following assumption is introduced: for any value of water ingestion, the adder 
becomes double the value varying linearly over an interval of time of about one 
minute and then, the original value applies. This is based entirely on trends in the 
results of various performance calculations conducted on a multi-stage 
compressor unit operating with several types of air-water mixtures. It is possible 
that other assumptions, not proposed here, may be preferred in different engines. 
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3.2 Predicted Results 

The methodology for predictions of performance for a bypass engine 
under various conditions of ingestion and operation, and specific results 
pertaining to an example engine are presented in Reference 8. 
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Figure 3.1. Schematic representation of a generic bypass fan engine. 
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Figure 3.2. (continued) 
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CHAPTER 4 
DISCUSSION 


The methodology developed in the WINCOF-I code can be utilized for the 
prediction of performance of an axial-flow fan or axial-flow compressor in a 
variety of cases of steady and varying ingestion of water. The type of cases that 
are significant in practice are discussed in Section 4.1. The WINCOF-I code can 
be incorporated into an engine performance code for the determination of the 
transient performance of an engine. Various cases of practical interest are 
discussed in Section 4.2. Finally, several recommendations for the manner in 
which the WINCOF-I code can be adapted for the improvement of design and 
operation of turbomachinery as well as a complete engine are presented in 
Section 4.3. 

4.1 Application of WINCOF-I Code 

The WINCOF-I code provides a means of determining the performance of 
an axial-flow fan or axial-flow compressor given (a) the geometry of the 
turbomachinery, (b) the aerodynamic rules governing the performance, and (c) 
the initial conditions of flow immediately upstream of the machine. The latter, 
namely the entry conditions, may vary in a variety of ways, and also differ 
substantially from one case to another of water ingestion into engines. 

The entry conditions may be characterized by the following parameters: 

(i) the variation of temperature and velocity of air as a function of blade 
radius (or height) and time; 

(ii) the mass fraction of water in air as a function of blade radius (or 
height) and time; 
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(iii) the ratios of water in film, droplet, and vapor form in different parts 
of the annulus formed by the blade height, including the casing and the hub walls, 
and their variation with respect to time; 

(iv) the droplet size distribution in the annulus as a function of time; and 

(v) the velocity of film and droplets across the annulus height as a function 
of time. 

In practice, the foregoing parameters are likely to vary on account of the 
following: 

(i) distortion introduced due to atmospheric conditions; 

(ii) distortion introduced due to inlet, spinner, and splitter plate geometry, 
or, equivalently, the scoop factor, as a function of time; 

(iii) film formation over inlet, spinner, and splitter plate; 

(iv) engine demand for air flow; 

(v) entry into and exit from environmental conditions giving rise to water 
ingestion; and 

(vi) changes in rainfall, atmospheric conditions, and engine demand in 
various combinations giving rise to modifications in ingestion as a function of 
time. 

In each of the foregoing cases, the WINCOF-I code can be utilized for 
determining turbomachinery performance. 

4.2 Applications in Engine Performance Estimation 

The WINCOF-I code can be readily incorporated into a code for the 
determination of the transient performance of an engine, as described in Chapter 
3. The engine performance is a function of the following parameters: 

(i) the basic design and matching of various components of the engine; 

(ii) the atmospheric conditions in the environment of operation; and 
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(iii) the thrust power demand from the engine with a given fuel flow 
control system. 

In practice, the latter two sets of parameters, (ii) and (iii), may vary for 
the same (practical) reasons as stated in section 4.1. These include time- 
dependent power demand changes during acceleration and deceleration of a flight 
vehicle with possible concurrent changes in flight altitude. In all such cases, the 
engine performance changes, relative to design point performance, can be 
established taking into account changes in turbomachinery performance utilizing 
the results of performance calculations performed with the WINCOF-I code. In 
addition, combustion performance changes also may be taken into account on a 
parametric basis, as stated in Chapter 3. 
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APPENDIX I 


RULES FOR ESTIMATION OF 
OFF-DESIGN PERFORMANCE OF COMPRESSORS 

The problem of determining the off-design performance of a multistage 
compressor may be stated as follows: Given the geometric details of the 

compressor, including the metal angles at inlet and outlet, and the aerodynamic 
design of the stages at the design point, to determine the performance of the 
compressor under conditions of operation away from the design point. The 
design point and, therefore, also off-design points for a stage may be designated 
by specifying the rotational angular speed of the rotor and the axial velocity or 
mass flux, as as function of radius across the span of the rotor in the stage, and 
the ambient conditions. For a multistage compressor, the axial velocity in the 
first stage may be specified as the reference value for the machine as a whole so 
long as there is no change in mass flux across the compressor. 

The main performance parameters of a compressor stage are pressure ratio 
and efficiency as functions of the operating rotational speed, and the axial 
velocity at entry to the stage. In the case of a multistage unit the pressure ratio 
and efficiency for the unit as a whole are functions of the operating rotational 
speed and the axial velocity at entry to the first stage, so long as the mass flux 
remains constant everywhere along the compressor. The axial velocity at entry 
to each stage succeeding the first is determined by the performance of the 
preceding stage. Thus, one can proceed to determine the performance of each 
stage. Then the performance of the multistage unit can be obtained by a stage 
stacking procedure, as done in the current case, or some other means of matching 
successive stages. 
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It may be noted here that, as far as ambient conditions are concerned, the 
operational parameters of a compressor stage may be specified in terms of 
angular speed and mass flux that are normalized with respect to reference values 
(for example, standard temperature and pressure). In the alternative, the 
performance of a stage can be given in a form that does not involve ambient 
conditions; this alternative procedure is discussed further in Appendix IV. 

Figure 1.1 provides a schematic representation of a part of a blade row, 
wherein various relevant parameters are shown. 

The geometry of a compressor stage is usually fixed. However, there are 
cases in which the stage setting (or its stagger angle) and the casing clearance are 
adjustable for various reasons. In such cases, one knows a priori, as part of 
design information, the manner in which blade geometry varies with off-design 
conditions of operation. 

The problem of predicting the performance of a compressor stage under 
off-design conditions thus becomes one of establishing the aerodynamic 
parameters that, in turn, can be utilized to calculate the performance 
corresponding to the given off-design operating conditions. 

The aerodynamic parameters may be chosen, for example, as the incidence 
angle, the deviation angle and a parameter by means of which the losses in the 
flow over the blades in the blade rows of the stage under consideration can be 
estimated (References 1 and 9). It is then necessary to evolve a set of functional 
relations, rules, as called here, between the changes in the aerodynamic 
parameters and the changes in the operating conditions, both with respect to 
design point conditions. 

The rules currently utilized are a modification of the rules given in 
Reference. They are provided in the following. 
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(1) Diffusion factor 

A key parameter in the estimation of aerodynamic performance of a 
compressor stage is the so called diffusion factor (Reference 9). Briefly, it 
denotes the diffusion of air as a function of momentum diffusion and the inlet and 
outlet air angles of the blade row. The diffusion factor has been based on 
equivalent diffusion ratio, defined as follows in the current investigation. 


Du = YZL . “102 

Vz2 COS pi 


1.12 + k (i - i*) 143 +0-61 


cos 2 pi 
a 


k 


• AK3 

( 1 . 1 ) 


Here Vz is the axial velocity, P the blade metal angle, i incidence angle, a 
solidity of blade and k, a coefficient that is specific to a blade form. The 
subscripts 1 and 2 refer to blade inlet and outlet conditions, while i* denotes 
design value. The numerical values may be treated as constants pertaining to a 
typical blade in a class of blades designed with a particular level of technology. 
The parameter AK3 then is of the nature of an adjustable factor by means of 
which the equivalent diffusion ratio for another blade (of the same class) can be 
obtained. Thus, at design, noting that (i - i*) is equal to zero, one can find D e q 
and adjust it with respect to a measured or otherwise determined value by means 
of a suitable value for the AK3 factor. The expression (1.1) itself can then be 
used for determining D e q corresponding to an off-design value for the angle of 

incidence. 

(2) Deviation angle 

The deviation angle is a measure of the departure of the air exit angle in a 
given blade row compared to the metal exit angle of a blade for given entry 
conditions of air flow including the incidence angle. Its value may be positive or 
negative, a large positive value indicating separation of flow over the blade 
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surface somewhere along the chord. While the metal turning angle over a blade 
is given by the metal inlet and outlet angles, the deflection of air over the blade is 
obtained by including the incidence and the deviation angles with the inlet and the 
outlet metal angles, respectively. The deflection of air is a parameter in the 
determination of the change in angular momentum and hence of the work 
absorbed by the air in a rotor. 

The deviation angle is defined as follows in the current investigation. 

5 = 8* + 6.40 - 9.45 (Mi - 0.6)(Deq - D^) • AK1 (L2 ) 

where Mi is the Mach number at inlet to blade row and AK1 is an adjustable 
factor with the same status as the factor AK3. The expression (1.2) contains the 
factor (Mi - 0.6) and therefore needs further adjustment when the inlet Mach 

number is over 0.6. 

(3) Nondimensional Wake Momentum Thickness 

The total pressure increases across a compressor stage, by design, 
due to input of work. However there are losses due to various causes, including 
frictional losses, over the surface of a blade. The losses are denoted by a total 
pressure loss factor and defined, according to Reference 9, as follows. 



Here (0 / c) denotes the wake momentum thickness, 0, non-dimensionalized with 
respect to chord c. It is calculated in the current investigation according to the 
following rule. 
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(i) For Deq > Deq* '• 


a] = (af +(o.827Mi - 2.69 M? + 2.675M?)(Deq - D^f • AK2 


c) \c 


(1.4) 


(ii) For D e q < Dgq • 


aj = (a) + (o.89 Ml - 8.71 M? + 9.36 M?) (d^ - • AK2) (1 5) 


Here AK2 has the same status as AK1 and AK3 earlier. 



The relations 1.1 - 1.5 involve three adjustable factors AK1, AK2 and AK3. 


In determining those factors at the design point such that the design pressure ratio 
and efficiency are obtained with the aerodynamic parameters and operating 
conditions of the design point, some trial-and-error, as well as iteration with 
respect to blade exit flow conditions, are involved. The reason is that the three 
adjustable factors themselves are not related to one another through an 


independent set of relations. 

Once a set of adjustable factors is chosen such that the design point 
performance is recovered at the design point operating conditions, the relations 
1.1 - 1.5 can be utilized with the chosen values of adjustable factors for 


determining performance under off-design operating conditions. 

It may be noted that in the case of a multistage machine the set of rules 
must be "tuned" with appropriate values of adjustable factors in each stage. As a 
consequence, following stage stacking, it is possible that an accumulation of small 
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differences in performance of individual stages may lead to a noticeable change in 
overall machine performance. Some ingenuity is required in the final evolution 
of rules with acceptable, adjustable factors. 

It may also be pointed out that the rules including die adjustable factors are 
specific to the local station along the span. Thus, referring to figure 2.1, for the 
case of a multistage machine, a separate set of rules for use in various stages is 
required along each of the calculation streamtubes through the machine. In the 
current investigation, the performance of the fan and the compressor has been 
obtained along streamtubes 2 and 5. 

1.2. Extension to the Case of Water Ingestion 

In applying the rules to the case of water ingestion in the WINCOF and the 
WINCOF-I codes, the following methodology is utilized. 

(i) The angle of incidence is based on the velocities calculated with 
respect to air-water mixture, without taking into account any differences in 
velocities of air and small and large droplets. The velocities involved in 
determining the entry angle are the axial velocity of air-water mixture and the 
rotational speed at entry, the latter unaffected by the flowfield. 

(ii) The angle of deviation is based on the blade surface with any film 
formed due to droplet impact. 

(iii) Finally, the inlet Mach number is based on acoustic speed in the local 
state of air-water incidence. 
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APPENDIX n 

ANALYSIS OF CENTRIFUGAL ACTION AND 
FILM FORMATION 

An important change introduced in the WINCOF-I code, relative to the 
WINCOF code (Reference 1), is the method of including film formation and 
motion in the casing clearance continuous with the determination of centrifugal 
action on water in the air-water mixture flowing through the stage of a 
compressor. 

The details on the calculation of motion of water due to centrifugal action 
over blade surfaces and in blade passages has been discussed in detail in reference 

I. The method of calculation remains the same in the WINCOF-I code. In 
essence, centrifugal force is balanced with change in radial momentum of water. 
Gravitational action, drag and buoyancy are neglected since they have been shown 
to have little effect on the motion. 

Details regarding film formation and motion are given in section IL1. The 
manner in which film growth and motion on the one hand and radial motion of 
water due to centrifugal action are coupled in a time-dependent calculation 
scheme is described in section II.2. 

II. 1 Film Formation. Growth and Motion 

The casing film may have its origin on the inlet surface, in which case the 
presence of film becomes a part of the initial conditions of water distribution at 
entry to a fan or compressor stage. 

In any case, a film can become formed in a stage in the casing clearance 
space as water is displaced radially outwards towards the casing wall. The film 
suffers a motion along the gas path due to the shearing action of the adjoining air- 
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water mixture flow. The velocity of film motion is affected, although to an 
unknown extent, due to the continuous deposition of water and the mixing losses 
between the original film (at any instant of time) and the new deposit. 

In order to calculate the growth and motion of film over a finite, but short 
duration of time, the following assumptions are introduced. 

(i) The film as well as the air-water mixture consist of incompressible, 
viscous fluids. 

(ii) The air-water mixture is saturated with water vapor at all times. 

(iii) Surface tension and gravity effects are small and therefore, a 
continuous film exists. 

(iv) the film is in laminar motion while the air-water mixture is in 
turbulent motion. Suitable velocity profiles may be assumed (Reference 7 ) at the 
wall and the interface on the film and the air-water mixture sides. 

Regarding the velocity profiles, they are described as follows. 

laminar profile 


£= Ci + C2 £ + C3 | + Q t OLD 

The velocity boundary conditions for film flow are given by: 
y = o: u = o; and 

y = 81: u = ui 

Here the u represents velocity, ui the velocity at interference, 81 the thickness of 
film and y the coordinate direction normal to the wall. The shear stress at any 
location y is given by x = ji (du / dy) , where |X is the molecular viscosity of 
water. 
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turbulent profile 


Uoo 


% )W 


(n.2) 


The velocity boundary conditions for the air-water mixture flow are given by: 


y = Si : u = ui and 

y = 81 + 5 g : u — Uoo 

Here 5 g is the thickness of the viscous layer of the air-water mixture at the 
interface and Uoo the velocity of the air-water mixture outside the viscous layer. 

The shear stress at the interface may be written as follows. 

T 0 = 8uf ^ 
dx 

where 0 is the momentum thickness of the air-water mixture viscous layer and 
x the direction of air-water mixture and film motion, parallel to the axis of the 
machine. 

The film motion can be analyzed based on a control volume, as shown in 
figure m.l, and invoking the laws of conservation of mass and momentum. The 
film gains momentum transferred across the interface from the air-water mixture 
due to the difference in kinematic viscosity between the two fluids, influenced 
further by the turbulent nature of air-water mixture flow. A gain in momentum 
gives rise to a velocity to the water centrifuged to the casing clearance. By 
conservation of mass, the film becomes thinner at the exit plane of the control 
volume. If the length of the control volume along the direction of flow is set 
equal to a row of blades, one thus obtains the thickness and the velocity of film 
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provided one knows the value of interfacial velocity; that velocity, ui , is 

common to the film and the air-water mixture. 

Referring to figure ELI, the transfer of momentum can be represented by 

the following equation. 

5 g ui 8 g - f 8 g u| dy = 

where 5 g represents the density of the air-water mixture. 

The transfer of mass into the film due to centrifugal action may be written 
as follows. 

ihy + 8 W §1 ui = rhy §2 U 2 (II.4) 

where (u2 - ui) is the gain in velocity due to shearing action. 

Equations (II.3) and (EL4) can be solved by guessing the interfacial velocity 
and using the so-called shooting method. 

n.2. Coupling of Centrifugal Action and Film Motion 

Entry conditions to a stage are assumed to be given with respect to 
distribution of water in the air-water mixture in the span of the blade row and the 
casing clearance space. It is then the objective to determine the distribution of 
air-water mixture at the stage exit, for given values of rotational speed and flow 
coefficient, as a result of centrifugal action and film formation and motion in the 
clearance space. There are two possibilities of interest with respect to the 
condition of film at entry to a stage: its thickness is equal to or less than the 
casing clearance height. In the former case, any centrifuged water must be 




8 W u 2 dy 


(113) 
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splashed back into the span. Thus the result of centrifugal action is a 
redistribution of water in the span. In the second case, with a partially filled 
casing clearance space, centrifugal action adds to the film thickness. In both 
cases, the film, with the given value of initial velocity of motion, undergoes shear 
(locally in the blade row) by the air adjoining air-water mixture and therefore 
attains a new value of velocity of motion at the stage exit. It is the combination 
of film growth and the velocity of its motion due to shear that determines the 
gradual attainment (in time) of quasi-equilibrium condition in the stage, as stated 
earlier. 

The span of blading is divided into a certain number, e.g. ten, of 
streamtubes. The number depends on local droplet size and also, the amount of 
water in air. All streamtubes should be larger than the diameter of drops and the 
largest expected separation distance between drops. The width of the streamtubes 
need not be equal. The casing clearance width is treated as a single streamtube 
whether it is partially or fully filled. 

In order to perform calculation of centrifugal action on water droplets and 
the film growth in successive discrete time intervals, an interval of time must be 
chosen. The mean residence time of air-water mixture in the blade row, equal to 
chord distance divided by axial velocity of flow is generally short enough in 
duration for obtaining an accurate solution. If the changes of film thickness and 
velocity of motion are not gradual during calculation, then the duration of time 
step may be the main cause and needs to be reduced. We denote the calculation 
step interval of time in a stage by AtRS . 

Centrifugal action on water droplets in each streamtube is calculated over 
the interval of time equal to Atxs divided by the number of streamtubes in the 
span. Accounting for accumulation and displacement of water, starting from the 
hub towards the tip of the blade, one can establish the addition to the casing 
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clearance space. The shearing action on the film is calculated over the same 
interval of time. The calculation is repeated the number of times equal to the 
number of streamtubes chosen for the blading; thus, at the end of the set of 
repeated calculations, the total length of time covered is AtRs . The results of the 
calculation yield the distribution of water in the span and the casing clearance 
space as well on the velocity of film motion. 

In a single stage machine , utilizing the initial conditions at entry one can 
calculate the distribution of water in the blade passage taking into account droplet 
impact and rebound processes. The resulting distribution provides the 

distribution in the blade passage for starting the calculation. 

In a single stage machine, AtRs is obviously equal to the length of time to 
cover the length of the machine and, therefore, can be set equal to 5tRs , the 
sweep time over the machine. In view of the fact that in the first sweep quasi- 
equilibrium conditions may not necessarily have been attained, additional sweeps 
are required. At the initiation of each succeeding step, one has to recognize the 
entry conditions and the local conditions in the stage with respect to water 
distribution. Regarding the latter, a part of a streamtube towards the hub may 
have become depleted in the first sweep while all other streamtubes have a small 
change (generally addition) in the amount of water. 

After a certain number of sweeps, one attains the quasi-equilibrium 
condition. For a single stage machine, the performance becomes periodic in time 
as calculations are continued. 

In a multi-stage machine the foregoing procedure is repeated over all of 
the stages in an interval of time that is designated as die sweep time, AtRM . All 
of the data is each of the stages pertaining to water distribution must be saved as 
part of the initial conditions for the next sweep of calculation. 
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The entry conditions to the first stage give rise, in the next sweep to a set 
of new conditions in the first stage blade passage taking into account the result of 
calculation over the first time interval AtRS . The balance of the calculation 
procedure during the second and succeeding sweeps, then, is the same as in the 
first sweep. 

At the end of a certain number of sweeps all of the stages in the multi-stage 
machine can be expected to have attained quasi-equilibrium conditions. At that 
instant of time, a second cycle of changes begins and thus, a periodic change 
occurs in time. 
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APPENDIX IE 


HEAT AND MASS TRANSFER PROCESSES 

Heat and mass transfer between the gaseous phase (air and water vapor) 
and the liquid phase (water, especially in droplet form) is an important process 
occurring in a fan-compressor unit during water ingestion. The temperature of 
air increases along a compressor as work is added in die stages of the unit while 
that of water increases only by orders of degrees. This difference combined with 
increase in pressure of the mixture drives heat and mass transfer processes. The 
basic method of estimating interphase heat and mass transfer in the WINCOF-I 
code remains the same as in the WINCOF code (Reference 1). However there are 
two additional features introduced in the WINCOF-I code as follows. 

(i) Changes due to the large numbers of droplets in a unit volume, that 
affects transfer processes; and 

(ii) duration of heat and mass transfer processes noting the difference in 
residence time of air and droplets in the compressor. 

There is obviously considerable empiricism in accounting for both of those since 
no data or adequately unambiguous methods of analysis are available in either 
case. However, a parametric study of those effects can be revealing, for example 
a study to establish conditions under which most of the liquid water may (or may 
not) become converted into vapor form at the exit plane of a multistage 
compressor. 


62 



III.l Transfer Parameter 

A transfer parameter is introduced as the ratio of actual (heat or mass) 
transfer occurring in a conglomeration of droplets to the ideal value of (heat or 
mass) transfer that can occur between a gas and a single droplet in it. In the 
current case of air-water droplet mixture flow in a compressor, the principal 
impediment to transfer is the presence of a number of droplets in any relevant 
volume. 

A droplet can be expected to be surrounded by some water vapor, and in a 
sudden vaporization of the droplet, there can arise a volume of vapor at the 
ambient pressure equal to that of water in a droplet. The heating and 
vaporization processes in a process spread over a finite period of time are 
complex, especially because of the flow environment prevalent in an axial-flow 
compressor. 

In the WINCOF-I code, a simple parameter has been introduced to correct 
heat and mass transfer rates for a single droplet to obtain values for a 
conglomeration of droplets. It has been named the heat and mass transfer factor. 
The basis for evolving the factor consists in the following assumptions. 

(a) Droplets of diameter, (10 microns), in an air-water mixture with 2 
per cent water by weight (approximately 10 diameters apart in vapor form under 
ambient conditions) may behave as single droplets, with no interference between 
droplets insofar as heat and mass transfer processes are concerned. But, similar 
droplets when they are 2 diameters apart in vapor form preclude heat and mass 
transfer completely. 

(b) The factor varies linearly between a value of zero when heat and mass 
transfer is entirely eliminated and a value of unity when droplets are 10 or more 
vapor-diameters apart. 

m.2. Residence Time of Droplets in a Stage 
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It has been assumed, as stated in Section 2 of the current report, that 
droplets of ingested water may be distinguished as (a) small, moving with 
velocity of air and (b) large, for which a velocity can be assigned independently 
of that of air, at entry to the fan-compressor unit and also at entry to any stage. 
However, this may only be done empirically as a parameter that can be varied in 
a specific calculation 

m.3 Calculation Procedure 

Given the composition of air-water mixture and the mean volumetric 
diameter of droplets at a specific location in the compressor one can calculate the 
number density of droplets per unit volume. Based on the local values of 
pressure and temperature, the vapor-diameter for a droplet can be found and 
thereby, the distance between the droplets in vapor form under local conditions. 
The heat and mass transfer factor then may be defined as follows. 

§ = (III.1) 

10 - 2.0 


Here Ndia represents the number of vapor diameters that the droplets are apart 
locally. 

The relative velocity factor accounting for the difference between air and 
droplet velocity, £ is assigned various values as desired. In determining 
residence time of water in a stage, the ratio of stage width to mean axial flow 
velocity is multiplied by C , , which can have any positive value less than unity. 

Now, heat and mass transfer calculations may be performed at the exit 
plane of each blade row or stage of a fan-compressor unit. It is therefore 
possible to account for £ and £ factors at all such calculation stations as desired. 
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APPENDIX IV 


MODIFICATIONS INTRODUCED IN WINCOF-I CODE 



IV.L Variable Names 

XPHI 

ICHORD 

JEND 

NWATER 

ICASE 

NS 

NSF 

NSLPC 

NSHPC 

RRHUB(I) 

RCA(IJ) 

RBLADE(I) 

STAGRA(IJ) 

SRHUB(I) 

SCA(IJ) 

SBLADE(I) 

STAGSA(IJ) 


inlet flow coefficient 

number of chord-wise steps in centrifugal 
calculation 

number of sweeps through code 

flag for dry/water cases 

inlet water distribution case, usually uniform, 

with = water fraction 

number of stages 

number of fan stages 

number of low pressure compressor stages 
number of high pressure compressor stages 
hub radius at Ith stage rotor 
chord length of Ith stage rotor 
number of blades for Ith stage rotor 
stagger angle for Ith stage rotor 
hub radius at Ith stage stator inlet 
chord length of Ith stage stator 
number of blades for Ith stage stator 
stager angle for Ith stage stator 
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SIGMRA(U) 
SIGMSA(IJ) 
BETSSA (U) 

FNF 

XDIN 

ICENT 

XDDIN 

IICENT 

ICENTV 

TOG 

TOW 

PO 

DIN 

DDIN 

FND 

TOID 

POID 

FNDLPC 

FNDHPC 

XCH4 

RHUMID 

FMWA 


solidity of Ith stage rotor 

solidity of Ith stage stator 

stator outlet absolute flow angle at design pt. for 

Ith stage 

fraction of design corrected rotor speed for a 
particular speed 

initial water content (mass fraction) of small 
droplets 

index for centrifugal calculation of small droplets 
initial water content (mass fraction) of large 
droplets 

index for centrifugal calculation of large droplets 
index for centrifugal calculation of water vapor 
total temperature of gas phase at compressor inlet 
temperature of droplet at compressor inlet 
total pressure at compressor inlet 
initial diameter of small droplets 
initial diameter of large droplets 
rotor corrected speed at design pt. 
compressor inlet temperature at design pt. 
compressor inlet pressure at design pt. 
rotor corrected speed of LPC at design pt. 
rotor corrected speed of HPC at design pt. 
initial methane content (mass fraction) 
initial relative humidity 
molecular weight of air 


66 



PREB 

DLIMT 

GPR(I) 

GAPSa) 

RRTIP(I) 

SRTIP(I) 

IRAD 

RT(I) 

RM(I) 

rhcd 

RFC(I) 

ST(I) 

SM(I) 

SH(I) 

SFC(I) 

BLOCK(I) 

BLOCKS(I) 

IDESIN 


percent of water droplets that rebound after 

impingement on blade surface 

maximum diameter for small droplets 

gap between Ith stage rotor 

gap between rotor blade and stator blade for Ith 

stage 

blade tip radius at Ith stage rotor inlet 
blade tip radius at Ith stage stator outlet 
index for radius at which calculation is earned 
out 

rotor inlet radius at which tip performance 
calculation is carried out 

rotor inlet radius at which mean line performance 

calculation is carried out 

rotor inlet radius at which hub performance 

calculation is carried out 

rotor inlet radius, streamline 3 

stator inlet radius at which tip performance 

calculation is carried out 

stator inlet radius at which mean line performance 

calculation is carried out 

stator inlet radius at which hub performance 

calculation is carried out 

stator inlet radius, streamline 3 

blockage factor for the Ith stage rotor 

blockage factor for Ith stage stator 

index for output 
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IDESPT 

JCENT 

BTIMRA(IJ) 

BT2MRA(I,J) 

BT1MSA(I,J) 

BT2MSA(I,J) 

DSMASS 

BYPASS 

PR12DA(I,J) 

PR13DA(IJ) 

ETARDA(I,J) 

XSAREA(I,J) 

XSAREAS(IJ) 

VAIRWA 

CLEAR(I) 

ZDIS(I) 

PREDES 

TEMDES 

AK1(I) 

AK2(I) 

AK3(I) 

DVZ1(K,J) 

DVZ2(IJ) 

DVZ3(IJ) 


index for design point output 

index for centrifugal calculation 

blade metal angle at Ith stage rotor inlet 

blade metal angle at Ith stage rotor outlet 

blade metal angle at Ith stage stator inlet 

blade metal angle at Ith stage stator outlet 

streamtube design mass flow for fan 

bypass ratio at design point 

total pressure ratio for the Ith stage rotor at 

design pt. 

total pressure ratio for the Ith stage at design pt. 

adiabatic efficiency for Ith stage 

stream tube area for Ith stage rotor inlet 

stream tube area for Ith stage stator inlet 

initial guess at the film-water interface speed 

clearance of rotor 

length of stage 

pressure on standard day 

temperature on standard day 

constant modifier of deviation angle calculation 

constant modifier of equivalent diffusion ratio 

calculation 

constant modifier of wake momentum thickness 
calculation 

gas axial velocity at rotor inlet at design pt. 
gas axial velocity at rotor outlet at design pt. 
gas axial velocity at stator outlet at design pt. 
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IV.1.1. Heat transfer routine 

The new variables used in the heat transfer routine are as follows: 


Input Variables: 


TGI 

temperature of gaseous phase at stage inlet 

TG3 

temperature of gaseous phase at stage outlet 

TW1 

temperature of droplet at stage inlet 

TW3 

temperature of droplet at stage outlet 

DAVEN2 

droplet nominal diameter at stage inlet 

DAVEN 

droplet nominal diameter at stage outlet 

DELZI 

length of stage 

VZ 

axial velocity 

TTIME 

time of residence of average droplet 

WMASS1 

mass flow of water 

VMASS1 

mass flow of vapor 

AMASS 

mass flow rate of dry air 

CHMASS 

mass flow rate of methane 

CPG 

specific heat constant pressure to gaseous phase 

CPW 

specific heat of water 

RE 

Reynolds number based on relative velocity 
between droplet and gaseous phase 


ii. Output Variables: 

DELTGH temperature drop in gaseous phase due to heat 

transfer between water droplet and gaseous phase 
DELTWH temperature rise in droplet due to heat transfer 

between water droplet and gaseous phase 
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IV. 1.2. Mass transfer routine 

The new variables used in the mass transfer routine are as follows: 


(i) Input Variables: 
XW 

XNP 

VTOT 

DDAVE 

(ii) Output Variables: 

FACMT 


mass fraction of water to air 
number of droplets in the span 
volume of span of stage 
average droplet size in span 

mass transfer factor 


rV.2. Subroutines WICRON.WICDL and WICCEN 
1 . Description: 

Subroutine WICRGN is called at end of rotor aerodynamic 
performance calculations to perform the centrifugal action calculations for all 
(10) streamtubes. The subroutines WICDML and WICCEN are called for 
individual streamtubes, and water is added and subtracted from each streamtube 
corresponding to addition (from a streamtube at lower radius) and depletion (to a 
streamtube of higher radius). 

RT 

RRHUB 

FMMASS 


radius of blade at tip 

radius of blade at hub 

mass of water in casing coming from 

previous stage 
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3. 


Output Variable: 
WATRGN 


mass of water in stream tubes 


4. Usage: 

CALL WICRGN (WATRGN, NREGON, ISTAGE, RT, RRHUB, 
FMMASS) 

IV.3. Subroutine WICFLM 
1 . Description 

Subroutine WICFLM is called after WICRGN and calculates a mean 
velocity of casing water film given the incoming mass of centrifuged water. 


Input Variables: 


VZ 

axial velocity 

FMMASS 

mass of water in casing entering stage 

XWT 

percentage of water in streamtube 

CLEAR 

casing clearance 

RCASE 

radial distance of casing 

NRADS 

flag for redistributed water 

REDMAS 

amount of redistributed water 

RHOM 

density of water 

RT 

radius of blade at tip 

WTMASS 

mass of water in streamtube 

HTOTL 

thickness of casing film 

WATRGN 

mass of water in individual streamtube 

RRHUB 

radius of blade at hub 

NREGON 

number of streamlines 

MMASS 

mass of gas in streamtube 
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3. 

4. 


UMLAST 


mean velocity of film coming from last 


CSTAREA 

FTLMM 

NS 

DELZZ 

Output Variables: 

UIF 

Usage: 


stage 

casing area 

momentum of casing water 
number of stages 
length of stage 

film-gas interface velocity 


CALL WICFLM (VZ, FMMASS XWT, CLEAR, RCASE, NRADS, 
REDMAS, RHOM, RT, WTMASS, HTOTL, WATRGN, RRHUB, 
NREGON, MMASS, UMLAST,CSTAREA, FILMM, UIF, NS, 
DELZZ) 


IV.4 Subroutine Wichet 

1. Description 

Subroutine WICHET is called to determine heat transfer between air 


and water at blade row or stage after the aerodynamic performance calculation 
has been completed. 

2. Input Variables 

These are given in Section IV. 1.1. 

3. Output Variables 

These are given in Section IV. 1.1. 

4. USAGE 

Call WICHET (TGI, TG3, TW1, TW3, DAVEN2, DAVEN, 

DELZI, VZ, TTIME, WMASSI, VMASS1, AMASS, CHMASS, 

CPG, CPW, DELTGH, DELTWH, RE). 
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IV.5. Subroutine WICMTM 

1 . Description 

Subroutine WICMTM is called to determine mass transfer between 
water and air at blade row or stage after performing the heat transfer 
calculations. 

2. Input Variables 

These are given in Section IV. 1.2. 

3. Output Variables 

These are given in Section IV. 1.2. 

4. Usage 

Call WICMTM (XW, XNP, VTOT, FACMT). 
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APPENDIX V 


ILLUSTRATIVE CASE 

In order to illustrate the application of the WINCOF-I code for 
determining the performance of a multi-stage compressor, a generic two-stage 
compressor has been chosen. 

The details of the WINCOF-I code can be found in Ref. 1 and Appendix IV 
of the current report. 

The three cases chosen for performance estimation are as follows. 

1. Dry case: Basic operation of the compressor with air flow at entry. 

2. Wet case 1 : Operation of the compressor with ingestion of 4.0 per cent 
of large droplets, and output at the end of the first "sweep" through the two 
stages. 

3. Wet case 2: Operation of the compressor as in Wet case 1, and output 
at the end of the first 10 "sweeps" through the two stages. 

The input and the output in each of the three cases are presented in the 
following. 
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Input Data - Dry Case 


0.450 

01 

01 

01 

00 

02000200 

06 . 9507.64 

2 . 5202.453 

2 . 1702.436 

1 . 9092.383 

36 . 0026.00 

49 . 2051.40 

37 . 2038.10 
18 . 1018.50 
07 . 3508 . 1008.73 
2 . 1421 . 8441.617 
1 . 8801 . 6651.484 
1 . 5631 . 4511.326 

36 . 0040 . 0046.00 
37 . 0037 . 8037.90 
28 . 7030 . 9031.80 

20 . 8023 . 9025.40 
1 . 0000.715 
1 . 1060.891 
1 . 5411.260 
0 . 8570 . 8340.853 
0 . 9580 . 9290.940 

I . 1981 . 1251.099 
32 . 0333 . 2032.39 

23 . 9125 . 8126.12 
12 . 5414 . 4616.02 

1.00 

0 . 00010.00010 

0602 . 000597 . 001944.00 

0020 . 00600.0 

08879 . 00602 . 001944 . 0008879 . 008879.0 

0 . 0000000.00000 

028 . 97018 . 00016.00 

050 . 000300.0 
0 . 577000.72800 
0 . 790000.86000 
14 . 46914.237 
14 . 36614 . 11613.913 
2 

14 . 4714.24 

II . 2811.30 
6 . 9477.639 
14 . 4714.24 

14 . 3714.12 
11 . 2111.36 

07 . 3508.10 

14 . 3714.12 
0 . 9850.950 
0 . 9450 . 9550.965 
112 

51 . 0056.15 

42 . 7045.60 
37 . 9031.85 

47 . 4046.65 

31 . 7030.60 
- 1 . 7005.150 

54 . 2254 . 0055.65 
40 . 4043 . 9546.25 

37 . 0041 . 7545.00 



26 . 2226 . 8025.05 
17 . 0017 . 8517.35 
04 . 6006 . 0505.80 

011 . 300000000.0000000 

1 . 2881.232 

1 . 2541.233 
1 . 2771.201 
1 . 2771.222 
1 . 2481.227 
1 . 2621.184 
0 . 8860.943 
0 . 9120.962 
0 . 9110.927 

00 . 375824900.3220294 

00 . 338605600.2975193 

00 . 349854100.2657488 

00 . 365026700 . 313852200.2768841 

00 . 309584900 . 283225100.2416014 

00 . 278055900 . 245356500.2246608 

080.00 

0 . 00100.0010 

04 . 50304.534 

2 

02116 . 200518.7 
2 . 30001.4600 
0 . 09000.0900 
1 . 00001.5000 

588 . 5576 . 6558.3 

653 . 9634 . 6625.3 

613 . 4725 . 9684.6 
526 . 5511 . 4502.8 

632 . 0606 . 5591.6 

722 . 3729 . 5680.0 

560 . 5543.5 

623 . 2610.3 
718 . 2676.2 
0 . 37590.3219 
0 . 36920.3241 
0 . 32900.2545 

1 . 47401 . 44871 . 42461 . 40181 . 38011 . 35951 . 33981 . 32121 . 30341 . 28651 . 27031.2549 

1 . 24031 . 22631 . 21301 . 20031 . 18821 . 17671 . 16561 . 15521 . 14511 . 13561 . 12651.1179 

1 . 10971 . 10181 . 09441 . 08731 . 08061 . 07421 . 06811 . 06241 . 05701 . 05191 . 04711.0425 

1 . 03821 . 03421 . 03051 . 02701 . 02371 . 02071 . 01791 . 01531 . 01291 . 01081 . 00891.0071 

1 . 00561 . 00431 . 00311 . 00211 . 00141 . 00081 . 00031 . 00011 . 00001 . 00011 . 00031.0007 
1 . 00131 . 00201 . 00291 . 00391 . 00511 . 00641 . 00791 . 00951 . 01131 . 01321.0153 
1 . 48431 . 45861 . 43411 . 41091 . 38881 . 36781 . 34791 . 31071 . 27701 . 24641 . 21861.1932 
1 . 17021 . 14921 . 13021 . 11291 . 09721 . 08311 . 07031 . 05891 . 04861 . 03951 . 03151.0245 
1 . 01851 . 01341 . 00921 . 00581 . 00321 . 00141 . 00041 . 00001 . 00031 . 00141 . 00301.0054 
1 . 00831 . 01191.0160 

done 

9.99999 

#eor 

#eof 
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Output - Dry Case 
JSWEEP = 1 


MEAN 

FRACTION OF DESIGN SPEED = 1.00000 
»»»»» LOOP NUMBER 2 ««««« 

IGV AREA= 0.2416014000000 
NHG NUMBER OF STREAMLINES = 10.22089 

1 ***************************** INPUT DATA ************************************’ 

HEAT TRANSFER AFTER ROTOR AND STATOR VERSION 
0 NUMBER OF STAGES” 2 (FAN 0, LPC 2, HPC 0) 

PERFORMANCE AT MEAN 
0 VAPOR IS CENTRIFUGED 

0 LARGE DROPLETS IN ROTOR FREE STREAM ARE NOT CENTRIFUGED 


STAGE 

1 

2 

3 

RRHUB(I) 

6.95 

7.64 


RC (I) 

2.170 

2.436 


RBLADE(I) 

36.00 

26.00 


STAGER (I) 

37.20 

38.10 


STAGES ( I ) 

28.70 

30.90 

31.80 

SRHUB(I) 

7.35 

8.10 

8.73 

SC(I) 

1.880 

1.665 

1.484 

SBLADE ( I ) 

36.00 

40.00 

46.00 

SIGUMR(I) 

1.106 

0.891 


SIGUMS (I) 

0.958 

0.929 

0.940 

BET2SS (I) 

23.91 

25.81 

26.12 

GAPR(I) 

0.577 

0.728 


GAPS (I) 

0.790 

0.860 


RRTIP(I) 

14.47 

14.24 


SRTIP(I) 

14.37 

14.12 

13.91 

RT (I) 

14.47 

14.24 


RM(I) 

11.28 

11.30 


RH(I) 

6.95 

7.64 


ST (I) 

14.37 

14.12 


SM(I) 

11.21 

11.36 


SH (I) 

7.35 

8.10 


BLOCK (I) 

0.985 

0.950 


BLOCKS (I) 

0.945 

0.955 

0.965 

BETIMR(I) 

42.70 

45.60 


BET2MR ( I ) 

31.70 

30.60 


BET IMS (I) 

40.40 

43.95 

46.25 

BET2MS (I) 

17.00 

17.85 

17.35 

PR12D (I) 

1.254 

1.233 


PR13D (I) 

1.248 

1.227 


ETARD(I) 

0.912 

0.962 


DVZl(I) 

653.9 

634.6 

625.3 

DVZ2 (I) 

632.0 

606.5 

591.6 

DVZ3 (I) 

623.2 

610.3 


AK1(I) 

2.300 1 

.460 


AK2 (I) 

0.090 0 

.090 


AK3 (I) 

1.000 1 

.500 



I ***************************** INPUT DATA ***************************** 

0 FNF (FRACTION OF DESIGN CORRECTED SPEED) -1.000 
0 XDIN( INITIAL WATER CONTENT OF SMALL DROPLET) -0 . 000 
XDDIN (INITIAL WATER CONTENT OF LARGE DROPLET) -0 . 000 
RHUMID (INITIAL RELATIVE HUMIDITY)- 0.00 PER CENT 
XCH4 (INITIAL METHANE CONTENT) =0 . 000 
0 TOG (COMPRESSOR INLET TOTAL TEMPRATURE OF GAS)- 602.00 
TOW (COMPRESSOR INLET TEMPERATURE OF DROPLRET) — 597.00 
P0 (COMPRESSOR INLET TOTAL PRESSURE) -1944 . 00 
0 DIN (INITIAL DROPLET DIAMETER OF SMALL DROPLET)- 20.0 
DDIN (INITIAL DROPLET DIAMETER OF LARGE DROPLET)- 600.0 
0 FND (DESIGN ROTATIONAL SPEED)- 8879.0 
nn 



0 DSMASS (DESIGN MASS FLOW RATE) = 11.3000 

0 BYPASS RATIO = 0.0000 

0 COMPRESSOR INLET TOTAL TEMPERATURE (GAS PHASE) 602.00 R 
0 COMPRESSOR INLET TOTAL PRESSURE=1944 . 00 LB/FT**2 

0 PREB (PERCENT OF WATER THAT REBOUND AFTER IMPINGEMENT) = 50.0 PERCENT 
0 ROTOR SPEED— 9565.4 RPM 

0 CORRECTED ROTOR SPEED- 8879.0 RPM( 100.0PER CENT OF DESIGN CORRECTED SPEED) 
0 MOLECULAR WEIGHT OF AIR- 28.9700 

0 MAXIMUM DIAMETER OF SMALL DROPLETS- 300.0 MICRONS 
0 ROTOR CORRECTED SPEED AT DESIGN POINT- 8879.0 

ROTOR CORRECTED SPEED OF LPC AT DESIGN POINT- 8879.0 
ROTOR CORRECTED SPEED OF HPC AT THE DESIGN POINT- 8879.0 
DESIGN FLOW COEFFICIENT AT INLET =0.7762894645834 

1*************** DESIGN POINT INFORMATION ************ *** 

0 ***** COMPRESSOR INLET ***** 

0 TOTAL TEMPERATURE AT COMPRESSOR INLET- 602.00000 
TOTAL PRESSURE AT COMPRESSOR INLET- 1944.00 
STATIC TEMPERATURE AT COMPRESSOR INLET- 557.29566 
STATIC PRESSURE AT COMPRESSOR INLET- 1483.22 
STATIC DENSITY AT COMPRESSOR INLET- 0.04988 
0 ACOUSTIC SPEED AT COMPRESSOR INLET-1156 ..87477 
AXIAL VELOCITY AT COMPRESSOR INLET- 625.30000 
MACH NUMBER AT COMPRESSOR INLET- -0.63408 
STREAMTUBE AREA AT COMPRESSOR INLET- 0.24160 
FLOW COEFFICIENT AT COMPRESSOR INLET- 0.77629 


1*************** design point information ************ *** 


0 

***** 

STAGE- 

1 ***** 





0 



TOTAL 

TOTAL 

STATIC 

STATIC 

STATIC 




TEMP 

PRESSURE 

TEMP 

PRESSURE 

DENSITY 

0 

ROTOR 

INLET 

602.000 

1944.000 

554.852 

1460.540 

0.049 


ROTOR 

OUTLET 

646.018 

2437.776 

582.467 

1694.070 

0.055 

0 



AXIAL 

ABSOLUTE 

RELATIVE 

TAN COMP 

TAN COMP 




VELOCITY 

VELOCITY 

VELOCITY 

OF ABS VEL 

OF REL VEL 

0 

ROTOR 

INLET 

653.90000 

728.27594 

901.76441 

320.62538 

620.96203 


ROTOR 

OUTLET 

632.00000 

875.68251 

762.34215 

606.13188 

329.61235 

0 



ROTOR 

ABS MACH 

REL MACH 

REL TOTAL 

REL TOTAL 




SPEED 

NUMBER 

NUMBER 

TEMP 

PRESSURE 

0 

ROTOR 

INLET 

941.587 

0.653 

0.781 

622.410 

2185.064 


ROTOR 

OUTLET 

935.744 

0.741 

0.645 

630.648 

5747.723 

0 



ABS FLOW 

REL FLOW 

STREAMTUBE 


FLOW 




ANGLE 

ANGLE 

AREA 

RADIUS 

COEFFICIENT 

0 

ROTOR 

INLET 

26.12000 

43.52001 

0.33861 

11.28000 

0.54140 


ROTOR 

OUTLET 

43.80310 

25.61793 

0.30958 

11.21000 

0.52327 

0 

STAGE 

TOTAL PRESSURE RATIO AT DESIGN 

POINT- 1.24800 



STAGE 

ADIABATIC EFFICIENCY AT DESIGN 

POINT- 0.89109 



ROTOR 

TOTAL PRESSURE RATIO AT DESIGN 

POINT- 1.25400 



ROTOR 

ADIABATIC EFFICIENCY AT DESIGN 

POINT- 0.91200 



ROTOR 

TOTAL TEMPERATURE 

RATIO AT DESIGN POINT- 

1.07312 


1*************** DESIGN point information ************ *** 


0 

***** 

STAGE- 

2 ***** 





0 



TOTAL 

TOTAL 

STATIC 

STATIC 

STATIC 




TEMP 

PRESSURE 

TEMP 

PRESSURE 

DENSITY 

0 

ROTOR 

INLET 

646.018 

2426.112 

606.280 

1941.056 

0.060 


ROTOR 

OUTLET 

687.268 

2991.396 

632.514 

2232.765 

0.066 

0 



AXIAL 

ABSOLUTE 

RELATIVE 

TAN COMP 

TAN COMP 




VELOCITY 

VELOCITY 

VELOCITY 

OF ABS VEL 

OF REL VEL 

0 

ROTOR 

INLET 

634.60000 

694.17156 

916.97320 

281.34853 

661.90836 


ROTOR 

OUTLET 

606.50000 

813.80738 

725.84181 

542.62344 

405.64189 

0 



ROTOR 

ABS MACH 

REL MACH 

REL TOTAL 

REL TOTAL 




SPEED 

NUMBER 

NUMBER 

TEMP 

PRESSURE 

0 

ROTOR 

INLET 

943.257 

0.574 

0.760 

675.988 

2845.160 


ROTOR 

OUTLET 

948.265 

0.661 

0.589 

676.086 

6382.675 

0 



ABS FLOW 

REL FLOW 

STREAMTUBE 


FLOW 




ANGLE 

ANGLE 

AREA 

RADIUS 

COEFFICIENT 

0 

ROTOR 

INLET 

23.91000 

46.20664 

0.29752 

11.30000 

0.53399 
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ROTOR OUTLET 41.81836 33.97680 0.28323 11.36000 

0 STAGE TOTAL PRESSURE RATIO AT DESIGN POINT= 1.22700 

STAGE ADIABATIC EFFICIENCY AT DESIGN POINT= 0.93777 

ROTOR TOTAL PRESSURE RATIO AT DESIGN POINT= 1.23300 

ROTOR ADIABATIC EFFICIENCY AT DESIGN POINT- 0.96200 

ROTOR TOTAL TEMPERATURE RATIO AT DESIGN POINT— 1.06385 
!*************** DESIGN POINT INFORMATION ************ *** 

0********** OVERALL PERFORMANCE AT DESIGN POINT **** ****** 

0 COMPRESSOR INLET TOTAL TEMPERATURE- 602.00 
0 COMPRESSOR INLET TOTAL PRESSURE- 1944.00 
0 CORRECTED MASS FLOW RATE- 135.446 
0 OVERALL TOTAL PRESSURE RATIO- 1.5313 
0 OVERALL TOTAL TEMPERATURE RATIO-1.1416 
0 OVERALL ADIABATIC EFFICIENCY-0.9102 
0 OVERALL TEMPERATURE RISE- 85.268 

0 123456789 

BETISR(I) 43.52 46.21 
BET2SR (I) 25.62 33.98 

AINCSR(I) 0.82 0.61 

ADEVSR ( I ) -6.08 3.38 

BET1SS (I) 43.80 41.82 

BET2SS (I) 23.91 25.81 26.12 

AINCSS(I) 3.40 -2.13 
ADEVSS ( I) 6.91 7.96 

TD (I) 602. 646. 

OMEGR(I) 0.063 0.026 
OMEGS(I) 0.016 0.019 
SITADR(I) .0399 .0174 

SITADS (I) .0120 .0136 

DEQR(I) 1.601 1.587 

DEQS(I) 1.660 1.523 

PHI DESIGN = 0.7762895 

INLET PHI -0.4500000 

1 FAI-0. 4500000 
XDDIN - 0.0000000000000 

NHG MAIN WS(1) TG(1) P(l) RHUMID - 0.00000 602.000001944.00000 0.00001 

NHG MAIN XV(1) XWT(l) XCH4 - 0.00000 0.00000 0.00000 

0 VZ AT IGV INLET = 543.50378 MACH NUMBER = 0.46153 

I XWT WATRGN 

1 0.0000000000000 0.0000000000000 
2 0.0000000000000 0.0000000000000 

3 0.0000000000000 0.0000000000000 

4 0.0000000000000 0.0000000000000 

5 0.0000000000000 0.0000000000000 

6 0.0000000000000 0.0000000000000 

7 0.0000000000000 0.0000000000000 

8 0.0000000000000 0.0000000000000 

9 0.0000000000000 0.0000000000000 

10 0.0000000000000 0.0000000000000 
XV(1) = 1 . 18 658 574928 04E-00 09 
WATRGT - 0.0000000000000 

0 I STAGE— 0 ( IGV) 

0 0.45000 543.50378 1.00000 

CLEAR ( 1 ) - 0.01880000000000 
NHG MAIN START CALCULATIONS FOR STAGE 1 
D1 DWAKEM, W2— 0.0000000000000 514.38689162770 
D2 DWAKEM, RDELV1- 0.0000000000000 10.000000000000 
D3 DWAKEM, RDELV2- 0.0000000000000 0.0000000000000 
B N,DDAVE(N-2) (N)= 3 0.0000000000000 0.0000000000000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 = 0.00000 0.00000 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL- 0.00000 0.00000 0.00000 0.00000 0.00000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 = 0.00000 0.00000 ' 

FILMAS (1) = 0.0000000000000 
UI = 0.0000000000000 


0.51034 


10 11 12 
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HHC = 0.0000000000000 
HTOTL = 0.0000000000000 

******************** INITIAL FLOW COEFFICIENT* 
STAGE TOTAL PRESSURE RATIO- 1.35248 

STAGE TOTAL TEMPERATURE RATIO- 1.10864 

STAGE ADIABATIC EFFICIENCY- 0.82642 

STAGE FLOW COEFFICIENT-O . 305 
AXIAL VELOCITY- 368.37 
ROTOR SPEED-1207.79 


0.450 (STAGE- 1 ) 


***★****★★★*★ 


*ROTOR INLET* *ROTOR OUTLET* *STATOR OUTLET* 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 

S TAT I C TEMP ERATURE ( GAS ) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 
AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 

TANG. COMP. OF ABS . VEL. 
TANG. COMP. OF REL. VEL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 
FLOW COEFFICIENT 
FLOW AREA 

ABSOLUTE FLOW ANGLE 
RELATIVE FLOW ANGLE 
INCIDENCE 
DEVIATION 
DIFFUSION RATIO 
MOMENTUM THICKNESS 
OMEGA (GAS) 

OMEGA (TOTAL) 


1944.0000 
1805.3335 
602.0000 


2636.5495 

2192.1522 

667.4007 


2629.2242 
2450.7589 
667.4007 


589.4268 

633.1831 

654.1579 

0.0574 

0.0649 

0.0702 

0.0574 

0.0649 

0.0702 

368.3678 

336.4388 

339.1701 

389.0255 

641.8765 

399.3158 

895.7528 

514.3869 


941.5874 

935.7442 

943.2569 

125.0838 

546.6391 


816.5036 

389.1051 


1189.7575 

1253.3951 

1253.3859 

0.3270 

0.5205 

0.3186 

0.7529 

0.4171 


0.3050 

0.2786 

0.2854 

0.3386 

0.3280 

0.3007 

18.7555 

58.3890 

31.8559 

65.7174 

49.1517 


23.0174 

17.9890 



17.4517 

14.8559 


3.6654 

2.6688 


0.1314 

0.17567 

0.17567 


0.0192 

0.01648 

0.01648 


0.00000 

0.00000 


D1 DWAKEM, V3— 0.0000000000000 399.31578652366 
D2 DWAKEM, SDELV1- 0.0000000000000 10.000000000000 
D3 DWAKEM, SDELV2- 0.0000000000000 0.0000000000000 
B N,DDAVE(N-2) (N)-4 0.0000000000000 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 = 

NHG DS DL DLGE DSLL AMLGE AMSLL- 0.00000 0.00000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, Dl, D2, D3 - 0 

x ******************** INITIAL FLOW COEFFICIENT- 0.450 (ISTAGE 
0 STAGE TOTAL PRESSURE RATIO- 1.35248 

STAGE TOTAL TEMPERATURE RATIO- 1.10864 

STAGE ADIABATIC EFFICIENCY- 0.82642 

STAGE 1 TOTAL ETA 0.82642DEL T 65.40066 
0PSI- 0.541560 PSI1- 0.447556 LOSS- 0.094004 


0.00000 

0.00000 


0.00000 

0.00000 


00000 

* 1 ) 


0.00000 I 
★★★★*★★*★** 


** STAGE INLET** 


XV- 

xw - 

xww - 

XF - 

XWT— 

XAIR— 

XMETAN— 

XGAS 

WMASS— 

WWMASS— 

FMMASS— 

WTMASS— 

AMASS- 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

7.16060 


** STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

7.16060 


** STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

7.16060 


OA 



CHMASS* 

0.00000 

0.00000 

0.00000 

VMASS* 

0.00000 

0.00000 

0.00000 

GMASS* 

7.16060 

7.16060 

7.16060 

TMASS* 

7.16060 

7.16060 

7.16060 

WS= 

0.00000 

0.00000 

0.00000 

RHOA* 

0.06054 

0.06333 

0.07015 

RHOM* 

0.05453 

0.06332 

0.07013 

RHOG* 

0.05741 

0.06332 

0.07013 

TG= 

602.00000 

667.40066 

667.40066 

TW= 

597.00000 

597.00000 

597.00000 

TWW* 

597.00000 

0.00000 

597.00000 

NHG: TRAGAS, 

TRAWAT = 1.10864 

1.00000 


P= 

1944.00000 

2636.54951 

2629.22421 

TB= 

667.26838 

0.00000 

682.17635 

TDEW* 

272.00755 

274.47523 

274.47523 

WRITING TO 

EXTERNAL PLOT FILES 



CLEAR (2) = 

0.01650000000000 



NHG MAIN START CALCULATIONS FOR 

STAGE 2 


D1 DWAKEM, W2= 0.0000000000000 

581.75628581708 



0.00000 0.00000 
3 DSLL, DLGE, D1 , D2, 


D2 DWAKEM, RDELV1= 0.0000000000000 10.000000000000 
D3 DWAKEM, RDELV2= 0.0000000000000 0.0000000000000 
B N,DDAVE(N-2) (N)=5 0.0000000000000 0.000.0000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 
NHG DS DL DLGE DSLL AMLGE AMSLL* 

NHG WICSIZ WMASSS WMASSL AMING1 2 
FILMAS (2) = 0.0000000000000 
UI = 0.0000000000000 
HHC = 0.0000000000000 
HTOTL = 0.0000000000000 

l ******************** INITIAL FLOW COEFFICIENT 
0 STAGE TOTAL PRESSURE RATIO* 1.18858 

STAGE TOTAL TEMPERATURE RATIO* 1.06192 

STAGE ADIABATIC EFFICIENCY* 0.81180 

0 STAGE FLOW COEFFICIENTS . 289 

AXIAL VELOCITY* 343.18 
ROTOR SPEED=1188 . 42 


0.00000 0.00000 0.00000 

0.00000 0.00000 0.00000 

D3 = 0.00000 0.00000 


0.450 (STAGE* 2 ) 




*ROTOR INLET* *ROTOR OUTLET* * STATOR OUTLET* 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 
STATIC TEMPERATURE (GAS) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 

0 AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 

TANG. COMP. OF ABS. VEL. 
TANG. COMP. OF REL. VEL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 

0 FLOW COEFFICIENT 
FLOW AREA 

0 ABSOLUTE FLOW ANGLE 
RELATIVE FLOW ANGLE 
INCIDENCE 
DEVIATION 
DIFFUSION RATIO 
MOMENTUM THICKNESS 
OMEGA (GAS) 

OMEGA (TOTAL) 


2629.2242 

2446.6287 

667.4007 


3143.2899 

2724.6550 

708.7291 


3125.0534 

2926.7483 

708.7291 


653.8834 

680.5190 

695.6472 

0.0701 

0.0750 

0.0789 

0.0701 

0.0750 

0.0789 

343.1789 

338.8467 

322.4618 

404.0355 

583.7812 

397.5421 

806.6541 

581.7563 


943.2569 

948.2653 

0.0000 

213.2438 

475.3771 


730.0131 

472.8882 


1252.2959 

1291.6200 

1291.6693 

0.3226 

0.4570 

0.3078 

0.6441 

0.4554 


0.2888 

0.2851 

0.2713 

0.2975 

0.2816 

0.2816 

31.8559 

54.5189 

35.7933 

64.8218 

54.3766 


19.2218 

10.5689 



23.7766 

17.9433 


3.9221 

3.4434 


0.0816 

0.0371 


0.13314 

0.13314 


0.04356 

0.04356 


D1 DWAKEM, V3= 0.0000000000000 397.54214855740 
D2 DWAKEM, SDELV1* 0.0000000000000 10.000000000000 
D3 DWAKEM, SDELV2* 0.0000000000000 0.0000000000000 



B N, DDAVE (N-2) (N)=6 0.0000000000000 0.0000000000000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS, Dl, D2,D3 = 0.00000 0.00000 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL- 0.00000 0.00000 0.00000 0.00000 0.00000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, Dl, D2, D3 = 0.00000 0.00000 < 

1 ******************** INITIAL FLOW COEFFICIENT” 0.450 (I STAGE” 2 ) *********** 


0 STAGE TOTAL PRESSURE RATIO” 1.18858 

STAGE TOTAL TEMPERATURE RATIO” 1.06192 

STAGE ADIABATIC EFFICIENCY” 0.81180 


STAGE 2 TOTAL ETA 0.81376DEL T 41.32848 


0PSI= 0.886120 PSI1- 0.721090 

LOSS- 0.165029 


0 

** STAGE INLET** 

**STAGE OUTLET** 

**STAGE OUTLET** 



(BEFORE INTER- 

(AFTER INTER- 



STAGE ADJUST- 

STAGE ADJUST- 



MENT) 

MENT) 

XV- 

0.00000 

0.00000 

0.00000 

XW= 

0.00000 

0.00000 

0.00000 

XWW= 

0.00000 

0.00000 

0.00000 

XF = 

0.00000 

0.00000 

0.00000 

XWT- 

0.00000 

0.00000 

0.00000 

XAIR” 

1.00000 

1.00000 

1.00000 

XMETAN- 

0.00000 

0.00000 

0.00000 

XGAS 

1.00000 

1.00000 

1.00000 

WMASS- 

0.00000 

0.00000 

0.00000 

WWMASS” 

0.00000 

0.00000 

0.00000 

FMMASS” 

0.00000 

0.00000 

0.00000 

WTMASS- 

0.00000 

0.00000 

0.00000 

AMASS- 

7.16060 

7.16060 

7.16060 

CHMASS” 

0.00000 

0.00000 

0.00000 

VMASS- 

0.00000 

0.00000 

0.00000 

GMASS- 

7.16060 

7.16060 

7.16060 

TMASS- 

7.16060 

7.16060 

7.16060 

WS- 

0.00000 

0.00000 

0.00000 

RHOA- 

0.07385 

0.07516 

0.07892 

RHOM- 

0.05453 

0.07515 

0.07890 

RHOG= 

0.07013 

0.07515 

0.07890 

TG= 

667.40066 

708.72914 

708.72914 

TW- 

597.00000 

597.00000 

597.00000 

TWW— 

597.00000 

0.00000 

597.00000 

NHG: TRAGAS, 

TRAWAT = 1.06192 1.00000 


P= 

2629.22421 

3143.28986 

3125.05337 

TB- 

682.17635 

0.00000 

692.03622 

TDEW- 

274.47523 

268.01367 

268.01367 


WRITING TO EXTERNAL PLOT FILES 
1********** OVERALL PERFORMANCE ********* 

0 INITIAL FLOW COEFFICIENTS . 450 

0 CORRECTED SPEED” 8879.0 1.000 FRACTION OF DESIGN CORRECTED SPEED 

0 INITIAL WATER CONTENT (SMALL DROPLET) -0 . 000 
INITIAL WATER CONTENT (LARGE DROPLET) -0 . 000 
INITIAL WATER CONTENT (TOTAL) -0 . 000 
INITIAL RELATIVE HUMIDITY” 0.0 PER CENT 
INITIAL METHANE CONTENTS . 000 
0 COMPRESSOR INLET TOTAL TEMPERATURE” 602.00 
0 COMPRESSOR INLET TOTAL PRESSURE- 1944.00 
0 CORRECTED MASS FLOW RATE OF MIXTURE” 120.29 
0 CORRECTED MASS FLOW RATE OF GAS PHASE 120.29 
0 OVERALL TOTAL PRESSURE RATIO- 1.6075 
0 OVERALL TOTAL TEMPERATURE RATIO-1.1773 
0 OVERALL ADIABATIC EFFICIENCY-0.8152 
0********** PERFORMANCE OF FAN,LPC,HPC ********** 

0 GAS PHASE STAGNATION STAGNATION ADIABATIC 

0 CORRECTED PRESSURE TEMPERATURE EFFICIENCY 

0 MASS FLOW RATIO RATIO 

0 FAN 0.0000 0.0000 0.0000 0.0000 

0 LPC 0.0000 0.0000 0.0000 0.0000 

0 HPC 0.0000 0.0000 0.0000 0.0000 

0PSI- 0.917751 PSI1” 0.748161 LOSS- 0.169590 





0 0.000000 708.7 3125.1 597.0 0.815 

1 = 68 

1 0.000000 708.7 3125.1 597.0 0.815 

NUMBER OF LOOPS = 1 

TOTAL MASS = 1 . 0000000000000E-0008 
0PSI= 0.917751 PSI1= 0.748161 LOSS= 0.169590 
1 = 68 

GEMACH = 0.2731038592486 


0.0000 

0.0000 
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Input Data - Wet Case 1 


0.450 

01 

01 

02 

04 

02000200 

06 . 9507.64 
2 . 5202.453 
2 . 1702.436 
1 . 9092.383 

36 . 0026.00 

49 . 2051.40 

37 . 2038.10 
18 . 1018.50 
07 . 3508 . 1008.73 
2 . 1421 . 8441.617 
1 . 8801 . 6651.484 
1 . 5631 . 4511.326 

36 . 0040 . 0046.00 
37 . 0037 . 8037.90 
28 . 7030 . 9031.80 

20 . 8023 . 9025.40 
1 . 0000.715 
1 . 1060.891 
1 . 5411.260 
0 . 8570 . 8340.853 
0 . 9580 . 9290.940 

I . 1981 . 1251.099 
32 . 0333 . 2032.39 

23 . 9125 . 8126.12 
12 . 5414 . 4616.02 

1.00 

0.00011.00010 

0602 . 000597 . 001944.00 

0020 . 00600.0 

08879 . 00602 . 001944 . 0008879 . 008879.0 

0 . 0000000.00000 

028 . 97018 . 00016.00 

050 . 000300.0 
0 . 577000.72800 
0 . 790000.86000 
14 . 46914.237 
14 . 36614 . 11613.913 
2 

14 . 4714.24 

II . 2811.30 
6 . 9477.639 
14 . 4714.24 

14 . 3714.12 
11 . 2111.36 

07 . 3508.10 

14 . 3714.12 
0 . 9850.950 
0 . 9450 . 9550.965 
112 

51 . 0056.15 

42 . 7045.60 
37 . 9031.85 

47 . 4046.65 

31 . 7030.60 
- 1 . 7005.150 

54 . 2254 . 0055.65 
40 . 4043 . 9546.25 

37 . 0041 . 7545.00 



26.2226.8025.05 
17.0017.8517.35 
04.6006.0505.80 

011.300000000.0000000 

1.2881.232 

1.2541.233 
1.2771.201 
1.2771.222 
1.2481.227 
1.2621.184 
0.8860.943 
0.9120.962 
0.9110.927 

00.375824900.3220294 

00.338605600.2975193 

00.349854100.2657488 

00.365026700.313852200.2768841 

00.309584900.283225100.2416014 

00.278055900.245356500.2246608 

080.00 

0 . 00100.0010 

04.50304.534 

2 

02116.200518.7 
2.30001.4600 
0.09000.0900 
1.00001.5000 

588.5576.6558.3 

653.9634.6625.3 

613.4725.9684.6 
526.5511.4502.8 

632.0606.5591.6 

722.3729.5680.0 

560.5543.5 

623.2610.3 
718.2676.2 
0.37590.3219 
0.36920.3241 
0.32900.2545 

1.47401.44871.42461.40181.38011.35951.33981.32121.30341.28651.27031.2549 

1.24031.22631.21301.20031.18821.17671.16561.15521.14511.13561.12651.1179 

1.10971.10181.09441.08731.08061.07421.06811.06241.05701.05191.04711.0425 

1.03821.03421.03051.02701.02371.02071.01791.01531.01291.01081.00891.0071 

1.00561.00431.00311.00211.00141.00081.00031.00011.00001.00011.00031.0007 
1.00131.00201.00291.00391.00511.00641.00791.00951.01131.01321.0153 
1.48431.45861.43411.41091.38881.36781.34791.31071.27701.24641.21861.1932 
1.17021.14921.13021.11291.09721.08311.07031.05891.04861.03951.03151.0245 
1.01851.01341.00921.00581.00321.00141.00041.00001.00031.00141.00301.0054 

l!00831. 01191. 0160 

done 

9.99999 

#eor 

#eof 
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Output - Wet Case 1 
JSWEEP = 1 
MEAN 

FRACTION OF DESIGN SPEED = 1.00000 
»»»»» LOOP NUMBER 2 ««««« 

IGV AREA= 0.2416014000000 
NHG NUMBER OF STREAMLINES = 10.22089 

i ***************************** input data ************************************ 

HEAT TRANSFER AFTER ROTOR AND STATOR VERSION 
0 NUMBER OF STAGES- 2 (FAN 0, LPC 2, HPC 0) 

PERFORMANCE AT MEAN 
0 VAPOR IS CENTRIFUGED 

0 LARGE DROPLETS IN ROTOR FREE STREAM ARE NOT CENTRIFUGED 


STAGE 

1 

2 

3 

RRHUB(I) 

6.95 

7.64 


RC(I) 

2.170 

2.436 


RBLADE ( I ) 

36.00 

26.00 


STAGER (I) 

37.20 

38.10 


STAGES (I) 

28.70 

30.90 

31.80 

SRHUB(I) 

7.35 

8.10 

8.73 

SC(I) 

1.880 

1.665 

1.484 

SBLADE(I) 

36.00 

40.00 

46.00 

SIGUMR(I) 

1.106 

0.891 


S I GUMS ( I ) 

0.958 

0.929 

0.940 

BET2SS (I) 

23.91 

25.81 

26.12 

GAPR(I) 

0.577 

0.728 


GAPS (I) 

0.790 

0.860 


RRTIP (I) 

14.47 

14.24 


SRTIP (I) 

14.37 

14.12 

13.91 

RT (I) 

14.47 

14.24 


RM(I) 

11.28 

11.30 


RH (I) 

6.95 

7.64 


ST (I) 

14.37 

14.12 


SM(I) 

11.21 

11.36 


SH (I) 

7.35 

8.10 


BLOCK (I) 

0.985 

0.950 


BLOCKS (I) 

0.945 

0.955 

0.965 

BETIMR(I) 

42.70 

45.60 


BET2MR ( I ) 

31.70 

30.60 


BET IMS (I) 

40.40 

43.95 

46.25 

BET2MS ( I ) 

17.00 

17.85 

17.35 

PR12D (I) 

1.254 

1.233 


PR13D (I) 

1.248 

1.227 


ETARD(I) 

0.912 

0.962 


DVZl(I) 

653.9 

634.6 

625.3 

DVZ2 (I) 

632.0 

606.5 

591.6 

DVZ3 (I) 

623.2 

610.3 


AK1 (I) 

2.300 1 

.460 


AK2 (I) 

0.090 0 

.090 


AK3 (I) 

1.000 1 

.500 



l ***************************** INPUT DATA ***************************** 

0 FNF (FRACTION OF DESIGN CORRECTED SPEED) -1.000 
0 XDIN (INITIAL WATER CONTENT OF SMALL DROPLET) =0 . 000 
XDDIN (INITIAL WATER CONTENT OF LARGE DROPLET) -1 . 000 
RHUMID (INITIAL RELATIVE HUMIDITY)- 0.00 PER CENT 
XCH4 (INITIAL METHANE CONTENT) =0.000 
0 TOG (COMPRESSOR INLET TOTAL TEMPRATURE OF GAS)- 602.00 
TOW (COMPRESSOR INLET TEMPERATURE OF DROPLRET) — 597.00 
P0 (COMPRESSOR INLET TOTAL PRESSURE) =1944 . 00 
0 DIN (INITIAL DROPLET DIAMETER OF SMALL DROPLET)- 20.0 
DDIN( INITIAL DROPLET DIAMETER OF LARGE DROPLET)- 600.0 
0 FND (DESIGN ROTATIONAL SPEED)- 8879.0 


r> s' 



DSMASS (DESIGN MASS FLOW RATE) = 11.3000 

BYPASS RATIO = 0.0000 

COMPRESSOR INLET TOTAL TEMPERATURE (GAS PHASE) 602,00 R 

=°-o **» 

SSSStSM 6 ^ 8879.0 RPM ( 100.0PER CENT OF DESIGN CORRECTED SPEED) 

MOLECULAR WEIGHT OF AIR* 28.9700 

MAXIMUM DIAMETER OF SMALL DROPLETS* 300.0 MICRONS 
ROTOR CORRECTED SPEED AT DESIGN POINT* 8879.0 
ROTOR CORRECTED SPEED OF LPC AT DESIGN POINT* 8879.0 
ROTOR CORRECTED SPEED OF HPC AT THE DESIGN POINT* 8879.0 
DESIGN FLOW COEFFICIENT AT INLET =0.7762894645834 

!*************** DESIGN POINT INFORMATION ************ *** 

0 ***** COMPRESSOR INLET ***** 

0 TOTAL TEMPERATURE AT COMPRESSOR INLET* 602.00000 
TOTAL PRESSURE AT COMPRESSOR INLET* 1944.00 
STATIC TEMPERATURE AT COMPRESSOR INLET* 557.29566 
STATIC PRESSURE AT COMPRESSOR INLET* 1483.22 
STATIC DENSITY AT COMPRESSOR INLET* 0.04988 
0 ACOUSTIC SPEED AT COMPRESSOR INLET=1156 .87477 
AXIAL VELOCITY AT COMPRESSOR INLET* 625.30000 
MACH NUMBER AT COMPRESSOR INLET* -0.63408 
STREAMTUBE AREA AT COMPRESSOR INLET* 0.24160 

FLOW COEFFICIENT AT COMPRESSOR INLET* 0.77629 
l*************** DESIGN POINT INFORMATION ************ *** 

** ★★ ★ 

TOTAL STATIC 

PRESSURE TEMP 

1944.000 554.852 

2437.776 582.467 

ABSOLUTE RELATIVE 

VELOCITY VELOCITY 

728.27594 901.76441 

875.68251 762.34215 

ABS MACH REL MACH 

NUMBER NUMBER 

0.653 0.781 

0.741 0.645 

REL FLOW STREAMTUBE 

ANGLE AREA 

43.52001 0.33861 

25.61793 0.30958 


0 


0 


***** STAGE* 1 


ROTOR INLET 
ROTOR OUTLET 


TOTAL 
TEMP 
602.000 
646.018 
AXIAL 
VELOCITY 

ROTOR INLET 653.90000 
ROTOR OUTLET 632.00000 
ROTOR 


ROTOR INLET 
ROTOR OUTLET 


SPEED 
941.587 
935.744 
ABS FLOW 
ANGLE 
26.12000 

43.80310 

STAGE TOTAL PRESSURE RATIO AT DESIGN POINT 
STAGE ADIABATIC EFFICIENCY AT DESIGN POINT 
ROTOR TOTAL PRESSURE RATIO AT DESIGN POINT 
ROTOR ADIABATIC EFFICIENCY AT DESIGN POINT 


ROTOR INLET 
ROTOR OUTLET 


STATIC 
PRESSURE 
1460.540 
1694.070 
TAN COMP 
OF ABS VEL 
320.62538 
606.13188 
REL TOTAL 
TEMP 
622.410 
630.648 

RADIUS 

11.28000 

11.21000 


STATIC 
DENSITY 
0.049 
0.055 
TAN COMP 
OF REL VEL 
620.96203 
329.61235 
REL TOTAL 
PRESSURE 
2185.064 
5747.723 
FLOW 

COEFFICIENT 

0.54140 

0.52327 


1.24800 

0.89109 

1.25400 

0.91200 


ROTOR TOTAL TEMPERATURE RATIO AT DESIGN POINT- 1.07312 
!*************** DESIGN POINT INFORMATION ************ *** 


***** STAGE* 


ROTOR INLET 
ROTOR OUTLET 


2 ★ ★ ★ * ★ 
TOTAL 
TEMP 
646.018 
687.268 
0 AXIAL 

VELOCITY 

0 ROTOR INLET 634.60000 
ROTOR OUTLET 606.50000 
0 ROTOR 

SPEED 
943.257 
948.265 
ABS FLOW 
ANGLE 
23.91000 


ROTOR INLET 
ROTOR OUTLET 


0 ROTOR INLET 


TOTAL 
PRESSURE 
2426.112 
2991.396 
ABSOLUTE 
VELOCITY 
694.17156 
813.80738 
ABS MACH 
NUMBER 
0.574 
0.661 
REL FLOW 
ANGLE 
46.20664 


STATIC 
TEMP 
606.280 
632.514 
RELATIVE 
VELOCITY 
916.97320 
725.84181 
REL MACH 
NUMBER 
0.760 
0.589 
STREAMTUBE 
AREA 
0.29752 


STATIC 
PRESSURE 
1941.056 
2232.765 
TAN COMP 
OF ABS VEL 
281.34853 
542.62344 
REL TOTAL 
TEMP 
675.988 
676.086 

RADIUS 

11.30000 


STATIC 
DENSITY 
0.060 
0.066 
TAN COMP 
OF REL VEL 
661.90836 
405.64189 
REL TOTAL 
PRESSURE 
2845.160 
6382.675 
FLOW 

COEFFICIENT 

0.53399 


07 



ROTOR OUTLET 41.81836 33.97680 0.28323 11.36000 

0 STAGE TOTAL PRESSURE RATIO AT DESIGN POINT= 1.22700 

STAGE ADIABATIC EFFICIENCY AT DESIGN POINT- 0.93777 

ROTOR TOTAL PRESSURE RATIO AT DESIGN POINT= 1.23300 

ROTOR ADIABATIC EFFICIENCY AT DESIGN POINT- 0.96200 

ROTOR TOTAL TEMPERATURE RATIO AT DESIGN POINT= 1.06385 
1 ****** ********* DESIGN POINT INFORMATION ************ *** 

0********** OVERALL PERFORMANCE AT DESIGN POINT **** ****** 

0 COMPRESSOR INLET TOTAL TEMP ERATURE= 602.00 
0 COMPRESSOR INLET TOTAL PRESSURE- 1944.00 
0 CORRECTED MASS FLOW RATE- 135.446 
0 OVERALL TOTAL PRESSURE RATIO- 1.5313 
0 OVERALL TOTAL TEMPERATURE RATIO-1.1416 
0 OVERALL ADIABATIC EFFICIENCY-0.9102 
0 OVERALL TEMPERATURE RISE- 85.268 

0 123456789 

BETISR(I) 43.52 46.21 
BET2SR (I) 25.62 33.98 

AINCSR(I) 0.82 0.61 

ADEVSR(I) -6.08 3.38 

BET1SS (I) 43.80 41.82 

BET2SS (I) 23.91 25.81 26.12 

AINCSS (I) 3.40 -2.13 

ADEVSS (I) 6.91 7.96 

TD (I) 602. 646. 

OMEGR(I) 0.063 0.026 
OMEGS(I) 0.016 0.019 
SITADR(I) .0399 .0174 

SITADS (I) .0120 .0136 

DEQR(I) 1.601 1.587 

DEQS(I) 1.660 1.523 

PHI DESIGN - 0.7762895 

INLET PHI =0.4500000 

1 FAI-0. 4500000 
XDDIN - 0.0000000000000 

NHG MAIN WS(1) TG(1) P (1) RHUMID- 0.00000 602.000001944.00000 0.00001 

NHG MAIN XV (1) XWT(l) XCH4 = 0.00000 0.00000 0.00000 

0 VZ AT IGV INLET = 543.50378 MACH NUMBER - 0.46153 

I XWT WATRGN 

1 0.0000000000000 0.0000000000000 
2 0.0000000000000 0.0000000000000 

3 0.0000000000000 0.0000000000000 

4 0.0000000000000 0.0000000000000 

5 0.0000000000000 0.0000000000000 

6 0.0000000000000 0.0000000000000 

7 0.0000000000000 0.0000000000000 

8 0.0000000000000 0.0000000000000 

9 0.0000000000000 0.0000000000000 

10 0.0000000000000 0.0000000000000 
XV(1) - 1 . 1865857492804E-0009 
WATRGT - 0.0000000000000 

0 ISTAGE— 0 (IGV) 

0 0.45000 543.50378 1.00000 

CLEAR ( 1 ) - 0.01880000000000 
NHG MAIN START CALCULATIONS FOR STAGE 1 
D1 DWAKEM, W2— 0.0000000000000 514.38689162770 
D2 DWAKEM, RDELV1- 0.0000000000000 10.000000000000 
D3 DWAKEM, RDELV2- 0.0000000000000 0.0000000000000 
B N,DDAVE(N-2) (N)=3 0.0000000000000 0.0000000000000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, Dl, D2, D3 - 0.00000 0.00000 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL— 0.00000 0.00000 0.00000 0.00000 0.00000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, Dl, D2, D3 = 0.00000 0.00000 

FILMAS ( 1) - 0.0000000000000 
UI = 0.0000000000000 


0.51034 


10 11 12 


no 



HHC = 0.0000000000000 

HTOTL = 0.0000000000000 , . 

I ******************** INITIAL FLOW COEFFICIENT- 0.450 (STAGE- 1 ) 

0 STAGE TOTAL PRESSURE RATIO- 1.35248 

STAGE TOTAL TEMPERATURE RATIO- 1.10864 

STAGE ADIABATIC EFFICIENCY- 0.82642 

0 STAGE FLOW COEFFICIENT-O . 305 

AXIAL VELOCITY- 368.37 
ROTOR SPEED-1207.79 


★*★****★*★★★*■ 


0 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 
STATIC TEMPERATURE (GAS) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 

0 AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 

TANG. COMP. OF ABS . VEL. 
TANG. COMP. OF REL. VEL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 

0 FLOW COEFFICIENT 
FLOW AREA 

0 ABSOLUTE FLOW ANGLE 
RELATIVE FLOW ANGLE 
INCIDENCE 
DEVIATION 
DIFFUSION RATIO 
MOMENTUM THICKNESS 
OMEGA (GAS) 

OMEGA (TOTAL) 


*ROTOR INLET* *ROTOR OUTLET* * STATOR OUTLET* 
1944.0000 2636.5495 2629.2242 

2192.1522 
667.4007 
633.1831 
0.0649 
0.0649 
336.4388 
641.8765 
514.3869 
935.7442 
. 546.6391 
389.1051 
1253.3951 
0.5205 
0.4171 
0.2786 
0.3280 
58.3890 
49.1517 
17.9890 
17.4517 
3.6654 
0.1314 
0.17567 
0.17567 


1805.3335 

602.0000 

589.4268 

0.0574 

0.0574 

368.3678 

389.0255 

895.7528 

941.5874 

125.0838 

816.5036 

1189.7575 

0.3270 

0.7529 

0.3050 

0.3386 

18.7555 

65.7174 

23.0174 


2450.7589 

667.4007 

654.1579 

0.0702 

0.0702 

339.1701 

399.3158 

943.2569 


1253.3859 

0.3186 

0.2854 

0.3007 

31.8559 


14.8559 

2.6688 

0.0192 

0.01648 

0.01648 


D1 DWAKEM, V3— 0.0000000000000 399.31578652366 
D2 DWAKEM, SDELV1- 0.0000000000000 10.000000000000 
D3 DWAKEM, SDELV2- 0.0000000000000 0.0000000000000 
B N, DDAVE (N-2) (N) =4 0.0000000000000 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, Dl, D2, D3 = 
NHG DS DL DLGE DSLL AMLGE AMSLL- 0.00000 0.00000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 
1 ******************** INITIAL FLOW COEFFICIENT— 0.450 
0 STAGE TOTAL PRESSURE RATIO- 1.35248 

STAGE TOTAL TEMPERATURE RATIO- 1.10864 

STAGE ADIABATIC EFFICIENCY- 0.82642 

STAGE 1 TOTAL ETA 0.82642DEL T 65.40066 

0PSI- 0.541560 PSI1- 0.447556 LOSS- 0.094004 


0.00000 0.00000 
0.00000 0.00000 
0.00000 

(I STAGE- 1 ) * 


** STAGE INLET** 


XV- 

XW- 

xww- 

XF = 

XWT— 

XAIR— 

XMETAN— 

XGAS 

WMASS— 

WWMASS— 

FMMASS— 

WTMASS— 

AMASS- 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

7.16060 


** STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

7.16060 


**STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

7.16060 


0.00000 
0.00000 
0.00000 
*★★★★**★* * 


on 



CHMASS= 

0.00000 

0.00000 

0.00000 

VMASS= 

0.00000 

0.00000 

0.00000 

GMASS= 

7.16060 

7.16060 

7.16060 

TMASS= 

7.16060 

7.16060 

7.16060 

WS= 

0.00000 

0.00000 

0.00000 

RHOA- 

0.06054 

0.06333 

0.07015 

RHOM= 

0.05453 

0.06332 

0.07013 

RHOG= 

0.05741 

0.06332 

0.07013 

TG= 

602.00000 

667.40066 

667.40066 

TW= 

597.00000 

597.00000 

597.00000 

TWW= 

597.00000 

0.00000 

597.00000 

TRAGAS , 

TRAWAT = 1.10864 

1.00000 


P= 

1944.00000 

2636.54951 

2629.22421 

TB*= 

6 67.26838 

0.00000 

682.17635 

TDEW= 

272.00755 

274.47523 

274.47523 


0.00000 0.00000 

3 DSLL, DLGE, Dl, D2, D3 


0 


WRITING TO EXTERNAL PLOT FILES 
CLEAR { 2 ) = 0.01650000000000 
NHG MAIN START CALCULATIONS FOR STAGE 2 
Dl DWAKEM, W2= 0.0000000000000 581.75628581708 
D2 DWAKEM, RDELVl* 0.0000000000000 10.000000000000 
D3 DWAKEM, RDELV2- 0.0000000000000 0.0000000000000 
B N, DDAVE (N-2) (N)=5 0.0000000000000 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 
NHG DS DL DLGE DSLL AMLGE AMSLL* 

NHG WICSIZ WMASSS WMASSL AMING1 2 
FILMAS (2) * 0.0000000000000 
UI - 0.0000000000000 
HHC = 0.0000000000000 
HTOTL = 0.0000000000000 

******************** INITIAL FLOW COEFFICIENT* 0.450 (STAGE* 2 ) 
STAGE TOTAL PRESSURE RATIO- 1.18858 

STAGE TOTAL TEMPERATURE RATIO* 1.06192 

STAGE ADIABATIC EFFICIENCY- 
STAGE FLOW COEFFICIENT—0 . 289 
AXIAL VELOCITY* 343.18 
ROTOR SPEED=1188 . 42 


0.00000 0.00000 0.00000 

0.00000 0.00000 0.00000 

- 0.00000 0.00000 


************* 


0.81180 


*ROTOR INLET* *ROTOR OUTLET* * STATOR OUTLET* 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 


2629.2242 

2446.6287 

667.4007 


3143.2899 

2724.6550 

708.7291 


OMEGA (GAS) 
OMEGA (TOTAL) 


3125.0534 

2926.7483 

708.7291 


STATIC TEMPERATURE (GAS) 

653.8834 

680.5190 

695.6472 

STATIC DENSITY (GAS) 

0.0701 

0.0750 

0.0789 

STATIC DENSITY (MIXTURE) 

0.0701 

0.0750 

0.0789 

AXIAL VELOCITY 

343.1789 

338.8467 

322.4618 

ABSOLUTE VELOCITY 

404.0355 

583.7812 

397.5421 

RELATIVE VELOCITY 

806.6541 

581.7563 


BLADE SPEED 

943.2569 

948.2653 

0.0000 

TANG. COMP. OF ABS . VEL. 

213.2438 

475.3771 


TANG. COMP. OF REL. VEL. 

730.0131 

472.8882 


ACOUSTIC SPEED 

1252.2959 

1291.6200 

1291.6693 

ABSOLUTE MACH NUMBER 

0.3226 

0.4570 

0.3078 

RELATIVE MACH NUMBER 

0.6441 

0.4554 


FLOW COEFFICIENT 

0.2888 

0.2851 

0.2713 

FLOW AREA 

0.2975 

0.2816 

0.2816 

ABSOLUTE FLOW ANGLE 

31.8559 

54.5189 

35.7933 

RELATIVE FLOW ANGLE 

64.8218 

54.3766 


INCIDENCE 

19.2218 

10.5689 


DEVIATION 


23.7766 

17.9433 

DIFFUSION RATIO 


3.9221 

3.4434 

MOMENTUM THICKNESS 


0.0816 

0.0371 


0.13314 

0.13314 


0.04356 


Dl DWAKEM, V3* 0.0000000000000 397.54214855740 
D2 DWAKEM, SDELV1* 0.0000000000000 10.000000000000 
D3 DWAKEM, SDELV2* 0.0000000000000 0.0000000000000 


r> a 



B N, DDAVE (N-2) (N)=6 0.0000000000000 0.0000000000000 nnnnn « nnnnn 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 = 0.00000 n'nnnnn 
NHG DS DL DLGE DSLL AMLGE AMSLL- 0.00000 0.00000 0.00000 0.00000 0.00000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, Dl, D2, D3 “ 0.00000 0.00000 


******************** INITIAL FLOW COEFFICIENT- 0.450 (ISTAGE- 2 ) 




0 STAGE TOTAL PRESSURE RAT 10= 

STAGE TOTAL TEMPERATURE RATIO- 
STAGE ADIABATIC EFFICIENCY- 
STAGE 2 TOTAL ETA 0.81376DEL T 

0PSI- 0.886120 PSI1- 0.721090 LOSS- 0 
0 ** STAGE INLET** 


XV- 0.00000 

XW— 0.00000 

XWW- 0.00000 

XF = 0.00000 

XWT— 0.00000 

XAIR— 1.00000 

XMETAN- 0.00000 

XGAS 1.00000 

WMASS— 0.00000 

WWMASS— 0.00000 

FMMASS— 0.00000 

WTMASS— 0.00000 

AMASS- 7.16060 

CHMASS— 0.00000 

VMASS— 0.00000 

GMASS- 7.16060 

TMASS— 7.16060 

WS- 0.00000 

RHOA- 0.07385 

RHOM- 0.05453 

RHOG- 0.07013 

TG— 667.40066 

TW- 597.00000 

TWW— 597.00000 

NHG: TRAGAS, TRAWAT = 1.06192 

P= 2629.22421 

TB- 682.17635 

TDEW- 274.47523 


18858 
1.06192 
0.81180 
41.32848 
165029 
** STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
1.00000 
0.00000 
1.00000 
0.00000 
0.00000 
0.00000 
0.00000 
7.16060 
0.00000 
0.00000 
7.16060 
7.16060 
0.00000 
0.07516 
0.07515 
0.07515 
708.72914 

597.00000 
0.00000 

1.00000 
3143.28986 

0.00000 

268.01367 


**STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

7.16060 

0.00000 

0.00000 

7.16060 

7.16060 

0.00000 

0.07892 

0.07890 

0.07890 

708.72914 

597.00000 

597.00000 

3125.05337 

692.03622 

268.01367 


WRITING TO EXTERNAL PLOT FILES 
1********** OVERALL PERFORMANCE ********* 

0 INITIAL FLOW COEFFICIENT-O . 450 

0 CORRECTED SPEED- 8879.0 1.000 FRACTION OF DESIGN CORRECTED SPEED 

0 INITIAL WATER CONTENT ( SMALL DROPLET) -0 . 000 
INITIAL WATER CONTENT (LARGE DROPLET) -0 . 000 
INITIAL WATER CONTENT (TOTAL) =0 . 000 
INITIAL RELATIVE HUMIDITY- 0.0 PER CENT 
INITIAL METHANE CONTENT-O.OOO 
COMPRESSOR INLET TOTAL TEMPERATURE- 602.00 
COMPRESSOR INLET TOTAL PRESSURE- 1944.00 
CORRECTED MASS FLOW RATE OF MIXTURE- 120.29 
CORRECTED MASS FLOW RATE OF GAS PHASE 120.29 
OVERALL TOTAL PRESSURE RATIO- 1.6075 
0 OVERALL TOTAL TEMPERATURE RAT IO- 1.1 77 3 
0 OVERALL ADIABATIC EFFICIENCY-0.8152 
0********** PERFORMANCE OF FAN,LPC,HPC ********** 


0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


FAN 

LPC 

HPC 


GAS PHASE 
CORRECTED 
MASS FLOW 
0.0000 
0.0000 
0.0000 


STAGNATION 
PRESSURE 
RATIO 
0.0000 
0.0000 
0.0000 


STAGNATION 

TEMPERATURE 

RATIO 

0.0000 

0.0000 

0.0000 


ADIABATIC 

EFFICIENCY 

0.0000 

0.0000 

0.0000 


0PSI- 0.917751 PSI1- 0.748161 LOSS- 0.169590 


m 



3125.1 597.0 0.815 0.0000 


0 0.000000 708.7 

1 = 68 

PHI DESIGN = 0.7762895 


1 FAI=0. 4500000 

XDDIN = 0.04000000000000 

NHG MAIN WS (1) TG(1) P(l) RHUMID - 0.01187 602.000001944.00000 100.00000 

NHG MAIN XV (1) XWT(l) XCH4 = 0.01126 0.04000 0.00000 

0 VZ AT IGV INLET = 543.50378 MACH NUMBER - 0.46995 

I XWT WATRGN 

1 0.04000000000000 0.3016412676757 

2 0.04000000000000 0.3016412676757 

3 0.04000000000000 0.3016412676757 

4 0.04000000000000 0.3016412676757 

5 0.04000000000000 0.3016412676757 

6 0.04000000000000 0.3016412676757 

7 0.04000000000000 0.3016412676757 

8 0.04000000000000 0.3016412676757 

9 0.04000000000000 0.3016412676757 

10 0.04000000000000 0.3016412676757 
XV(1) = 0.01125764163510 

WATRGT = 3.0164126767578 

DO DWAKEM, W2* 600.00000000000 543.50378458031 
0 ISTAGE=0 (IGV) 

0 0.45000 543.50378 0.94874 

CLEAR ( 1 ) = 0.01880000000000 

NHG MAIN START CALCULATIONS FOR STAGE 1 
D1 DWAKEM, W2= 0.0000000000000 513.95749648846 
D2 DWAKEM, RDELV1- 600.00000000000 18.867096152077 
D3 DWAKEM, RDELV2= 436.40800073893 0.0000000000000 
B N,DDAVE(N-2) (N)*3 600.00000000000 0.0000000000000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 - 0.00000 0.14756 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL= 0.00000 0.00000 473.99116 0.00000 0.14756 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 = 0.00000 0.14756 

FILMAS(l) = 0.3114420393203 
UI = 319.86671351019 
HHC = 0.001830860297993 
HTOTL = 3 . 4 42 0173 60228 5E-0 005 
N, DDAVE (N-l) (N)-3 600.00000000000 473.99116130186 
N, DDAVE (N-l) (N)=3 600.00000000000 473.99116130186 

1 ******************** INITIAL FLOW COEFFICIENT* 0.450 (STAGE* 1 ) ************* 

0 STAGE TOTAL PRESSURE RATIO* 1.34569 

STAGE TOTAL TEMPERATURE RATIO* 1.10846 

STAGE ADIABATIC EFFICIENCY* 0.81082 

0 STAGE FLOW COEFFICIENT*0 . 304 

AXIAL VELOCITY* 367.58 
ROTOR SPEED=1207 .79 


0 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 
STATIC TEMPERATURE (GAS) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 
AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 
TANG. COMP. OF ABS. 
TANG. COMP. OF REL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 
FLOW COEFFICIENT 


* ROTOR INLET* 
1944.0000 
1801.3200 
602.0000 


VEL. 

VEL. 


*ROTOR OUTLET* 
2624.4055 
2162.9053 
667.2949 


♦STATOR OUTLET* 
2616.0249 
2431.1640 
667.2949 


589.0906 

632.9487 

654.1192 

0.0569 

0.0636 

0.0692 

0.0593 

0.0662 

0.0720 

367.5771 

339.5677 

340.5856 

388.1904 

646.3270 

400.9942 

895.6727 

513.9575 


941.5874 

935.7442 

943.2569 

124.8153 

549.9385 


816.7721 

385.8057 


1168.8862 

1211.6172 

1230.8882 

0.3321 

0.5334 

0.3258 

0.7663 

0.4242 


0.3043 

0.2811 

0.2866 


Q? 


Q-3. 



FLOW AREA. 0.3386 0.3280 

0 ABSOLUTE FLOW ANGLE 18.7555 58.3061 

RELATIVE FLOW ANGLE 65.7705 48.6473 

INCIDENCE 23.0705 17.9061 

DEVIATION 16.9473 

DIFFUSION RATIO 3.6741 

MOMENTUM THICKNESS 0.1392 

OMEGA (GAS) °- 17 ^6 

OMEGA (TOTAL) 0.18603 

D1 DWAKEM, V3= 600.00000000000 400.99420772009 
D2 DWAKEM, SDELV1= 600.00000000000 18.867096152074 
D3 DWAKEM, SDELV2= 399.87962081246 0.0000000000000 
B N, DDAVE (N-2) (N)=4 600.00000000000 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, Dl, D2, D3 - 

H A A n A n A A A A A A 


0.3007 

31.8587 


14.8587 

2.6690 

0.0195 

0.01679 

0.02515 


NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, Dl, D2, D3 - 0.08549 0.12973 0.06207 

NHG DS DL DLGE DSLL AMLGE AMSLL= 0.00000 0.00000 464.52000 4.47702 0.01783 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 - 0.08549 0.12973 < 

N, DDAVE (N-l) (N) =4 473.99116130186 464.51999818913 
N, DDAVE (N-l) (N)=4 473.99116130186 464.51999818913 

XNP, TG (3), P(3) =32007 53 6709. 313 667.21481348588 2616.0249275938 ************ 

x ******************** INITIAL FLOW COEFFICIENT- 0.450 (ISTAGE- 1 ) 

0 STAGE TOTAL PRESSURE RATIO- 1.34569 

STAGE TOTAL TEMPERATURE RATIO- 1.10833 

STAGE ADIABATIC EFFICIENCY- 0.81182 

STAGE 1 TOTAL ETA 0.81182DEL T 65.21481 

0PSI- 0.545538 PSI1- 0.442878 LOSS- 0.102660 


** STAGE INLET** 


XV- 

xw— 

xww- 

XF = 

XWT— 

XAIR— 

XMETAN- 

XGAS 

WMASS- 

WWMASS— 

FMMASS— 

WTMASS— 

AMASS- 

CHMASS— 

VMASS— 

GMASS- 

TMASS— 

WS- 

RHOA- 

RHOM- 

RHOG- 

TG- 

TW- 

TWW— 


0.01126 

0.00000 

0.04000 

0.00140 

0.04000 

0.94874 

0.00000 

0.96000 

0.00000 

0.29512 

0.00000 

0.29512 

6.99988 

0.00000 

0.08306 

7.08293 

7.37806 

0.01187 

0.06054 

0.05619 

0.05691 

602.00000 

597.00000 

597.00000 


** STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.01126 

0.00000 

0.04000 

0.00140 

0.04000 

0.94874 

0.00000 

0.96000 

0.00000 

0.29512 

0.00000 

0.29512 

6.99988 

0.00000 

0.08306 

7.08293 

7.37806 

0.01187 

0.06256 

0.06469 

0.06211 

667.29487 

597.00000 

597.00000 


**STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.01128 

0.00000 

0.03998 

0.00413 

0.03998 

0.94874 

0.00000 

0.96002 

0.00000 

0.29498 

0.31144 

0.29498 

6.99988 

0.00000 

0.08320 

7.08307 

7.37806 

0.01189 

0.06553 

0.07203 

0.06915 

667.21481 

597.00000 

597.44952 


NHG: TRAGAS, TRAWAT - 1.10833 1.00075 n .. a _ 

P= 1944.00000 2624.40546 2616.02493 

TB— 667.26838 0.00000 681.92240 

TDEW- 518.09930 526.98696 526.94014 

WRITING TO EXTERNAL PLOT FILES 
CLEAR ( 2 ) - 0.01650000000000 
NHG MAIN START CALCULATIONS FOR STAGE 2 
Dl DWAKEM, W2- 0.0000000000000 581.54035901182 
D2 DWAKEM, RDELV1- 600.00000000000 18.739086560407 
D3 DWAKEM, RDELV2- 363.78484232851 0.0000000000000 
B N, DDAVE (N-2) (N) -5 473.99116130186 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, Dl, D2, D3 = 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL- 0.00000 0.00000 411.14024 


0.14749 

0.00000 





0.00000 0.14749 < 


NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, Dl, D2,D3 
FILMAS (2) = 0.6208438817958 
UI = 286.46803716888 
HHC = 0.004811309023339 
HTOTL = 7 . 9386598885093E-0005 
N, DDAVE (N-l) (N)=5 464.51999818913 411.14023858397 
N,DDAVE(N-1) (N)— 5 464.51999818913 411.14023858397 


0 


******************** INITIAL FLOW COEFFICIENT* 
STAGE TOTAL PRESSURE RAT 10= 1.18312 

STAGE TOTAL TEMPERATURE RAT IO= 1.06183 

STAGE ADIABATIC EFFICIENCY- 0.78832 

STAGE FLOW COEFFICIENT-O . 290 
AXIAL VELOCITY- 344.27 
ROTOR SPEED-1188.42 


0.450 (STAGE- 2 ) ************* 


*ROTOR INLET* *ROTOR OUTLET* *STATOR OUTLET* 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 
STATIC TEMPERATURE (GAS) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 
AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 

TANG. COMP. OF ABS . VEL. 
TANG. COMP. OF REL. VEL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 
FLOW COEFFICIENT 
FLOW AREA 

ABSOLUTE FLOW ANGLE 
RELATIVE FLOW ANGLE 
INCIDENCE 
DEVIATION 
DIFFUSION RATIO 
MOMENTUM THICKNESS 
OMEGA (GAS) 

OMEGA (TOTAL) 


2616.0249 
2427.0784 
667.2148 


3118.1280 

2683.6525 

708.4709 


3095.0801 

2889.3275 

708.4709 


653.1883 

680.1123 

695.3672 

0.0692 

0.0734 

0.0773 

0.0720 

0.0765 

0.0805 

344.2742 

342.5236 

325.2771 

405.3371 

588.2975 

400.4374 

806.4845 

581.5404 


943.2569 

948.2653 

0.0000 

213.9473 

478.3007 


729.3096 

469.9647 


1230.0358 

1255.1304 

1268.4592 

0.3295 

0.4687 

0.3157 

0.6557 

0.4633 


0.2897 

0.2882 

0.2737 

0.2975 

0.2816 

0.2816 

31.8587 

54.3926 

35.6783 

64.7301 

53.9144 


19.1301 

10.4426 



23.3144 

17.8283 


3.9115 

3.4311 


0.0850 

0.0370 


0.13510 

0.14237 


0.04353 

0.05711 


D1 DWAKEM, V3= 600.00000000000 400.43739382676 
D2 DWAKEM, SDELV1- 600.00000000000 18.739086560403 
D3 DWAKEM, SDELV2- 341.41060899562 0.0000000000000 
B N, DDAVE(N-2) (N)=6 464.51999818913 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, D1,D2, D3 


0.09935 


0.13015 

3.46783 


NHG DS DL DLGE DSLL AMLGE AMSLL- 0.00000 0.00000 419.44960 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 
N, DAVE(N-l) (N)= 6 0.0000000000000 3.4678316406546 
N,DDAVE(N-1) (N)-6 411.14023858397 419.44960239536 
N,DAVE (N-l) (N) -6 0.0000000000000 3.4678316406546 
N, DDAVE (N-l ) (N)=6 411.14023858397 419.44960239536 
N, DAVE (N-l) (N)=6 0.0000000000000 3.4678316406546 
N, DDAVE (N-l) (N)=6 411.14023858397 419.44960239536 
XNP, TG (3) , P (3) =32007536709 .313 708.33570034424 3095.0801326574 
XNP, TG (3) , P (3) -31992431033 .269 708.33570034424 3095.0801326574 
l ******************** INITIAL FLOW COEFFICIENT- 0.450 (ISTAGE- 2 ) 
0 STAGE TOTAL PRESSURE RATIO- 1.18312 

STAGE TOTAL TEMPERATURE RATIO- 1.06163 

STAGE ADIABATIC EFFICIENCY- 0.79091 

STAGE 2 TOTAL ETA 0.79645DEL T 41.12089 
0PSI= 0.891954 PSI1- 0.710396 LOSS- 0.181558 

0 ** STAGE INLET** ** STAGE OUTLET** ** STAGE OUTLET** 

(BEFORE INTER- (AFTER INTER- 
STAGE ADJUST- STAGE ADJUST- 
MENT) MENT) 


0.04814 

0.01734 


0.09935 0.13015 


I 


*★★****★*★* 


A A 



XV* 

0.01128 

0.01128 

xw= 

0.00000 

0.00000 

xww* 

0.03998 

0.03998 

XF = 

0.00413 

0.00413 

XWT* 

0.03998 

0.03998 

XAIR* 

0.94874 

0.94874 

XMETAN* 

0.00000 

0.00000 

XGAS 

0.96002 

0.96002 

WMASS* 

0.00000 

0.00000 

WWMASS* 

0.29498 

0.29498 

FMMASS* 

0.31144 

0.31144 

WTMASS* 

0.29498 

0.29498 

AMASS* 

6.99988 

6.99988 

CHMASS* 

0.00000 

0.00000 

VMASS* 

0.08320 

0.08320 

GMASS* 

7.08307 

7.08307 

TMASS* 

7.37806 

7.37806 

WS= 

0.01189 

0.01189 

RHOA* 

0.07350 

0.07422 

RHOM* 

0.05619 

0.07675 

RHOG* 

0.06915 

0.07368 

TG= 

667.21481 

708.47090 

TW= 

597.00000 

597.00000 

TWW* 

597.44952 

597.44952 

TRAGAS, 

TRAWAT = 1.06163 

1.00127 

P= 

2616.02493 

3118.12802 

TB* 

681.92240 

0.00000 

TDEW* 

526.94014 

525.79479 


0.01133 

0.00000 

0.03993 

0.00823 

0.03993 

0.94874 

0.00000 

0.96007 

0.00000 

0.29458 

0.62084 

0.29458 

6.99988 

0.00000 

0.08360 

7.08348 

7.37806 

0.01194 

0.07377 

0.08067 

0.07745 

708.33570 

597.00000 

598.21106 

3095.08013 

691.50661 

525.71037 


WRITING TO EXTERNAL PLOT FILES 
1********** OVERALL PERFORMANCE ********* 

0 INITIAL FLOW COEFFICIENTS . 450 „„„„„„ 

0 CORRECTED SPEED= 8879.0 1.000 FRACTION OF DESIGN CORRECTED SPEED 

0 INITIAL WATER CONTENT (SMALL DROPLET) =0 . 000 
INITIAL WATER CONTENT (LARGE DROPLET) =0 . 040 
INITIAL WATER CONTENT (TOTAL) *0 . 040 
INITIAL RELATIVE HUMIDITY=100 . 0 PER CENT 
INITIAL METHANE CONTENT =0 .000 
0 COMPRESSOR INLET TOTAL TEMPERATURE* 602.00 
0 COMPRESSOR INLET TOTAL PRESSURE* 1944.00 
0 CORRECTED MASS FLOW RATE OF MIXTURE* 123.94 
0 CORRECTED MASS FLOW RATE OF GAS PHASE 118.99 
0 OVERALL TOTAL PRESSURE RATIO* 1.5921 
0 OVERALL TOTAL TEMPERATURE RATIO=l .1766 
0 OVERALL ADIABATIC EFFICIENCY-0.7979 
0********** PERFORMANCE OF FAN, LPC f HPC ********** 


V 

0 

GAS PHASE 

STAGNATION 

STAGNATION 

ADIABATIC 

0 

CORRECTED 

PRESSURE 

TEMPERATURE 

EFFICIENCY 

0 

MASS FLOW 

RATIO 

RATIO 

0.0000 

0 

FAN 0.0000 

0.0000 

0.0000 

0 

LPC 0.0000 

0.0000 

0.0000 

0.0000 

0 

HPC 0.0000 

0.0000 

0.0000 

0.0000 

0PSI* 0.914367 PSI1- 0.729576 
1= 68 

LOSS* 0.184791 


1 0.001948 

708.3 3095 

.1 598.2 0. 

798 0.0048 


NUMBER OF LOOPS = 1 
TOTAL MASS = 0.3781818163733 
0PSI* 0.914367 PSI1* 0.729576 LOSS* 0.184791 


I* 68 

GEMACH = 0.2731038592486 
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Input Data - Wet Case 2 


0.450 

01 

10 

02 

04 

02000200 

06.9507.64 
2.5202.453 
2.1702.436 
1.9092.383 

36.0026.00 

49.2051.40 

37.2038.10 
18.1018.50 
07.3508.1008.73 
2.1421.8441.617 
1.8801.6651.484 
1.5631.4511.326 

36.0040.0046.00 
37.0037.8037.90 
28.7030.9031.80 

20.8023.9025.40 
1.0000.715 
1.1060.891 
1.5411.260 
0.8570.8340.853 
0.9580.9290.940 

I. 1981.1251.099 
32.0333.2032.39 

23.9125.8126.12 
12.5414.4616.02 

1.00 

0.00011.00010 

0602.000597.001944.00 

0020.00600.0 

08879.00602.001944.0008879.008879.0 

0 . 0000000.00000 

028.97018.00016.00 

050.000300.0 
0.577000.72800 
0.790000.86000 
14.46914.237 
14.36614.11613.913 
2 

14.4714.24 

II. 2811.30 
6.9477.639 
14.4714.24 

14.3714.12 
11.2111.36 

07.3508.10 

14.3714.12 
0.9850.950 
0.9450.9550.965 
112 

51.0056.15 

42.7045.60 
37.9031.85 

47.4046.65 

31.7030.60 
-1.7005.150 

54.2254.0055.65 
40.4043.9546.25 

37.0041.7545.00 


26 . 2226 . 8025.05 
17 . 0017 . 8517.35 
04 . 6006 . 0505.80 

011 . 300000000.0000000 

1 . 2881.232 

1 . 2541.233 
1 . 2771.201 
1 . 2771.222 
1 . 2481.227 
1 . 2621.184 
0 . 8860.943 
0 . 9120.962 
0 . 9110.927 

00 . 375824900.3220294 

00 . 338605600.2975193 

00 . 349854100.2657488 

00 . 365026700 . 313852200.2768841 

00 . 309584900 . 283225100.2416014 

00 . 278055900 . 245356500.2246608 

080.00 

0 . 00100.0010 

04 . 50304.534 

2 

02116 . 200518.7 
2 . 30001.4600 
0 . 09000.0900 
1 . 00001.5000 

588 . 5576 . 6558.3 

653 . 9634 . 6625.3 

613 . 4725 . 9684.6 
526 . 5511 . 4502.8 

632 . 0606 . 5591.6 

722 . 3729 . 5680.0 

560 . 5543.5 

623 . 2610.3 
718 . 2676.2 
0 . 37590.3219 
0 . 36920.3241 
0 . 32900.2545 

1 . 47401 . 44871 . 42461 . 40181 . 38011 . 35951 . 33981 . 32121 . 30341 . 28651 . 27031.2549 

1 . 24031 . 22631 . 21301 . 20031 . 18821 . 17671 . 16561 . 15521 . 14511 . 13561 . 12651.1179 

1 . 10971 . 10181 . 09441 . 08731 . 08061 . 07421 . 06811 . 06241 . 05701 . 05191 . 04711.0425 

1 . 03821 . 03421 . 03051 . 02701 . 02371 . 02071 . 01791 . 01531 . 01291 . 01081 . 00891.0071 

1 . 00561 . 00431 . 00311 . 00211 . 00141 . 00081 . 00031 . 00011 . 00001 . 00011 . 00031.0007 
1 . 00131 . 00201 . 00291 . 00391 . 00511 . 00641 . 00791 . 00951 . 01131 . 01321.0153 
1 . 48431 . 45861 . 43411 . 41091 . 38881 . 36781 . 34791 . 31071 . 27701 . 24641 . 21861.1932 
1 . 17021 . 14921 . 13021 . 11291 . 09721 . 08311 . 07031 . 05891 . 04861 . 03951 . 03151.0245 
1 . 01851 . 01341 . 00921 . 00581 . 00321 . 00141 . 00041 . 00001 . 00031 . 00141 . 00301.0054 

1 . 00831 . 01191.0160 

done 

9.99999 

#eor 

#eof 
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Ouput — Wet Case 2 
JSWEEP =10 
MEAN 

FRACTION OF DESIGN SPEED = 1.00000 
»»»»» LOOP NUMBER 2 ««««« 

IGV AREA= 0.2416014000000 
NHG NUMBER OF STREAMLINES = 10.22089 

1 ***************************** INPUT DATA ************************************' 

HEAT TRANSFER AFTER ROTOR AND STATOR VERSION 
0 NUMBER OF STAGES- 2 (FAN 0, LPC 2, HPC 0) 

PERFORMANCE AT MEAN 
0 VAPOR IS CENTRIFUGED 

0 LARGE DROPLETS IN ROTOR FREE STREAM ARE NOT CENTRIFUGED 


STAGE 

1 

2 

3 

RRHUB(I) 

6.95 

7.64 


RC (I) 

2.170 

2.436 


RBLADE(I) 

36.00 

26.00 


STAGER (I) 

37.20 

38.10 


STAGES (I) 

28.70 

30.90 

31.80 

SRHUB(I) 

7.35 

8.10 

8.73 

SC(I) 

1.880 

1.665 

1.484 

SBLADE(I) 

36.00 

40.00 

46.00 

SIGUMR (I) 

1.106 

0.891 


SIGUMS(I) 

0.958 

0.929 

0.940 

BET2SS (I) 

23.91 

25.81 

26.12 

GAPR(I) 

0.577 

0.728 


GAPS (I) 

0.790 

0.860 


RRTIP (I) 

14.47 

14.24 


SRTIP (I) 

14.37 

14.12 

13.91 

RT (I) 

14.47 

14.24 


RM(I) 

11.28 

11.30 


RH(I) 

6.95 

7.64 


ST (I) 

14.37 

14.12 


SM(I) 

11.21 

11.36 


SH (I) 

7.35 

8.10 


BLOCK (I) 

0.985 

0.950 


BLOCKS (I) 

0.945 

0.955 

0.965 

BET1MR (I) 

42.70 

45.60 


BET2MR (I) 

31.70 

30.60 


BET IMS (I) 

40.40 

43.95 

46.25 

BET2MS (I) 

17.00 

17.85 

17.35 

PR12D (I) 

1.254 

1.233 


PR13D (I) 

1.248 

1.227 


ETARD (I) 

0.912 

0.962 


DVZ1 (I) 

653.9 

634.6 

625.3 

DVZ2 (I) 

632.0 

606.5 

591.6 

DVZ3 (I) 

623.2 

610.3 


AK1 (I) 

2.300 1 

.460 


AK2 (I) 

0.090 0 

.090 


AK3 (I) 

1.000 1 

.500 



1 ***************************** INPUT DATA ***************************** 

0 FNF (FRACTION OF DESIGN CORRECTED SPEED) -1.000 
0 XDIN (INITIAL WATER CONTENT OF SMALL DROPLET) -0 . 000 
XDD IN (INITIAL WATER CONTENT OF LARGE DROPLET) =1 . 000 
RHUMID( INITIAL RELATIVE HUMIDITY)- 0.00 PER CENT 
XCH4 (INITIAL METHANE CONTENT) =0 . 000 
0 TOG (COMPRESSOR INLET TOTAL TEMPRATURE OF GAS)- 602.00 
TOW (COMPRESSOR INLET TEMPERATURE OF DROPLRET) = 597.00 
P0 (COMPRESSOR INLET TOTAL PRESSURE) -1944 . 00 
0 DIN (INITIAL DROPLET DIAMETER OF SMALL DROPLET)- 20.0 
DDIN( INITIAL DROPLET DIAMETER OF LARGE DROPLET)- 600.0 
0 FND (DESIGN ROTATIONAL SPEED)- 8879.0 


rvo 


0 DSMASS (DESIGN MASS FLOW RATE) = 11.3000 

0 BYPASS RATIO = 0.0000 

0 COMPRESSOR INLET TOTAL TEMPERATURE ( GAS PHASE) 602.00 R 
0 COMPRESSOR INLET TOTAL PRESSURE-1944.00 LB/FT**2 

0 PREB (PERCENT OF WATER THAT REBOUND AFTER IMPINGEMENT)- 50.0 PERCENT 
0 ROTOR SPEED= 9565.4 RPM 

0 CORRECTED ROTOR SPEED- 8879.0 RPM( 100.0PER CENT OF DESIGN CORRECTED SPEED) 
0 MOLECULAR WEIGHT OF AIR- 28.9700 

0 MAXIMUM DIAMETER OF SMALL DROPLETS- 300.0 MICRONS 
0 ROTOR CORRECTED SPEED AT DESIGN POINT- 8879.0 

ROTOR CORRECTED SPEED OF LPC AT DESIGN POINT- 8879.0 
ROTOR CORRECTED SPEED OF HPC AT THE DESIGN POINT- 8879.0 
DESIGN FLOW COEFFICIENT AT INLET =0.7762894644891 

l*************** DESIGN POINT INFORMATION ************ *** 

0 ***** COMPRESSOR INLET ***** 

0 TOTAL TEMPERATURE AT COMPRESSOR INLET- 602.00000 
TOTAL PRESSURE AT COMPRESSOR INLET- 1944.00 
STATIC TEMPERATURE AT COMPRESSOR INLET- 557.29566 
STATIC PRESSURE AT COMPRESSOR INLET- 1483.22 
STATIC DENSITY AT COMPRESSOR INLET- 0.04988 
0 ACOUSTIC SPEED AT COMPRESSOR INLET-1156.87477 
AXIAL VELOCITY AT COMPRESSOR INLET- 625.30000 
MACH NUMBER AT COMPRESSOR INLET- -0.63408 
STREAMTUBE AREA AT COMPRESSOR INLET- 0.24160 
FLOW COEFFICIENT AT COMPRESSOR INLET- 0.77629 
1*************** DESIGN POINT INFORMATION ************ *** 

0 ***** STAGE- 1 ***** 

0 TOTAL TOTAL STATIC STATIC STATIC 

TEMP PRESSURE TEMP PRESSURE DENSITY 

0 ROTOR INLET 602.000 1944.000 554.852 1460.540 0.049 

ROTOR OUTLET 646.018 2437.776 582.467 1694.070 0.055 

0 AXIAL ABSOLUTE RELATIVE TAN COMP TAN COMP 

VELOCITY VELOCITY VELOCITY OF ABS VEL OF REL VEL 

0 ROTOR INLET 653.90000 728.27594 901.76441 320.62538 620.96203 

ROTOR OUTLET 632.00000 875.68251 762.34215 606.13188 329.61235 

0 ROTOR ABS MACH REL MACH REL TOTAL REL TOTAL 

SPEED NUMBER NUMBER TEMP PRESSURE 

0 ROTOR INLET 941.587 0.653 0.781 622.410 2185.064 

ROTOR OUTLET 935.744 0.741 0.645 630.648 5747.723 

0 ABS FLOW REL FLOW STREAMTUBE FLOW 

ANGLE ANGLE AREA RADIUS COEFFICIENT 

0 ROTOR INLET 26.12000 43.52001 0.33861 11.28000 0.54140 

ROTOR OUTLET 43.80310 25.61793 0.30958 11.21000 0.52327 

0 STAGE TOTAL PRESSURE RATIO AT DESIGN POINT- 1.24800 

STAGE ADIABATIC EFFICIENCY AT DESIGN POINT- 0.89109 

ROTOR TOTAL PRESSURE RATIO AT DESIGN POINT- 1.25400 

ROTOR ADIABATIC EFFICIENCY AT DESIGN POINT- 0.91200 

ROTOR TOTAL TEMPERATURE RATIO AT DESIGN POINT- 1.07312 
1*************** DESIGN POINT INFORMATION ************ *** 

0 ***** STAGE- 2 ***** 

0 TOTAL TOTAL STATIC STATIC STATIC 

TEMP PRESSURE TEMP PRESSURE DENSITY 

0 ROTOR INLET 646.018 2426.112 606.280 1941.056 0.060 

ROTOR OUT LE T 687.268 2991.396 632.514 2232.765 0.066 

0 AXIAL ABSOLUTE RELATIVE TAN COMP TAN COMP 

VELOCITY VELOCITY VELOCITY OF ABS VEL OF REL VEL 

0 ROTOR INLET 634.60000 694.17156 916.97320 281.34853 661.90836 

ROTOR OUTLET 606.50000 813.80738 725.84181 542.62344 405.64189 

0 ROTOR ABS MACH REL MACH REL TOTAL REL TOTAL 

SPEED NUMBER NUMBER TEMP PRESSURE 

0 ROTOR INLET 943.257 0.574 0.760 675.988 2845.160 

ROTOR OUTLET 948.265 0.661 0.589 676.086 6382.675 

0 ABS FLOW REL FLOW STREAMTUBE FLOW 

ANGLE ANGLE AREA RADIUS COEFFICIENT 

0 ROTOR INLET 23.91000 46.20664 0.29752 11.30000 0.53399 
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ROTOR OUTLET 41.81836 33.97680 0.28323 11.36000 0.51034 

0 STAGE TOTAL PRESSURE RATIO AT DESIGN POINT- 1.22700 

STAGE ADIABATIC EFFICIENCY AT DESIGN POINT- 0.93777 

ROTOR TOTAL PRESSURE RATIO AT DESIGN POINT- 1.23300 

ROTOR ADIABATIC EFFICIENCY AT DESIGN POINT- 0.96200 

ROTOR TOTAL TEMPERATURE RATIO AT DESIGN POINT- 1.06385 
1*************** DESIGN POINT INFORMATION ************ *** 

0********** OVERALL PERFORMANCE AT DESIGN POINT **** ****** 

0 COMPRESSOR INLET TOTAL TEMPERATURE- 602.00 
0 COMPRESSOR INLET TOTAL PRESSURE- 1944.00 
0 CORRECTED MASS FLOW RATE- 135.446 
0 OVERALL TOTAL PRESSURE RATIO- 1.5313 
0 OVERALL TOTAL TEMPERATURE RATIO-1.1416 
0 OVERALL ADIABATIC EFFICIENCY-0.9102 
0 OVERALL TEMPERATURE RISE- 85.268 

0 123456789 10 11 12 


BETISR(I) 

43.52 

46.21 

BET2SR(I) 

25.62 

33.98 

AINCSR ( I ) 

0.82 

0.61 

ADEVSR(I) 

-6.08 

3.38 

BET1SS (I) 

43.80 

41.82 

BET2SS (I) 

23.91 

25.81 2 

AINCSS (I) 

3.40 

-2.13 

ADEVSS (I) 

6.91 

7.96 

TD (I) 

602. 

646. 

OMEGR(I) 

0.063 

0.026 

OMEGS (I) 

0.016 

0.019 

SIT ADR (I) 

.0399 

.0174 

SITADS(I) 

.0120 

.0136 

DEQR(I) 

1.601 

1.587 

DEQS(I) 

1.660 

1.523 


1 FAI-0. 4500000 


XDDIN - 0.0000000000000 

NHG MAIN WS(1) TG(1) P (1) RHUMID- 0.00000 602.000001944.00000 0.00001 

NHG MAIN XV(1) XWT(l) XCH4 - 0.00000 0.00000 0.00000 

0 VZ AT IGV INLET - 543.50378 MACH NUMBER - 0.46153 

I XWT WATRGN 

1 0.0000000000000 0.0000000000000 

2 0.0000000000000 0.0000000000000 

3 0.0000000000000 0.0000000000000 

4 0.0000000000000 0.0000000000000 

5 0.0000000000000 0.0000000000000 

6 0.0000000000000 0.0000000000000 

7 0.0000000000000 0.0000000000000 

8 0.0000000000000 0.0000000000000 

9 0.0000000000000 0.0000000000000 

10 0.0000000000000 0.0000000000000 
XV(1) = 1 . 1865857492804E-0009 
WATRGT - 0.0000000000000 

0 I STAGE- 0 (IGV) 

0 0.45000 543.50378 1.00000 

CLEAR ( 1 ) - 0.01880000000000 

NHG MAIN START CALCULATIONS FOR STAGE 1 
D1 DWAKEM, W2- 0.0000000000000 514.38689162816 
D2 DWAKEM, RDELV1- 0.0000000000000 10.000000000000 
D3 DWAKEM, RDELV2- 0.0000000000000 0.0000000000000 

B N, DDAVE(N-2) (N)— 3 0.0000000000000 473.99116130186 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 = 0.00000 0.00000 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL— 0.00000 473.99116 0.00000 0.00000 0.00000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, Dl, D2, D3 - 0.00000 0.00000 

FILMAS (1) - 2.8017523868103 
UI = 180.80262815912 
HHC - 0.03580445235613 
HTOTL - 6 . 7312370 42 9536E-00 04 

1 ******************** INITIAL FLOW COEFFICIENT- 0.450 (STAGE- 1 ) ************* 

0 STAGE TOTAL PRESSURE RATIO- 1.35248 


1rtT\ 



STAGE TOTAL TEMPERATURE RATIO= 
STAGE ADIABATIC EFFICIENCY” 
STAGE FLOW COEFFICIENT=0 . 305 
AXIAL VELOCITY” 368.37 
ROTOR SPEED=1207 . 79 


1.10864 

0.82642 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 
STATIC TEMPERATURE (GAS) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 

0 AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 

TANG. COMP. OF ABS . VEL. 
TANG. COMP. OF REL. VEL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 

0 FLOW COEFFICIENT 
FLOW AREA 

0 ABSOLUTE FLOW ANGLE 
RELATIVE FLOW ANGLE 
INCIDENCE 
DEVIATION 
DIFFUSION RATIO 
MOMENTUM THICKNESS 
OMEGA (GAS) 

OMEGA (TOTAL) 


*ROTOR INLET* *ROTOR OUTLET* * STATOR OUTLET* 
1944.0000 2636.5495 2629.2242 

1805.3335 2192.1522 

602.0000 667.4007 


2450.7589 

667.4007 


589.4268 

633.1831 

654.1579 

0.0574 

0.0649 

0.0702 

0.0574 

0.0649 

0.0702 

368.3678 

336.4388 

339.1701 

389.0255 

641.8765 

399.3158 

895.7528 

514.3869 


941.5874 

935.7442 

943.2569 

125.0838 

546.6391 


816.5036 

389.1051 


1189.7575 

1253.3951 

1253.3859 

0.3270 

0.5205 

0.3186 

0.7529 

0.4171 


0.3050 

0.2786 

0.2854 

0.3386 

0.3280 

0.3007 

18.7555 

58.3890 

31.8559 

65.7174 

49.1517 


23.0174 

17.9890 



17.4517 

14.8559 


3.6654 

2.6688 


0.1314 

0.0192 


0.17567 

0.17567 


D1 DWAKEM, V3— 0.0000000000000 399.31578652344 
D2 DWAKEM, SDELV1= 0.0000000000000 10.000000000000 
D3 DWAKEM, SDELV2= 0.0000000000000 0.0000000000000 
B N, DDAVE (N-2) (N)-4 0.0000000000000 464.51999818913 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, D1,D2, D3 = 
NHG DS DL DLGE DSLL AMLGE AMSLL= 0.00000 464.52000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 
l ******************** INITIAL FLOW COEFFICIENT” 0.450 
0 STAGE TOTAL PRESSURE RATIO” 1.35248 

STAGE TOTAL TEMPERATURE RATIO” 1.10864 

STAGE ADIABATIC EFFICIENCY” 0.82642 

STAGE 1 TOTAL ETA 0.82642DEL T 65.40066 

0PSI= 0.541560 PSI1= 0.447556 LOSS” 0.094004 


0.01648 

0.01648 


0.00000 0.00000 0.00000 

0.00000 0.00000 0.00000 

- 0.00000 0.00000 
(ISTAGE” 1 ) *********** 


** STAGE INLET** 


XV- 

XW= 

XWW= 

XF = 

XWT= 

XAIR” 

XMETAN” 

XGAS 

WMASS” 

WWMASS” 

FMMASS— 

WTMASS” 

AMASS- 

CHMASS— 

VMASS- 

GMASS” 

TMASS” 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

2.80175 

0.00000 

7.16060 

0.00000 

0.00000 

7.16060 

7.16060 


** STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
1.00000 
0.00000 
1.00000 
0.00000 
0.00000 
2.80175 
0.00000 
7.16060 
0.00000 
0.00000 
7.16060 
7.16060 
ini 


**STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.00000 

0.00000 

0.00000 

0.03828 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

2.80175 

0.00000 

7.16060 

0.00000 

0.00000 

7.16060 

7.16060 



ws= 

0.00000 

0.00000 

0.00000 

RHOA= 

0.06054 

0.06333 

0.07015 

RHOM= 

0.05453 

0.06332 

0.07013 

RHOG= 

0.05741 

0.06332 

0.07013 

TG= 

602.00000 

667.40066 

667.40066 

TW* 

597.00000 

597.00000 

597.00000 

TWW= 

597.00000 

0.00000 

597.00000 


1.10864 


1.00000 

2636.54951 

0.00000 

274.47523 


2629.22421 

682.17635 

274.47523 


0 


NHG: TRAGAS, TRAWAT = 

P= 1944.00000 

TB= 667.26838 

TDEW= 272.00755 

WRITING TO EXTERNAL PLOT FILES 
CLEAR (2) - 0.01650000000000 
NHG MAIN START CALCULATIONS FOR STAGE 2 
D1 DWAKEM, W2* 0.0000000000000 581.75628581723 
D2 DWAKEM, RDELV1= 0.0000000000000 10.000000000000 
D3 DWAKEM, RDELV2* 0.0000000000000 0.0000000000000 
B N,DDAVE(N-2) (N)-5 0.0000000000000 411.14023858397 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS, D1,D2,D3 = 

NHG DS DL DLGE DSLL AMLGE AMSLL= 

NHG WICSIZ WMASSS WMASSL AMING1 2 
FILMAS (2) = 2.8017523868103 
UI = 180.80262815912 
HHC = 0.05352129167392 
HTOTL = 8 . 8310131261979E-0004 

******************** INITIAL FLOW COEFFICIENT* 0.450 (STAGE* 2 ) 
STAGE TOTAL PRESSURE RATIO* 1.18858 

STAGE TOTAL TEMPERATURE RATIO* 1.06192 

STAGE ADIABATIC EFFICIENCY* 0.81180 

STAGE FLOW COEFFICIENT-O . 289 
AXIAL VELOCITY* 343.18 
ROTOR SPEED=1188 . 42 


0.00000 411.14024 
3 DSLL, DLGE, D1,D2,D3 


0.00000 

0.00000 


0.00000 

0.00000 

0.00000 


0.00000 
0.00000 
0.00000 1 




0 

*ROTOR INLET* 

*ROTOR OUTLET* *STATOR OUTLET* 

TOTAL PRESSURE 

2629.2242 

3143.2899 

3125.0534 

STATIC PRESSURE 

2446.6287 

2724.6550 

2926.7483 

TOTAL TEMPERATURE (GAS) 

667.4007 

708.7291 

708.7291 

STATIC TEMPERATURE (GAS) 

653.8834 

680.5190 

695.6472 

STATIC DENSITY (GAS) 

0.0701 

0.0750 

0.0789 

STATIC DENSITY (MIXTURE) 

0.0701 

0.0750 

0.0789 

0 AXIAL VELOCITY 

343.1789 

338.8467 

322.4618 

ABSOLUTE VELOCITY 

404.0355 

583.7812 

397.5421 

RELATIVE VELOCITY 

806.6541 

581.7563 


BLADE SPEED 

943.2569 

948.2653 

0.0000 

TANG. COMP. OF ABS. VEL 

. 213.2438 

475.3771 


TANG. COMP. OF REL. VEL 

. 730.0131 

472.8882 


ACOUSTIC SPEED 

1252.2959 

1291.6200 

1291.6693 

ABSOLUTE MACH NUMBER 

0.3226 

0.4570 

0.3078 

RELATIVE MACH NUMBER 

0.6441 

0.4554 


0 FLOW COEFFICIENT 

0.2888 

0.2851 

0.2713 

FLOW AREA 

0.2975 

0.2816 

0.2816 

0 ABSOLUTE FLOW ANGLE 

31.8559 

54.5189 

35.7933 

RELATIVE FLOW ANGLE 

64.8218 

54.3766 


INCIDENCE 

19.2218 

10.5689 


DEVIATION 


23.7766 

17.9433 

DIFFUSION RATIO 


3.9221 

3.4434 

MOMENTUM THICKNESS 


0.0816 

0.0371 

OMEGA (GAS) 


0.13314 

0.04356 

OMEGA (TOTAL) 


0.13314 

0.04356 


D1 DWAKEM, V3* 0.0000000000000 397.54214855734 
D2 DWAKEM, SDELV1* 0.0000000000000 10.000000000000 
D3 DWAKEM, SDELV2* 0.0000000000000 0.0000000000000 
B N, DDAVE (N-2) (N)=6 0.0000000000000 419.44960239536 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 = 0.00000 0.00000 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL* 0.00000 419.44960 0.00000 0.00000 0.00000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 = 0.00000 0.00000 ' 


i s\r\ 



I ******************** INITIAL FLOW COEFFICIENT= 0.450 (ISTAGE= 2 ) 
0 STAGE TOTAL PRESSURE RAT 10= 1.18858 

STAGE TOTAL TEMPERATURE RAT 10= 1.06192 

STAGE ADIABATIC EFFICIENCY= 0.81180 

STAGE 2 TOTAL ETA 0.81376DEL T 41.32848 

0PSI= 0.886120 PSI1= 0.721090 LOSS= 0.165029 




** STAGE INLET** 


XV= 0.00000 

XW= 0.00000 

xww= o.ooooo 

XF = 0.03828 

XWT= 0.00000 

XAIR= 1.00000 

XMETAN= 0.00000 

XGAS 1.00000 

WMASS= 0.00000 

WWMASS= 0.00000 

FMMASS= 2.80175 

WTMASS= 0.00000 

AMASS = 7.16060 

CHMASS= 0.00000 

VMASS= 0.00000 

GMASS= 7.16060 

TMASS= 7.16060 

WS= 0.00000 

RHOA= 0.07385 

RHOM= 0.05453 

RHOG= 0.07013 

TG= 667.40066 

TW= 597.00000 

TWW= 597.00000 

TRAGAS , TRAWAT = 1.06192 


** STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.00000 
0.00000 
0.00000 
0.03828 
0.00000 
1.00000 
0.00000 
1.00000 
0.00000 
0.00000 
2.80175 
0.00000 
7.16060 
0.00000 
0.00000 
7.16060 
7.16060 
0.00000 
0.07516 
0.07515 
0.07515 
708.72914 
597.00000 
0.00000 
! 1.00000 
3143.28986 
0.00000 
268.01367 


**STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.00000 

0.00000 

0.00000 

0.03828 

0.00000 

1.00000 

0.00000 

1.00000 

0.00000 

0.00000 

2.80175 

0.00000 

7.16060 

0.00000 

0.00000 

7.16060 

7.16060 

0.00000 

0.07892 

0.07890 

0.07890 

708.72914 

597.00000 

597.00000 


NHG: TRAGAS, TRAWAT = 1.06192 1.00000 

P= 2629.22421 3143.28986 3125.05337 

TB- 682.17635 0.00000 692.03622 

TDEW- 274.47523 268.01367 268.01367 

WRITING TO EXTERNAL PLOT FILES 
1********** OVERALL PERFORMANCE ********* 

0 INITIAL FLOW COEFFICIENTS . 450 

0 CORRECTED SPEED* 8879.0 1.000 FRACTION OF DESIGN CORRECTED SPEED 

0 INITIAL WATER CONTENT (SMALL DROPLET) =0 . 000 
INITIAL WATER CONTENT (LARGE DROPLET) =0 . 000 
INITIAL WATER CONTENT (TOTAL) =0 . 000 
INITIAL RELATIVE HUMIDITY* 0 . 0 PER CENT 
INITIAL METHANE CONTENTS .000 
0 COMPRESSOR INLET TOTAL TEMPERATURE” 602.00 
0 COMPRESSOR INLET TOTAL PRESSURE- 1944.00 
0 CORRECTED MASS FLOW RATE OF MIXTURE- 120.29 
0 CORRECTED MASS FLOW RATE OF GAS PHASE 120.29 
0 OVERALL TOTAL PRESSURE RATIO- 1.6075 
0 OVERALL TOTAL TEMPERATURE RATIO-1.1773 
0 OVERALL ADIABATIC EFFICIENCY-0.8152 
0********** PERFORMANCE OF FAN,LPC,HPC ********** 


STAGNATION STAGNATION 


0 GAS PHASE STAGNATION STAG! 

0 CORRECTED PRESSURE TEMP: 

0 MASS FLOW RATIO RA 

0 FAN 0.0000 0.0000 

0 LPC 0.0000 0.0000 

0 HPC 0.0000 0.0000 

0PSI- 0.917751 PSI1- 0.748161 LOSS- 

I- 68 

1 FAI-0. 4500000 

XDDIN = 0.04000000000000 
NHG MAIN WS(1) TG(1) P(l) RHUMID = 


PRESSURE 

RATIO 

0.0000 

0.0000 

0.0000 


TEMPERATURE 

RATIO 

0.0000 

0.0000 

0.0000 


ADIABATIC 

EFFICIENCY 

0.0000 

0.0000 

0.0000 


0.748161 LOSS- 0.169590 


0.01187 602.000001944.00000 100.00000 


im 



NHG MAIN XV(1) XWT (1) XCH4 = 0.01126 0.04000 0.00000 

0 VZ AT IGV INLET = 543.50378 MACH NUMBER = 0.46995 

I XWT WATRGN 

1 0.04000000000000 0.3016412676757 

2 0.04000000000000 0.3016412676757 

3 0.04000000000000 0.3016412676757 

4 0.04000000000000 0.3016412676757 

5 0.04000000000000 0.3016412676757 

6 0.04000000000000 0.3016412676757 

7 0.04000000000000 0.3016412676757 

8 0.04000000000000 0.3016412676757 

9 0.04000000000000 0.3016412676757 

10 0.04000000000000 0.3016412676757 
XV(1) = 0.01125764163510 

WATRGT = 3.0164126767578 

DO DWAKEM, W2= 600.00000000000 543.50378458031 
0 ISTAGE=0 (IGV) 

0 0.45000 543.50378 0.94874 

CLEAR ( 1 ) = 0.01880000000000 

NHG MAIN START CALCULATIONS FOR STAGE 1 
D1 DWAKEM, W2- 0.0000000000000 513.95749648891 
D2 DWAKEM, RDELV1= 600.00000000000 18.867096152029 
D3 DWAKEM, RDELV2= 436.40800073631 0.0000000000000 
B N, DDAVE (N-2) (N)=3 600.00000000000 0.0000000000000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, Dl, D2, D3 - 0.00000 0.14756 0.00000 

NHG DS DL DLGE DSLL AMLGE AMSLL= 0.00000 0.00000 473.99116 0.00000 0.14756 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, Dl,D2,D3 - 0.00000 0.14756 < 

FILMAS(l) ” 3.1127784005606 
UI = 137.33246072363 
HHC = 0.04420138993732 
HTOTL - 8 . 3098613082176E-0004 
N, DDAVE (N-l) (N)»3 600.00000000000 473.99116129979 
N, DDAVE (N-l) (N)”3 600.00000000000 473.99116129979 

1 ******************** INITIAL FLOW COEFFICIENT” 0.450 (STAGE” 1 ) ************* 

0 STAGE TOTAL PRESSURE RATIO” 1.34569 

STAGE TOTAL TEMPERATURE RATIO” 1.10846 

STAGE ADIABATIC EFFICIENCY” 0.81082 

0 STAGE FLOW COEFFICIENT-O . 304 

AXIAL VELOCITY” 367.58 
ROTOR SPEED=1207 . 79 


0 


0 


0 

0 


*ROTOR INLET* *ROTOR OUTLET* * STATOR OUTLET* 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 
STATIC TEMPERATURE (GAS) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 
AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 

TANG. COMP. OF ABS. VEL. 
TANG. COMP. OF REL. VEL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 
FLOW COEFFICIENT 
FLOW AREA 

ABSOLUTE FLOW ANGLE 
RELATIVE FLOW ANGLE 
INCIDENCE 
DEVIATION 
DIFFUSION RATIO 
MOMENTUM THICKNESS 
OMEGA (GAS) 


1944.0000 

1801.3200 

602.0000 

589.0906 

0.0569 

0.0593 

367.5771 

388.1904 

895.6727 

941.5874 

124.8153 

816.7721 

1168.8862 

0.3321 

0.7663 

0.3043 

0.3386 

18.7555 

65.7705 

23.0705 


104 


2624.4055 

2162.9053 

667.2949 

632.9487 

0.0636 

0.0662 

339.5677 

646.3270 

513.9575 

935.7442 

549.9385 

385.8057 

1211.6172 

0.5334 

0.4242 

0.2811 

0.3280 

58.3061 

48.6473 

17.9061 

16.9473 

3.6741 

0.1392 

0.17986 


2616.0249 

2431.1640 

667.2949 

654.1192 

0.0692 

0.0720 

340.5856 

400.9942 

943.2569 


1230.8882 

0.3258 

0.2866 

0.3007 

31.8587 


14.8587 

2.6690 

0.0195 

0.01679 



0.12973 
4.47702 
08549 


OMEGA (TOTAL) 0.18603 0.02515 

D1 DWAKEM, V3— 600.00000000000 400.99420771989 
D2 DWAKEM, SDELV1= 600.00000000000 18.867096152027 
D3 DWAKEM, SDELV2- 399.87962081002 0.0000000000000 
B N,DDAVE(N-2) (N)=4 600.00000000000 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL, DS, Dl, D2, D3 = 0.08549 

NHG DS DL DLGE DSLL AMLGE AMSLL= 0.00000 0.00000 464.52000 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 = 0 

N, DDAVE (N-l) (N) =4 473.99116129979 464.51999818754 
N, DDAVE (N-l) (N)=4 473.99116129979 464.51999818754 
XNP, TG (3), P(3) =32007 53 67 09. 3 13 667.21481348613 2616.0249275946 
l ******************** INITIAL FLOW COEFFICIENT- 0.450 (ISTAGE- 1 ) 
0 STAGE TOTAL PRESSURE RAT IO= 1.34569 

STAGE TOTAL TEMPERATURE RAT IO= 1.10833 

STAGE ADIABATIC EFFICIENCY- 0.81180 

STAGE 1 TOTAL ETA 0.81180DEL T 65.21481 

0PSI= 0.545629 PSI1= 0.442940 LOSS- 0.102689 


0.06207 
0.01783 
0.12973 i 


0 


** STAGE INLET** 


NHG: 


XV- 0.01126 

XW= 0.00000 

XWW- 0.04000 

XF = 0.00140 

XWT— 0.04000 

XAIR— 0.94874 

XMETAN- 0.00000 

XGAS 0.96000 

WMASS— 0.00000 

WWMASS— 0.29512 

FMMASS— 2.80175 

WTMASS— 0.29512 

AMASS- 6.99988 

CHMASS— 0.00000 

VMASS— 0.08306 

GMASS- 7.08293 

TMASS— 7.37806 

WS- 0.01187 

RHOA- 0.06054 

RHOM- 0.05619 

RHOG- 0.05691 

TG- 602.00000 

TW- 597.00000 

TWW- 597.00000 

TRAGAS, TRAWAT = 1.10833 

P= 1944.00000 

TB- 667.26838 

TDEW- 518.09930 


** STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.01126 
0.00000 
0.04000 
0.00140 
0.04000 
0.94874 
. 0.00000 
0.96000 
0.00000 
0.29512 
2.80175 
0.29512 
6.99988 
0.00000 
0.08306 
7.08293 
7.37806 
0.01187 
0.06256 
0.06469 
0.06211 
667.29487 
597.00000 
597.00000 
1.00075 
2624.40546 
0.00000 
526.98696 


** STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.01146 

0.00000 

0.03979 

0.04128 

0.03979 

0.94874 

0.00000 

0.96021 

0.00000 

0.29360 

3.11278 

0.29360 

6.99988 

0.00000 

0.08458 

7.08445 

7.37806 

0.01208 

0.06552 

0.07200 

0.06914 

667.21481 

597.00000 

597.44952 

2616.02493 

681.92240 

527.42653 


WRITING TO EXTERNAL PLOT FILES 
CLEAR ( 2 ) = 0.01650000000000 
NHG MAIN START CALCULATIONS FOR STAGE 2 
Dl DWAKEM, W2- 0.0000000000000 581.57288050930 
D2 DWAKEM, RDELV1- 600.00000000000 18.736622631016 
D3 DWAKEM, RDELV2- 363.92477277446 0.0000000000000 
B N, DDAVE (N-2) (N)— 5 473.99116129979 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 
NHG DS DL DLGE DSLL AMLGE AMSLL- 
NHG WICSIZ WMASSS WMASSL AMING1 2 
FILMAS (2) = 3.4204451329624 
UI = 169.65473048130 
HHC = 0.04545078593931 
HTOTL = 7.4993796799873E-0004 
N, DDAVE (N-l) (N)=5 464.51999818754 411.26397697453 
N, DDAVE (N-l) (N) =5 464.51999818754 411.26397697453 
1 ******************** INITIAL FLOW COEFFICIENT- 0.450 (STAGE- 2 ) 

1 AC 


**★**★*★*+*• 


0.00000 

0.00000 0.00000 411.26398 

3 DSLL, DLGE, D1,D2,D3 


0.14680 0.00000 

0.00000 0.14680 

0.00000 0.14680 


*★*★**★*★★★** 



0 


STAGE TOTAL PRESSURE RAT 10= 
STAGE TOTAL TEMPERATURE RATIO 
STAGE ADIABATIC EFFICIENCY* 
STAGE FLOW COEFFICIENTS . 290 
AXIAL VELOCITY* 344.39 
ROTOR SPEED=1188 . 42 


1.18315 

1.06181 

0.78865 


*ROTOR INLET* *ROTOR OUTLET* * STATOR OUTLET* 


TOTAL PRESSURE 
STATIC PRESSURE 
TOTAL TEMPERATURE (GAS) 
STATIC TEMPERATURE (GAS) 
STATIC DENSITY (GAS) 
STATIC DENSITY (MIXTURE) 

0 AXIAL VELOCITY 
ABSOLUTE VELOCITY 
RELATIVE VELOCITY 
BLADE SPEED 

TANG . COMP . OF ABS . VEL . 
TANG. COMP. OF REL. VEL. 
ACOUSTIC SPEED 
ABSOLUTE MACH NUMBER 
RELATIVE MACH NUMBER 

0 FLOW COEFFICIENT 
FLOW AREA 

0 ABSOLUTE FLOW ANGLE 
RELATIVE FLOW ANGLE 
INCIDENCE 
DEVIATION 
DIFFUSION RATIO 
MOMENTUM THICKNESS 
OMEGA (GAS) 

OMEGA (TOTAL) 


2616.0249 

2427.0139 

667.2148 


3118.1626 

2683.6917 

708.4574 


3095.1560 

2889.3858 

708.4574 


653.1841 

680.0961 

695.3513 

0.0691 

0.0734 

0.0773 

0.0720 

0.0765 

0.0805 

344.3882 

342.6169 

325.3643 

405.4713 

588.3744 

400.5076 

806.4691 

581.5729 


943.2569 

948.2653 

0.0000 

214.0181 

478.3284 


729.2388 

469.9369 


1230.2106 

1255.2980 

1268.6291 

0.3296 

0.4687 

0.3157 

0.6556 

0.4633 


0.2898 

0.2883 

0.2738 

0.2975 

0.2816 

0.2816 

31.8587 

54.3867 

35.6709 

64.7206 

53.9053 



19.1206 


10.4367 

23.3053 

3.9100 

0.0849 

0.13491 

0.14214 


17.8209 

3.4303 

0.0370 

0.04351 

0.05701 


D1 DWAKEM, V3= 600.00000000000 400.50763293193 
D2 DWAKEM, SDELV1* 600.00000000000 18.736622631013 
D3 DWAKEM, SDELV2- 341.52501255762 0.0000000000000 
B N,DDAVE(N-2) (N)-6 464.51999818754 0.0000000000000 
NHG WICSIZ WMASSS WMASSL AMING1 2 3 DL,DS,D1,D2,D3 ■ 


NHG DS DL DLGE DSLL AMLGE AMSLL* 0.00000 0.00000 419.54008 

NHG WICSIZ WMASSS WMASSL AMING1 2 3 DSLL, DLGE, D1,D2,D3 = 0 

N, DAVE(N-l) (N)*6 0.0000000000000 3.4665741052753 
• N,DDAVE(N-1) (N)=6 411.26397697453 419.54007507640 
N,DAVE(N— 1) (N) *6 0.0000000000000 3.4665741052753 
N, DDAVE (N-l) (N)=6 411.26397697453 419.54007507640 
N,DAVE (N-l) (N) -6 0.0000000000000 3.4665741052753 
N, DDAVE (N-l) (N) =6 411.26397697453 419.54007507640 
XNP, TG (3) , P (3) =32007536709 . 313 708.32307639771 3095.1559980129 
XNP, TG (3),P(3) =31842809389 .321 708.32307639771 3095.1559980129 
1 ******************** INITIAL FLOW COEFFICIENT- 0.450 (ISTAGE 


0.09890 0.12954 0.04' 

L9. 54008 3.46657 0.01' 

= 0.09890 0.12954 


STAGE TOTAL PRESSURE RATIO* 
STAGE TOTAL TEMPERATURE RATIO- 
STAGE ADIABATIC EFFICIENCY- 
STAGE 2 TOTAL ETA 0.79648DEL T 


1.18315 

1.06161 

0.79117 

41.10826 


0PSI* 0.892426 PSI1- 0.710803 LOSS- 0.181623 


** STAGE INLET** 


XV* 

XW- 

xww* 

XF = 

XWT* 

XAIR* 

XMETAN* 

XGAS 


0.01146 

0.00000 

0.03979 

0.04128 

0.03979 

0.94874 

0.00000 

0.96021 


**STAGE OUTLET** 
(BEFORE INTER- 
STAGE ADJUST- 
MENT) 

0.01146 

0.00000 

0.03979 

0.04128 

0.03979 

0.94874 

0.00000 

0.96021 


** STAGE OUTLET** 
(AFTER INTER- 
STAGE ADJUST- 
MENT) 

0.01206 

0.00000 

0.03920 

0.04536 

0.03920 

0.94874 

0.00000 

0.96080 


i rue 



WMASS” 

0.00000 

0.00000 

WWMASS” 

0.29360 

0.29360 

FMMASS” 

3.11278 

3.11278 

WTMASS” 

0.29360 

0.29360 

AMASS” 

6.99988 

6.99988 

CHMASS” 

0.00000 

0.00000 

VMASS” 

0.08458 

0.08458 

GMASS” 

7.08445 

7.08445 

TMASS” 

7.37806 

7.37806 

WS= 

0.01208 

0.01208 

RHOA- 

0.07350 

0.07423 

RHOM= 

0.05619 

0.07673 

RHOG= 

0.06914 

0.07368 

TG= 

667.21481 

708.45743 

TW= 

597.00000 

597.00000 

TWW= 

597.44952 

597.44952 

TRAGAS, 

TRAWAT = 1.06161 

1.00127 

P= 

2616.02493 

3118.16261 

TB= 

681.92240 

0.00000 

TDEW= 

527.42653 

526.30787 


0.00000 

0.28921 

3.42045 

0.28921 

6.99988 

0.00000 

0.08897 

7.08885 

7.37806 

0.01271 

0.07374 

0.08057 

0.07742 

708.32308 

597.00000 

598.21027 


3095.15600 

691.50795 

527.65598 


WRITING TO EXTERNAL PLOT FILES 
1********** OVERALL PERFORMANCE ********* . 

0 INITIAL FLOW COEFFICIENTS . 450 

0 CORRECTED SPEED= 8879.0 1.000 FRACTION OF DESIGN CORRECTED SPEED 

0 INITIAL WATER CONTENT (SMALL DROPLET) =0 . 000 
INITIAL WATER CONTENT (LARGE DROPLET) -0 . 040 
INITIAL WATER CONTENT (TOTAL) =0 . 040 


INITIAL RELATIVE HUMIDITY=100 . 0 PER CENT 
INITIAL METHANE CONTENTS .000 
0 COMPRESSOR INLET TOTAL TEMPERATURE= 602.00 
0 COMPRESSOR INLET TOTAL PRESSURE” 1944.00 
0 CORRECTED MASS FLOW RATE OF MIXTURE” 123.94 
0 CORRECTED MASS FLOW RATE OF GAS PHASE 118.99 
0 OVERALL TOTAL PRESSURE RATIO” 1.5922 
0 OVERALL TOTAL TEMPERATURE RAT I 0=1 . 1766 
0 OVERALL ADIABATIC EFFICIENCY-0.7980 
0********** PERFORMANCE OF FAN,LPC,HPC ********** 

0 GAS PHASE STAGNATION STAGNATION ADIABATIC 

0 CORRECTED PRESSURE TEMPERATURE EFFICIENCY 

0 MASS FLOW RATIO RATIO 

0 FAN 0.0000 0.0000 0.0000 0.0000 

0 LPC 0.0000 0.0000 0.0000 0.0000 

0 HPC 0.0000 0.0000 0.0000 0.0000 

0PSI= 0.914259 PSI1= 0.729541 LOSS” 0.184717 
1 = 68 

10 0.019479 708.3 3095.2 598.2 0.798 0.0455 

NUMBER OF LOOPS =10 
TOTAL MASS = 0.3781818163733 
0PSI= 0.914259 PSI1- 0.729541 LOSS- 0.184717 


I” 68 

GEMACH = 0.2731038592486 
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APPENDIX VI 


PERFORMANCE MAPS FOR ENGINE SIMULATION 

The engine simulation code utilized (Reference 2) permits calculation of 
transient performance of a bypass engine for given ambient, flight (altitude and 
Mach number), and power demand conditions with a standard aviation fuel. The 
code is written in a form in which the working fluid (air) is taken through the 
engine from inlet to nozzle while fuel is added in the burner. Thus the 
performance of any of the engine components can be calculated by using a 
substitute method (or sub-routine in the code) when so desired in place of the 
original method or subroutine. However, it is obvious that, for simplicity of 
engine simulation code operation, it is best if the component performance is 
obtained in any substitute subroutine in the same manner as in the original. 

In the current investigation, wherein the interest is in transient engine 
simulation under conditions of water ingestion, one of the primary components 
affected and on which investigations have been carried out specifically and in 
detail is the fan-compressor unit. The performance of that unit can be established 
utilizing the WINCOF-I code. In the following, a brief description is provided 
on the manner in which the output from the WINCOF-I code for the fan- 
compressor unit of the generic engine is processed such that it is compatible with 
the engine simulation code. 

Traditionally, performance data for a compressor are presented in the 
form of standard compressor maps; these are curves showing variation of overall 
pressure ratio and adiabatic efficiency as a function of corrected mass flow, at a 
range of corrected operating speeds. For use in an engine simulation code, these 
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curves would have to be stored in tabular data form, which would require very 
large amounts of memory. Also, when the engine is operating in conditions 
between available data points, interpolation would be required. However, since 
these curves are highly non-linear, any interpolation would very likely be in 
error. Therefore, in order to run an engine simulation code efficiently and 
accurately, a better method of storing compressor performance data is required. 

One such method is based on the loss characteristics of a compressor 
(Reference 10). This method is used to generate a set of parametrized 
compressor performance maps. 

The maps are presented in two ways: a compressor efficiency 

representation and a compressor flow representation. The compressor efficiency 
representation consists of three maps: (a) minimum work coefficient versus 
corrected rotational speed, (b) minimum loss versus speed, and (c) loss minus 
minimum loss versus work coefficient minus minimum-loss work coefficient 
squared, with the sign maintained. These quantities are defined as follows. 

(a) work coefficient. 


V = 


AH 

U 2 / 2G J 


(VI. 1) 


(b) pressure coefficient 

Vi 


AHi 

U 2 / 2GJ 


(VI.2) 


and 

(c) loss 


LOSS = \jr - \|fi 


(VI.3) 


In the three definitions, referring to figure VI. 1, AH represents the work 
actually done on air during compression and AHi the amount of work that, 
ideally, would have resulted in the same pressure ratio across the compressor. 
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They are non-dimensionalized with respect to the kinetic energy of rotation at 
compressor wheel speed. 

The compressor flow representation consists of: (d) minimum loss flow 
coefficient versus corrected rotational speed, and (e) pseudo-Mach number versus 
work coefficient minus minimum loss work coefficient. The latter may be 
explained as follows. When the flow in the compressor is near choking, or a 
maximum value, it is assumed that there is a Mach number of unity somewhere 
along the gas path. From this a critical flow area can be calculated. If this area 
can be assumed to be constant for operation at a given rotational speed, a pseudo- 
Mach number can be defined at each point along the speed line as follows. 


rcicorr 

n^corr t max 




JLt-L 

2 ( 7 - 1 ) 


(VI.4) 


It may be pointed out that the ratio of mass fluxes is nothing but the ratio of flow 
coefficients. 

Examples of the five maps (a) - (e) are given in figure VI.2. 

Utilizing the five curves, one can generate two sets of curves: (1) loss 
minus minimum loss as a function of (V " Yml) * KV - Ym] for various values 
of operating rotor speeds and (2) psendo-Mach number as a function of (\y - 
VML)- In both cases it is interesting that the result is linear with breaks 
corresponding to stalling and choking of compressor. This is the main feature 
that one needed to realize so that the maps may be inserted in a simple table look- 
up in the engine simulation code. 
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Figure VL1. 
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ENTROPY. S 


=» AXIAL VELOCITY COMPONENT AT ROTOR INLET 
- WHEEL SPEED 

I. WORK COEFFICIENT. * » AH/(v 2 /2g o J) 

2- PRESSURE COEFFICIENT, ^ - AH ( / ( v 2 /2g Q J) 

3. FLOW COEFFICIENT. <*> = C /v 

Z 1 

k. EFFICIENCY. EFF - 

5. LOSS. XLS =» (AH - AH ( )/(v 2 /2g o J) - ^ 

Compressor stage performance parameters: definitions. 
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Figure VL2. 
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Compressor stage performance parameters: 4> for minimum loss 
and pseudo-Mach number. 
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(PSI-PSI ML > |PSI-PSI ML | 


Figure VI.2. Compressor stage performance parameters: loss and 
(loss— minimum loss), (continued) 
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